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SUMMARY 


Computerized  simulation  of  meshing  and  bearing  contact  for  spiral  bevel  gears  and  hypoid  gears 
[1,2]  is  a  significant  achievement  that  could  improve  substantially  the  technology  and  the  quality  of 
the  gears.  This  report  covers  a  new  approach  to  the  synthesis  of  face-milled  spiral  bevel  gears  and 
their  tooth  contact  analysis.  The  proposed  approach  is  based  on  the  following  ideas  proposed  in  [3] 
(i)  application  of  the  principle  of  local  synthesis  that  provides  optimal  conditions  of  meshing  and 
contact  at  the  mean  contact  point  M  and  in  the  neighborhood  of  M;  (ii)  appUcation  of  relations 
between  principle  directions  and  curvatures  for  surfaces  being  in  line  contact  or  in  point  contact. 

The  developed  local  synthesis  of  gears  provides  (i)  the  required  gear  ratio  at  M;  (ii)  a  localized 
bearing  contact  with  the  desired  direction  of  the  tangent  to  the  contact  path  on  gear  tooth  surface 
and  the  desired  length  of  the  major  axis  of  contact  ellipse  at  M;  (iii)  a  predesigned  parabolic  function 
of  a  controlled  level  (8-10  arc  seconds)  for  transmission  errors;  such  a  function  of  transmission  errors 
enables  to  absorb  linear  functions  of  transmission  errors  caused  by  misalignment  [3]  and  reduce  the 
level  of  vibrations. 

The  proposed  approach  does  not  require  either  the  tilt  of  the  head-cutter  for  the  process  of 
generation  or  modified  roll  for  the  pinion  generation.  Improved  conditions  of  meshing  and  contact 
of  the  gears  can  be  achieved  without  the  above  mentioned  parameters.  The  report  is  complemented 
with  a  computer  program  for  determination  of  basic  machine-tool  settings  and  tooth  contact  anal- 
ysi  for  the  designed  gears.  The  approach  is  illustrated  with  a  numerical  example. 

The  contents  of  the  following  sections  cover  the  following  topics; 

(1) .  Basic  ideas  of  local  synthesis  of  gears  and  the  mathematical  concept  of  this  approach 
(Chapter  1).  The  local  synthesis  discussed  in  this  chapter  is  applicable  for  all  types  of  gears  and 
provides  the  optimal  conditions  of  meshing  and  contact  at  the  mean  point  of  tangency  of  gear  tooth 
surfaces. 

(2) .  Methods  for  generation  of  the  pinion  and  the  gear  and  basic  machine-tool  settings  that  are 


X 


necessary  for  gear  generation  (Chapter  2). 

(3) .  Determination  of  geometry  of  gear  tooth  surface,  the  gear  mean  contact  point  and  .he 
principal  directions  and  curvatures  at  this  point  (Chapter  3). 

(4) .  Application  of  basic  principles  of  local  synthesis  for  spiral  bevel  gears  (Chapter  4). 

(5) .  Determination  of  pinion  machine-tool  settings  considering  as  given:(i)  the  gear  geometry, 
and  (ii)  the  conditions  of  meshing  and  contact  at  the  mean  contact  point  obtained  from  the  local 
synthesis  (Chapter  5). 

(6) .  Computerized  simulation  of  meshing  and  contact  (Tooth  Contact  Analysis)  for  spiral  bevel 
gears  that  have  been  synthesized  in  the  previous  chapters  (Chapter  6). 

(7) .  Analysis  of  the  shift  of  bearing  contact  caused  by  the  misalignment  of  gears  (Chapter  7). 

(8) .  The  theory  of  modified  roll  (variation  of  cutting  ratio  in  the  process  for  generation)  and 
mechanisms  used  for  application  of  modified  roll  (Appendix  A). 

(9) .  Description  of  developed  computer  programs  and  numerical  examples  that  illustrates  the 
application  of  those  programs. 
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1  Local  Synthesis  of  Gears  (General  Concept) 


1.1  Introduction 

The  main  goals  of  local  synthesis  are  to  provide:  (i)  contact  of  gear  tooth  surfaces  at  the  mean  point 
of  contact  of  gear  tooth  surfaces,  and  (ii)  improved  conditions  of  meshing  within  the  neighborhood 
of  the  mean  contact  point.  The  local  synthesis  is  the  first  stage  of  the  global  synthesis  with  a  goal 
to  provide  improved  conditions  of  meshing  for  the  entire  area  of  meshing.  The  criteria  of  ci  ditions 
of  meshing  are  the  transmission  errors  and  the  bearing  contact.  The  principles  of  local  synthesis 
that  are  discussed  in  this  chapter  for  face-milled  spiral  bevel  gears  can  be  appbed  for  other  types 
of  gears  as  well. 

1.2  Basic  Linear  Equations 

Consider  two  right-handed  trihedrons  Sai^f  and  S(,(Cf,eg,n)  (Fig.  1.2.1).  The  common 

origin  of  the  trihedrons  coincides  with  the  contact  point  M,  the  n-axis  represents  the  direction  of 
the  surface  unit  normal,  e/  and  are  the  unit  vectors  of  the  principal  directions  of  surface  Ei,  e, 
and  Cq  represent  the  principzd  directions  of  surface  Ej,  and  is  the  angle  formed  between  ey 
and  e,  (measured  clockwise  from  to  Cf  and  coimterclockwise  from  ey  to  e,  ).  In  reference  [4] 
three  linear  equations  were  deriv'd  that  relate  the  velocity  of  the  contact  point  over  suii  'e  Ei 
with  the  principal  curvatures  and  directions  of  contacting  surfaces  and  the  transfer  components  of 
velocities.  These  equations  are; 


(1)  ,  (1) 
a-iiv),  +  ai2Vq  -  013 

ai2t’i^’  +  =  023  ’ 

0131-1^’  +  0234'^  =  ^33  , 


(1.2.1) 


Here  (see  the  designations  in  [4]  ) 


1 


Qii 

2  ( 

-  Ky  COS 

12) 

-  K/i  sin^ 

U12 

= 

Kf  -  K/, 

021=  2 

sin ! 

2cr('2' 

Ul3 

= 

031  = 

)  _  1 

>ne„] 

022 

= 

-  Kf  sin^ 

12)  _ 

-  COS^ 

(1.2.2) 

U23 

= 

03 

(12 

2  = 

)  _ 

[ji” 

>7ie,] 

033 

K.« 

-1 

■n 

'  -  In  ^12)1 

^  j 

-ii;’) 

1 

to 

X 

i 

- 

{! 

m'sj 

{ii  X  h2) 

(r-A)} 

,(.) 

r<‘' 

= 

Equations  (1.2.1)  and  (1.2.2)  can  be  applied  for  two  cases  where:  (i)  surfaces  Li  and  ^2  are 
in  line  contact,  and  (ii)  the  surfaces  are  in  point  contact.  The  instantaneous  line  of  contact  is 
typical  for  the  case  when  the  gear  tooth  surface  (Ej)  is  generated  by  the  tool  surface  (E2).  The 
instantaneous  point  of  contact  is  typical  for  gears  with  localized  bearing  contact. 

Line  Contact 

When  the  gear  tooth  surfaces  are  in  line  contact,  the  direction  of  velocity  can  be  varied,  and 
equations  (1.2.1)  can  not  provide  a  unique  solution  for  the  unknowns  and  This  results 

in  that  the  rank  of  the  augmented  matrix 


On 

Ol2 

Ol3 

012 

O22 

023 

Ol3 

023 

O33 

must  be  less  than  2.  This  requirement  yields 
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(1.2.4) 


Equivalent  equations  are 


012033  —  013O23 


On  — 


012  - 


013O23 


(1.2.5) 
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Using  equations  (1.2.5)  and  (1.2.2)  we  obtain  equations  that  will  enable  us  to  determine 
Kj  and  K/i  for  Ej  considering  as  given  k,  and  k,  for  surface  E2’  The  equations  are: 


tan2o<'')  = 


_ 2013Q23 _ 

O23  —  flj3  +  (k,  —  Kq)a33 


(1.2.6) 


Kf  -  Kh 


2013023 
033  sin2(T('^' 


(1.2.7) 


Kf  +  Kli  —  («.,  +  Kq)  — 


2  2 
®13  "t"  ®23 


(1.2.8) 


(12)  (12)  (12) 

Equation  (1.2.6)  provides  two  solutions  ;  oj  and  cTj  =  Oi  +71/  2  and  bom  of  them  can  be 
used  for  computations  of  kj  and  k/,  that  are  represented  by  equations  (12. 7)  and  (1.2.8).  Fig. 1.2.2 
shows  the  orientation  of  two  couples  of  unit  vectors  (i  =  1,2),  with  respect  to  unit  vector 

c,.  The  magnitude  of  principal  curvature  for  the  direction  with  collinear  vectors  and  is  the 


same  although  the  notation  for  the  unit  vectors  has  been  changed.  Similarly,  we  can 

say  that  . 

Knowing  the  angle  and  the  unit  vectors  e,  and  Cq,  the  principal  directions  on  surface  Ej 

can  be  determined  with  the  following  equations, 


Cj  =  cos  CT 


(i; 


) -  (12'  ^ 
-  Sin  CT  t 
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Ch  —  sin<r*^'’e„ 


cos  (T 


(12) 


(1.2.10) 


Point  Contact 

In  the  case  of  instantaneous  point  of  contact,  the  direction  of  motion  of  the  contact  point  over 
the  surface  is  definite,  equations  (1.2.2)  for  the  unknowns  can  provide  a  unique  solution  for  the 
unknowns  t'i^*  and  Ug"  and  the  rank  of  matrix  [A]  is  2.  This  yields  that 
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Equation  (1.2.11)  yields  the  following  relation 


'‘/i- ”*12)  -  0  (1.2.12) 


Our  goal  is  to  determine  and  (the  principal  curvatures  and  directions  of  Ej)  and 

provide  at  the  mean  contact  point  (i)  a  certain  direction  of  the  tangent  to  contact  path  on  surface 
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E2  ,  (ii)  a  desired  length  of  the  major  axis  of  instantaneous  contact  ellipse,  and  (iii)  a  parabolic 
function  of  transmission  errors.  For  these  purpose  we  have  to  derive  extra  equations  in  addition  to 
equation  (1.2.12) 


Determination  of 

The  derivative  m2j(c)i)  is  the  second  derivative  of  function  02(01)  that  is  taken  at  the  mean 
contact  point;  0i  and  O2  are  the  angles  of  rotation  of  gears  1  and  2.  In  the  case  of  an  ideal  gear 
tram,  function  02(01)  is  linear  and  is  represented  by 


-V. 

<^2  =  <?1  Tr 


(1.2.13) 


However,  due  to  misalignment  between  the  meshing  gears  the  real  function  02(0])  becomes  a 
piecewise  periodic  function  wdth  the  period  equal  to  the  cycle  of  meshing  of  a  pair  of  teeth  (Fig. 
1.2.3).  Due  to  the  jump  of  angular  velocity  at  the  junction  of  cycles,  the  acceleration  approaches 
to  an  infinitely  large  value  and  this  can  cause  large  vibration  and  noise.  For  this  reason  it  is 
necessary  to  predesign  a  parabolic  function  of  transmission  error  that  can  absorb  a  linear  function 
of  transmission  error  and  reduce  the  jump  of  angular  velocity  and  acceleration  3  .  This  goal(the 
predesign  of  a  parabolic  function)  can  be  achieved  with  certain  relations  between  the  principad 
curvatures  of  contacting  surfaces  . 

Fig.  1.2.4  shows  the  predesigned  transmission  function  for  the  gear  convex  side  (Fig.  1.2.4(a)) 
and  gear  concave  side  (Fig.  1.2.3(b)).  Both  functions-c‘‘2(<i'i )  and  02  *(oi)-  are  in  tangency  at  the 
mean  contact  point  and  have  the  same  derivative  rn^i  .  at  this  po  nt. 

Consider  now  that  the  predesigned  transmission  function  is  represented  as 


O2  - 


(0! 


F{cn 


lOi 


(1.2.14) 
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Here:  and  6^^''  are  the  initial  angles  of  rotation  of  gears  1  and  2  that  provide  the  tangency 

of  gear  tooth  surfaces  at  the  mean  contact  point  M . 

Using  the  Taylor  expansion  up  to  the  members  of  second  order,  vve  obtain 
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where  m2i(Oi)  is  equal  to  A'l  / A'2  Uie  mean  contact  point  and  is  the  to  be  chosen  constant 
value:  positive  for  the  gear  concave  side,  and  negative  for  the  gear  convex  side.  The  synthesized 
gears  rotates  with  a  parabolic  function  of  transmission  errors  represented  by 


where 


^02(01)  =  ^ni2](oi  - 


(1.2.16) 
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i\ 1  1 

Equation  (1.2.16)  enables  the  determination  of  mU  considering  as  known  the  expected  values 
of  transmission  errors. 


Relation  between  Directions  of  Paths  of  Contact 

We  recall  that  velocities  tv'*  and  tv"’'  are  related  by  the  equation  4  . 


,v(  1  ' 


•  f 


(1.2.17) 


Directions  of  velocities  tv'  *  and  tv'*  coincide  with  the  tangents  to  the  contact  path  that  form 


angles  t]i  and  7^2  with  the  unit  vector  c,  (Eig  1.2.5).  Equations  (1.2.17)  vield 
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(1.2.18) 


.,(2)  ^  j,(l)  ,  j,(12) 


j,(2)^  (1)  (12) 

<]  9  ’  9 


According  to  Fig.  1.2.5 


t-'**  =  vi'han^i 


(1.2.19) 


Third  equation  of  system  (12.1)  and  equations  (1,2.18)  and  (1.2.19)  yield 


ta.nt]i  = 


n2)  ,  (12)^ . 

-0311-,,  (033  Q3lt'»  )  tan  7)2 

^33  -  «32(t'9  -  r,  tan  7/2) 


(1.2.20) 


ai3  +  a23tanT)i 


(1.2.21) 


^  (1)  _  Q33  tan 

013  +  a23tanT)i 


(1.2.22) 


Prescribing  a  certain  value  for  t}2  (choosing  the  direction  for  path  of  contact  on  Eo),  we  can 
determine  tanr/i,  t'i’*  and  I'y'*.  We  recall  that  coefficients  031,032  and  033  do  not  depend  on  the 
to-be  determined  principal  curvatures  and  k/,  and 

Relations  between  the  Magnitude  of  Major  Axis  of  Contact  Ellipse,  Its  Orientation  and 
Principal  Curvatures  and  Directions  of  Contacting  Surfaces 

Our  goal  is  to  relate  parameters  and  of  the  pinion  surface  Ei  with  the  length  of  the 

major  axis  of  the  instantaneous  contact  ellipse.  This  ellipse  is  considered  at  the  mean  contact  point 
and  the  elastic  approach  S  of  contacting  surfaces  is  considered  as  known  from  the  experimental 
data.  The  derivation  of  the  above  mentioned  relations  is  based  on  the  following  procedure 
Step  1:  Using  equations  (1.2.2),  we  obtain 
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Oil  +  022  —  ^ 

On  -  022  =  92  -  g\  cos2(T<^^) 

(ail  -  022)^  +  4o^2  =  9\-  25152  cos  2o(^2>  ^  5^ 


Step  2:  It  is  known  from  [4]  that 


0 


A  = 


25152  cos  2<t 


(1.2.23) 


(1.2.24) 


(1.2.25) 


Equation  (1.2.25)  yields 


[(oii  ^  022  4^1)^  -  (oii  -  022)^  ^  4oj2 


(1.2.26) 


Step  3:  We  may  consider  now  a  system  of  three  linear  equations  in  unknowns  011.012  and  022 


(I)  (1)  _ 

1’,  011  +  ai2  —  ai3 

ci  'an  +  i4^'a22  =  a23 


(1.2.27) 


Q]  1  -f-  0-22  —  ^  ^ 


Step  4:  The  solution  of  equation  system  (1.2,27)  for  the  unknowns  aii.ai2  and  an:  allows  to 
express  these  unknowns  in  terms  of  013,023,  A'e- I'l.'*  and  rl,’*.  Then,  using  equation  (1.2.25)  we 
can  get  the  following  equation  for  A'e 
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(1.2.28) 


Ki:  = 


AA^-{nl+nl) 

2 A  -  (ni  cos  2t)i  +  n2  sin  2ifi ) 


gi3  -  ajatan^T^i 
(1  -t  tan^  7^1)033 


ai3  tan  t?i  +  Q23)(qi3  +  Q23  tan  t?i  ) 
(1  tan^  J?i)a33 


-7 


(1.2.29) 


The  advantage  of  equation  (1.2.28)  is  that  we  are  able  to  determine  knowing  the  major  axis 
2a  of  the  contact  ellipse  and  the  elastic  approach  6. 

Step  5:  The  sought  for  principal  curvatures  and  directions  for  the  pinion  identified  with  k/.k/, 
and  can  be  determined  from  the  following  equations 


=  A'^^’  -  A'l 


(1.2.30) 


tan  2a' 


2022  _  2n2  -  A's  sin2T;i 

92  -  (an  -  a22)  92  -  2ni  +  A'j;  co5  2t7i 


T.2.31) 


2ai2  _  2n2  -  A'i;sin2qi 
sin2(T'^^'  sin2a<*^> 


(1.2.32) 


U)  . 

..  _  ..(1)  _  Si 

«/  =  -  2 


(1.2.33) 


K/i  _  — 


(1)  _  ~  9l 


(1.2.3}; 
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Step  6:  The  orientation  of  unit  vector  cy  and  T/.  is  represented  with  equations  (1.2,9)  and 
(1.2.10).  The  orientation  of  the  contact  ellipse  with  respect  to  ey  is  determined  with  angle  a*'' 
(Fig.  1.2.6)  that  is  represented  with  the  equations 


cos  2a^  ^  *  = 

sin2Q***  = 


9\  -  92  cos2tT<^^’ 

{9l  -  2^152  cos  2(7(12)  +  gl)'2 

g2  sin 

{9i  -  25132  cos  2(7(12)  ^  ^2)5 


(1.2,35) 


(1.2.36) 


The  minor  axis  of  the  contact  (26)  ellipse  is  determined  with  the  equations 


(1.2.37) 


+  yjgi  -  25152  cos  2(7  +  5? 


(1.2,38) 


Local  Synthesis  Computational  Procedure  ; 

The  following  is  an  overview  of  the  computational  procedure  that  is  to-be  used  for  the  local 
synthesis. 

The  input  data  are:  k,,  k,,  e,,  c",,,  1^*  and  6.  The  to-be  chosen  parameters  are: 

7^2,77221  and  2a.  The  output  data  are:  Ky.  k/,,  (7(12), ey  and  e/,. 

Step  1:  Choose  92  and  determine  tj^  from  equation  (1.2.20) 

Step  2:  Determine  el'*  and  r,*,'*  from  equations  (1.2.21)  and  (1.2.22) 
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Step  3:  Determine  A  from  equation  (1.2.29) 


Step  4:  Determine  A's  from  equation  (1.2.28) 

Step  5:  Determine  and  k/,  by  using  the  set  of  equations  from  (1.2.30)  to  (1.2.34) 

Step  6:  Determine  the  orientation  of  the  contact  ellipse  and  its  minor  axis  by  using  equations 
from  (1.2.35)  to  (1.2.37) 

1.3  Conclusion 

The  contact  of  tooth  surfaces  is  considered  for  two  cases;  line  contact  and  point  contact.  For  line 
contact,  the  principal  directions  and  curvatures  of  one  surface  can  be  determined  in  terms  of  the 
other’s  knowing  the  relative  motion  between  the  two  .  For  point  contact,  we  proposed  an  approach 
for  local  synthesis  of  spiral  bevel  gears  which  enables:  (i)  to  provide  a  limited  level  of  transmission 
errors,  (ii)  optimal  direction  for  the  path  of  contact  on  gear  surface  E2,and  (iii)  the  guaranteed 
length  of  the  major  axis  of  contact  ellipse. 

The  output  data  obtained  from  the  procedure  of  local  synthesis  are:  k/,  k/,.  (t* c/  and  t^,. 
The  machine-tool  settings  for  the  generation  of  the  gear  tooth  surfaces  must  be  carefully  chosen  to 
guarantee  the  above  mentioned  conditions  of  local  meshing  and  contact. 
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2  Pinion  and  Gear  Generation 


2.1  Pinion  Generation 

To  describe  the  pinion  generation  we  will  use  the  following  coordinate  system  (Fig. 2. 1.1);  (i)  5mi” 
a  fixed  coordinate  system  that  is  rigidly  connected  to  the  cutting  machine;  (ii)  5,1  -a  movable 
coordinate  system  that  is  rigidly  connected  to  the  cradle  and  performs  rotation  with  the  cradle 
about  the  Z„i\-  axis;  initially,  S^i  coincides  witu  S„n  (Fig. 2. 1.1  (b));  angle  Of  determines  the 
current  position  of  Sd  (Fig  2.1.1  (c))  :  (iii)  Coordinate  systems  5„  and  Sf,  that  are  rigidly  connected 
to  the  cradle  and  its  coordinate  system  5,:i;  systems  5„  and  Sf,  are  used  to  describe  the  installment 
of  the  head-cutter  on  the  cradle.  Angle  qi  determines  the  orientation  of  S,,  with  respect  to  Sri', 
(iv)  Coordinate  system  Sf  that  is  rigidly  connected  to  the  head-cutter  (not  shown  in  Fig.2.1.1); 
the  head-cutter  in  the  process  for  generation  performs  rotation  with  the  cradle  (transfer  motion) 
and  relative  motion  with  respect  to  the  cradle  about  an  axis  that  passes  through  0„:  (v)  Auxiliary 
coordinate  systems  5^  and  Sp  are  used  to  describe  the  installment  of  the  pinion  on  the  cutting 
machine  (Fig.2.1.1  and  Fig.2.1.2);  the  pinion  axis  forms  angle  ■)„.i  with  axis  Xd  that  is  parEillel  to 
-^mi-  (vi)  A  movable  coordinate  system  5;  that  is  rigidly  connected  to  the  being  generated  pinion; 
the  pinion  rotates  about  the  axis  Xp  and  (t>\  is  the  current  angle  of  pinion  rotation  (Fig.2.1.2). 

Henceforth,  we  have  to  differentiate  the  parameter  of  motions  that  are  performed  in  the  process 
for  generation  and  the  parameters  of  installment  of  the  head  cutter  and  the  pinion  on  the  cutting 
machine. 

In  the  process  for  generation  the  cradle  of  the  cutting  machine  with  the  mounted  head-cutter 
performs  rotation  with  angular  velocity  (Fig.2.1.2).  The  head-cutter  performs  rotational 

motion  with  respect  to  the  cradle  but  this  motion  is  not  related  with  the  process  for  generation 
and  just  provides  the  desired  velocity  of  cutting.  The  being  generated  pinion  performs  rotational 
motion  with  angular  velocity  (Fig.2.1.2)  that  is  related  with  . 

The  parameters  of  installment  of  the  head-cutter  are:  (i)  the  swivel  angle  j  (Fig.2.1.1)  and  the 
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tilt  angle  i  that  is  the  turn  angle  of  St  about  Yf,  (Fig. 2. 1.3);  5,1  =  10^0^1 1  is  the  radial  setting;  gi 
is  the  cradle  angle. 

The  parameters  of  installment  of  the  pinion  are;  £’,„i-the  shortest  machine  center  distance  (Fig. 
2.1.1,  Fig. 2. 1.2);  root  angle  sliding  base  XbiI  machine  center  to  back  Agi- 

2.2  Gear  Generation 

While  describing  the  gear  generation,  we  will  consider  the  following  coordinate  systems;  (i)  5„i2 
that  is  rigidly  connected  to  the  cutting  machine;  (ii)  Sc2  that  is  rigidly  connect  to  the  cradle,  (iii)  Sp2 
that  is  rigidly  connected  to  the  head-cutter  and  5^2i  (i'')  5-12  that  is  an  additional  fixed  coordinate 
system  rigidly  connected  to  S„l2  '■  and  (v)  52  that  is  rigidly  connected  to  the  being  generated  gear. 

The  cradle  performs  rotation  about  the  Z,„2  axis  with  angular  velocity  (Fig.2.2.1).  The 
initial  and  current  positions  of  coordinate  systems  5^2  and  5p2  with  respect  to  5^2  are  shown  in 
Fig. 2, 2.1  (a)  and  Fig.2.2.1  (b),  respectively. 

Coordinate  system  S,i2  ('t  is  rigidly  connected  to  5„,2  )  is  used  to  describe  the  installment  of 
the  gear  at  the  cutting  machine  (Fig. 2. 2, 2(a)).  In  the  general  case  apices  O2/:  and  O2  of  the  gear 
root  cone  and  pitch  cone  do  not  coincide.  Apex  O2/?  is  located  on  aous  ^^^2  of  the  cutting  machine. 
The  origin  0,^2  of  5^2  coincides  with  the  apex  O2  of  the  gear  pitch  core.  Axes  Xd2  and  Xmi  form 
angle  ■)rri2  which  is  the  gear  machine  root  angle. 

Coordinate  system  52  is  rigidly  connected  to  the  gear  that  in  the  process  of  generation  performs 
rotation  about  Xd2  with  angular  velocity  (Fig.2.2.2(b)).  Angle  4>2  is  the  current  angle  of 
rotation  of  gear  2. 

2.3  Gear  Machine  Tool  Settings 
Gear  Cutting  Ratio 

Fig. 2. 3.1  shows  the  sketch  of  the  gear  with  noncoinciding  apexes  of  the  root  and  pitch  cones. 
In  the  process  for  generation  the  pitch  line  O2P  is  the  instantaneous  axis  of  rotation.  It  is  evident 
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that  the  angular  velocity  of  rotation  in  relative  motion,  ,  must  lie  in  the  plane  that  is  formed 
by  vectors  and  (Fig.2.3.2) 


-:(p2)  _  -(2) 

W  —  U-  —  W 


(2.3.1) 


The  cutting  gear  ratio  is: 


cos  6(7  _  cos(r2  -  72) 

sinr2  sinr2 


(2.3.2) 


Gear  Settings 

Fig. 2. 3. 3  shows  the  installment  of  the  head-  cutter.  We  designate  the  mean  j'ltch  cone  distance 
O2P  (Fig. 2. 3.1,  Fig.2.3.3)  by  Am.  Then  we  obtain  (Fig. 2. 3. 3) 

Hg  =  Am  cos  8g  -  RuisinriG  (2.3.3) 


Vg  =  Rn2  cos  Xi’G  (2.3.4) 

.9,2  =  (5,2  =  0m20j,2)  (2.3.5) 

92  =  sin”’  ^  (2.3.6) 

Or2 


Here:  i'g  is  the  spiral  angle  on  the  root  cone,  R,^2  is  the  mean  radius  of  the  head  cutter.  The 


sliding  base  '0„,202  is 


E2  =  2/7  sm  ■)„,2 


(2.3.7) 
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Here:  7,^2  is  the  same  as  the  gear  root  cone  angle  72.  and  Zn  is  the  distance  bei^ween  O2/?  and 
O2,  which  are  the  apexes  of  the  root  cone  and  the  pitch  cone,  respectively. 
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3  Gear  Geometry 


3.1  Gear  Surface 

The  gear  tooth  surface  is  the  envelope  to  the  family  of  generating  surfaces.  We  recall  that  the 
cradle  carries  the  head-cutter  that  is  provided  with  finishing  blades.  The  blades  are  rotated  about 
the  axis  of  the  head-cutter  and  generate  two  cone  surfaces.  Fig.3.1.1  shows  one  of  the  cones. 

The  family  of  a  generating  surface  (the  cone  surface)  is  generated  in  S2  while  the  cradle  and 
being  generated  gear  perform  related  rotations,  about  the  Z„i2'axis  and  X2-  axis  (Fig. 2. 2. 2). 

The  derivation  of  the  gear  tooth  surface  is  based  on  the  following  procedure: 

Step  1;  W’e  represent  the  cone  surface  and  its  unit  normal  in  system  5y,2  (Fig.3.1.1)  as  follows 


= 


(Pr,  -  SG  sin  og)  cos  Og 
(fc  ~  •SG  sin  ac)  sin  9g 
—sg  cos  qg 
1 


^ 


P2 


N, 


P2t 


,V  -  ^^P^  V  ^P^ 

09  G  OsG 


(3.1.1) 


(3.1.2) 


i.e.. 


np2  = 


-  cos  QG  COS  9g 

-  cos  Qg  sin  9g 

sin  QG 


(3.1.3) 


Here:  sg  and  9g  are  the  surface  coordinates;  qg  is  the  blade  angle;  is  the  radius  of  the 
head-cutter  that  is  measured  at  the  bottom  of  the  blades.  It  is  evident  (Fig. 3. 1.2)  that 


Pc  =  Ru2  ± 


PM 

'Y 


(3.1.4) 
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Here:  iZu2  is  the  nominal  radius,  PW  is  the  so  called  point  width;  the  positive  sign  in  (3.1.4) 
corresponds  to  the  gear  concave  side  and  the  negative  sign  corresponds  to  the  gear  convex  side. 

Equations  (3.1.1)  and  (3.1.3)  represent  both  generating  cones  with  ac  >  0  for  the  gear  convex 
side  and  ag  <  0  for  the  gear  concave  side. 

Step  2;  The  family  of  generating  surfaces  that  is  generated  in  .^2  is  represented  by  the  following 
matrix  equation 

^2  (  ][-^^T7t2C2  l['^fr^2P2  ]^p2  (3.1.5) 

Here  (Fig. 2. 2. 2,  Fig.2.2.1); 

■  1  0  0  0  ■ 

0  cos  <t>2  sin  (t>2  0 

0  -  sin  <t>2  cos  02  0 

0  0  0  1 

cos'>„,2  0  sin-)„i2  -A'B2sin')m2 

0  10  0 
-sin-),,, 2  0  cosy„i2  -A'b2COS7„2 
0  0  0  1 


(3.1.6) 

(3.1.7) 


17 


[^C2P3  ]  — 


10  0  Sr2  COS  q2 

0  10  Sr2  sin  q2 

0  0  1  0 

0  0  0  1 


(3.1.8) 


cos  4>p 

-  sin  0p 

0 

0 

sin  4>j, 

cos  dp 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

(3.1.9) 


The  machine  root  angle  '),n2  in  equation  (3.17)  is  equal  to  gear  root  cone  angle  ■) . 
Step  3:  The  derivation  of  the  equation  of  meshing  is  based  on  the  equation 


^m2  ■  ^m2  ^ 


(3.1.10) 


The  subscript  “m2”  means  that  vectors  in  equation  (3.1.10)  are  represented  in  coordinate  system 
Srri2  ;  nr/iS  is  the  unit  normal  to  the  generating  surface;  is  the  relative  (sliding) 

velocity  Vector  n„,2  is  represented  by  the  matrix  equation 


-  cosaGCOs(0G  +  0J.) 

rim2  =  [Lm2P2]np2  =  -  cos  OG  sin(^G  +  0p) 


(3.1.11) 


sin  og 


where  3x3  submatrix  [Mrnjn]- 

We  consider  that  the  axes  of  rotation  of  the  cradle  and  the  gear  intersect  each  other  (Fig. 2. 2. 2(a)), 
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(3.1.12) 


,-T<p2) 

‘■m2 


.  m2 


ni2 


5(;'2 ) 
^rn2 


X  r. 


where 

=  ;- COS72  0  - (3.1.13) 

^aG 

We  assumed  that 
following  relation 


Here 


1  in  equation  (3.1.13).  Equations  from  (3.1.10)  to  (3.1.13)  yield  the 


SG 


M^G-<t>j,) 

B{OG,0r) 


(3.1.14) 


A{6g.  Oj,) 


n„,jx:--4i(sin'V2 
+  nrn2zAi  cos  72 


^— )'  -  cos  72  *  .42(sin72 

BaG 


Y-y- 

B-aG 

(3.1.15) 


B{6G,<i>p) 


-n„i2x  sin  OG  sin(0G  -  <?;.)  +  n,„j„isin72  -  — — )sinaGCos(0G  Op) 

^aG 

-  cos  og  cos  72,  cos 72  sin og  sin(^G  O,,)  (3.1.16) 


.4i  =  r,.  sin(%  o,,)  -  .fr2  sin(^2  -  O,,) 


(3.1.17) 


.42  =  rr,  cos(0G  *  0p)  -  Sr2  cos((72  -  Op) 


(3.1.18) 


Step  4:  Equations  (3.1.5)  and  (3.1.14)  considered  simultaneously  represent  the  gear  surface  in 


three-  parametric  form  but  with  related  parameters.  Since  parameter  sq  in  equation  of  meshing 
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(3.1.14)  is  linear,  it  can  be  eliminated  in  equation  (3.1.5),  and  then  the  gear  tooth  surface  will  be 
represented  in  two-parametric  form,  by  the  vector  function  ^2(^0-©?;)  • 

3.2  Mean  Contact  Point  and  Gear  Principal  Directions  and  Curvatures 

The  mean  contact  point  M  is  shown  in  Fig.2.3.1.  Usually,  M  is  chosen  in  the  middle  of  the  tooth 
surface.  The  gear  tooth  surface  and  the  pinion  tooth  surface  must  contact  each  other  at  M . 

The  procedure  of  local  synthesis  discussed  in  section  2.1  is  directed  at  providing  improved 
conditions  of  meshing  and  contact  at  M  and  in  the  neighborhood  of  .1/.  The  location  of  point 
M  is  determined  with  parameters  XL  and  RL  (Fig.2.3.1)  that  are  represented  by  the  following 
equations 

h.  c 

XL  -  .4,,,  cos r2  -  (/>G - - — )sinF2  (3.2.1) 

RL  =  .4„.sinr2  -  (fcc  -  ^^^^^)cosF2  (3.2.2) 


Here:  -4,n  is  the  pitch  cone  mean  distance;  /i,„  is  the  mean  whole  depth;  he  is  the  gear  mean 
dedendum:  c  is  the  clearance  Equations  (3.2.1).  (3.2.2)  and  vector  equation  f2{0c.Oj.)  for  the  gear 
tooth  surface  allows  to  determine  the  surface  parauneters  and  o*  for  the  mean  contact  point 
from  the  equations 


-V2(^?G-«-\'.)  =  -VI 


(3.2.3) 


z^ob-o;.)  =  {Rif 


(3.2.4) 
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Gear  Principal  Directions  and  Curvatures 


The  gear  principal  directions  and  curvatures  can  be  expressed  in  terms  of  principal  curvatures 
and  directions  of  the  generating  surface  (see  chapter  (13)  in  [4]),  that  is  the  cone  surface. 

Step  1:  The  cone  principcil  directions  are  represented  in  Sp2  by  the  equations  (see  (3.1.1)) 


df. 


P2 


Mp) 

*P2 


I 

'dffc 


=  [-sin0G  cos^G  0]^ 


df 


dp)  — 
qp2 


dsg  _ 
dfn 
dsg 


=  [- sin  og  COS  ^G  -sinoGsin^G  -cosog]^ 


(3.2.5) 


(3.2.6) 


The  superscript  “p”  indicates  that  the  cone  surface  Ep  is  considered.  Unit  vector  is  directed 
along  the  cone  generatrix  and  unit  vector  etpl  is  perpendicular  to  The  unit  vectors  of  cone 
principal  directions  are  represented  in  5^2  by  the  equations 


=  [-  sin(6lG  +  <t>p)  cos(flG  +  <t>p)  0]^ 


(3.2.7) 


=  [- sin  ag  cos{6g  +  <t>p)  ~  sin  qg  sin(^G  +  c!>p)  -  cosqg’^  (3.2.8) 


The  cone  principal  curvatirres  are; 


^(P)  ^ 

Tr  -Sg  sin  QG 


and  =  0 


(3.2.9) 


Step  2:  The  determination  of  principal  curvatures  and  directions  for  gear  tooth  surface  E2  is 
based  on  equations  from  (1.2.6)  to  (1.2.8).  The  superscript  “2”  in  these  equations  must  be  changed 
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for  “p”  and  superscript  “1”  for  “2”.  The  second  derivative  of  cutting  ratio,  s  m'^2 

because  the  cutting  ratio  is  constant.  The  principal  curvatures  of  the  gear  tooth  surface  will  be 
determined  as  /cy  and  K^^.  The  principcil  directions  on  gear  tooth  surface  will  be  represented  in  by 
ej  and  and  they  can  be  determined  from  equations  (1.2.9)  and  (1.2.10).  To  represent  in  S2  the 
principal  directions  on  gear  tooth  surface  T,2  and  its  unit  normal  we  use  the  matrix  equation  that 
describe  the  coordinate  transformation  from  5,„2  to  S2-  This  equation  is 

^2  “  ,^2d2  (3.2.10) 

Here:  a,„2  stands  for  vectors  n„i2-  ^hrn2:  ^2  stands  for  n2,  and 
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4  Local  Synthesis  of  Spiral  Bevel  Gears 


4.1  Conditions  of  Synthesis 

The  basic  principles  of  local  synthesis  of  gear  tooth  surfaces  discussed  in  Section  1  will  enable  us 
to  determine  the  principle  curvatures  and  directions  of  the  being  synthesized  pinion.  Thus,  we  will 
be  able  to  determine  the  required  machine-tool  settings  for  the  pinion.  While  solving  the  problem 
of  local  synthesis,  we  will  consider  as  known: 

(i)  The  location  of  the  mean  contact  point  A1  in  a  fbced  coordinate  system,  and  the  orientation 
of  the  normal  to  gear  surface  Lj- 

(ii)  The  principle  curvatures  and  directions  on  L2  at  M.  The  local  synthesis  of  gear  tooth 
surfaces  must  satisfy  the  following  requirements: 

(1)  The  pinion  and  gear  tooth  surfaces  must  be  in  contact  at  M . 

(2)  The  tangent  to  the  contact  path  on  the  gear  tooth  surface  must  be  of  the  prescribed  direction. 

(3)  Function  of  gear  ratio  in  the  neighborhood  of  mean  contact  point  must  be  a  linear 

one,  be  of  prescribed  value  at  M  and  have  the  prescribed  value  for  the  derivative  m2j(d>i)  at  A/. 
The  satisfaction  of  these  requirements  provides  a  parabolic  type  of  function  for  transmission  errors 
of  the  desired  value  at  each  cycle  of  meshing. 

(4)  The  major  axis  of  the  instantaneous  contact  ellipse  must  be  of  the  desired  value  (with  the 
given  elastic  approach  of  tooth  surfaces). 

4.2  Procedure  of  Synthesis 

We  will  consider  in  this  section  the  following  steps  of  the  computationEil  procedure;  (i)  representa¬ 
tion  of  gear  mean  contact  point  in  a  fixed  coordinate  system  5/,  ;  (ii)  satisfication  of  equation  of 
meshing  of  the  pinion  and  gear  at  the  mean  contact  point  ;  (iii)  representation  of  principle  directions 
on  gear  tooth  surface  E2  in  5/,;  (iv)  observation  of  the  desired  derivative  T/i2](<hi).  (v)  observation 
at  the  mean  contact  point  of  the  desired  direction  of  the  tangent  to  the  path  contact  on  gear  tooth 
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surface  ;  (vi)  observation  at  the  mean  contact  point  of  the  desired  length  of  the  major  axis  of  the 
contact  ellipse;  (vi)  determination  of  principal  directions  and  curvatures  on  pinion  tooth  surface 
El  at  the  meam  contact  point. 

Step  1:  We  set  up  a  fixed  coordinate  system  S/,  that  is  rigidly  connected  to  the  gear  mesh 
housing  (Fig.4.2.1(a)).  In  addition  to  S/i,  we  will  use  coordinate  systems  52  (Fig. 4. 2. 1(a))  and  5i 
(Fig.4. 2.1(b))  that  are  rigidly  connected  to  gears  2  and  1,  respectively  We  designate  with  c!>2  and 
<t)j  the  angles  of  rotation  of  gears  being  in  mesh.  We  have  to  emphasize  that  with  this  designation 
(p^{i  =  1,2)  we  differentiate  the  angle  of  gear  rotation  in  meshing  from  the  angle  cs,  of  gear  rotation 
in  the  process  of  generation. 

The  orientation  of  coordinate  system  5/,  is  based  on  following  considerations:  (i)  The  axes  of 
rotation  of  the  pinion  and  the  gear  in  a  drive  of  spiral  bevel  gears  intersect  each  other.  Taking  into 
account  the  possible  gear  misaligmnent,  we  will  consider  that  the  pinion-gear  axes  are  crossed  at 
angle  F  and  the  shortest  distance  is  E.  (ii)  We  will  choose  that  Xh  coincides  with  the  pinion  axis 
and  Oh  is  located  on  the  shortest  distance  (Fig. 4. 2. 1(a)).  (iii)  Considering  as  given  the  shaft  angle 
r,  we  will  define  j/,-  the  unit  vector  of  V/,  -  as  follows 


jh  = 


ih  X  dh 

ii/,  X  dh\ 


(4.2.1) 


where  a/,  is  the  unit  vector  of  gear  axis  that  is  parallel  to  plane  {Xh^Yh)- 
The  coordinate  treinsformation  from  S2  to  5/,  is  based  on  matrix  equation 

(4.2.2) 


where  5j  (Fig. 4. 2.1)  is  an  auxiliary  fixed  coordinate  system. 


The  unit  normal  to  E2  is  repre- 
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sented  in  S/,  as 


Here  (Fig.  4.2.1) 


=  [Fh<f][Fd2]"2(^G,  0p) 


{^'^<12]  - 


10  0  0 

0  -  cos  02  <?2  0 

0  -  sin  4>2  ~  ® 

0  0  0  1 


Muu]  = 


cos  r  0  sin  r  0 


1  0 


-  sin  r  0  cos  r  0 


0  0  0  1 


(4.2.3) 


(4.2.4) 


(4.2.5) 


where  F  is  the  shaft  angle. 

Equations  (4.2.3),  (4.2.2)  and  (4.2.3)  enable  to  represent  in  5/,  the  position  vector  and  unit 
contact  normal  at  M  by 


h 


(4.2.6) 


where  (0q,  0*)  are  the  surface  coordinates  for  the  mean  contact  point  at  112  :  the  angle  01,  of  rotation 
of  gear  2  will  be  determined  from  the  equation  of  meshing  (see  below). 
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Step  2:  The  equation  of  meshing  of  pinion  and  gear  at  the  mean  contact  point  is 


(4.2.7) 


Here  (Fig.  4.2.1) 


4^^’  =  [(4^»-4^’)xrr]-(Fx.<^>) 

=  [-1  0  of  (i4"i  =  i) 

^  ®  -sinF]^ 

since  at  point  M  the  angular  velocity  ratio  is 

u;<^>  _  Ny 

~  Y2 


(4.2.8) 


(4.2.9) 


(4.2.10) 


(4.2.11) 


Substituting  equations  (4.2.3),  (4.2.8)-  (4.2.11)  in  equation  (4.2.7),  we  can  solve  equation  (4.2.7) 
for  <^2-  Usually  equation  (4.2.7)  yields  two  solutions  for  <^>'2  but  the  smaller  one,  say  (<^>2)*,  should 
be  chosen. 

Step  3:  We  consider  as  known  the  principal  curvatures  and  directions  on  Ej  at  any  point  of  E2, 
including  the  mean  contact  point  (see  section  3).  To  represent  in  S;.  the  principal  directions  at  the 
mean  contact  point,  we  use  the  matrix  equation 

4^>  =  \L,J\Ld2Vi2  (4.2,12) 
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where  0,2  is  the  unit  vector  of  principal  directions  on  E2  that  is  represented  in  82-  The  following 
steps  of  computationtil  procedure  are  exactly  the  same  that  have  been  described  in  section  1.2. 
This  procedure  permit  determination  of  the  pinion  principal  directions  and  curvatures  at  the  mean 
contact  point. 
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5  Pinion  Machine-Tool  Settings 


5.1  Introciuci,ion 


We  consider  at  this  stage  of  investigation  as  known: 

(i)  the  common  position  vector  and  imit  normal  nj*’  at  the  point  of  contact  point  M  of  E2 
and  El 

(ii)  pinion  surface  principal  directions  and  curvatures  at  M . 

The  goal  is  to  determine  the  settings  of  the  pinion  and  the  head-  cutter  that  will  satisfy  the 
conditions  of  local  synthesis.  We  consider  that  the  pinion  surface  and  the  generating  surface  are 
in  hne  contact.  Henceforth,  we  will  consider  two  types  of  the  generating  surface:  (a)  a  cone 
surface,  and  (b)  a  surface  of  revolution.  We  consider  that  each  side  of  the  pinion  tooth  is  generated 
separately  and  two  head-cutters  must  be  apphed  for  the  pinion  generation. 


5.2  Head-Cutter  Surface 
Cone  Surface 

The  cone  surface  is  generated  by  straight  blades  being  rotated  about  the  rp-axis  (Fig.  5.2.1(a)). 
The  Ef  equations  are  represented  in  roordinate  system  5/r  that  is  rigidl  •  connected  to  the  head- 
cutter  as  following: 


ff 


+  Sf  sin  Qp)  cos 
{Rep  +  sin  Qf  ^  sin  Op 

-Sf  cos  Qf 

1  J 


(5.2.1) 


Here:  sp  and  Op  are  the  surface  coordinates;  ap  and  R^p  are  the  blade  angle  and  the  radius  of  the 
cone  in  plane  zp  =  0.  The  blade  angle  Qf  is  standardized  and  is  considered  as  known.  Partimeter 
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Q/r  is  considered  as  negative  for  the  pinion  convex  side  and  af  is  positive  for  the  pinion  concave 
side.  The  point  radius  is  considered  as  unknown  and  must  be  determined  later. 

The  unit  normal  to  pinion  tooth  surface  is  represented  as 


77F  = 


and 


Off 

do F  ds F 


(5.2.2) 


i.e., 


fir  -  -icos  Of  cos^^f  cosQ/rsint?f  sino/r'^  (5.2.3) 


The  principal  directions  on  the  cone  surface  are: 


(5.2.4) 


(5.2.5) 


The  corresponding  principal  curvatures  are 


■4F) 


dvF 

d*  F 
dfF 
d.*!  F 


dfF 
SOf  _ 

dfp 

dOF 


—  I  -  sin  Of  cos  Of 


vT 


=  [sin  Of  cos  sina/rsint?p 


cos  OF 


JF)  COS  Of  1 

K j  —  - - — ~ —  3nd  K  j  — 

^  sin  a/r  p 


(5.2.6) 


Surface  of  Revolution 
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We  consider  that  the  head-cutter  surface  D/r  is  generated  by  a  circular  arc  of  radius  p  by 
rotation  about  the  ::o-axis  that  coincides  with  the  zjr-axis  of  the  head-cutter  (Fig.  5.2.1(b))  and 
(Fig. 5. 2. 1(c)).  The  shape  of  the  blade  is  represented  in  So  by  the  vector  equations 


OoN  =  OoC  +  CN  ^  (A;(">  +  PCOS  X)io  +  (4'^>  +  psin  X)k„ 


(5.2.7) 


Here:  (Xq'^\  Zq'^^)  are  algebraic  values  that  represent  in  S„  the  location  of  center  C  of  the  arc; 
p  =  |CA^j  is  the  radius  of  the  circular  arc  and  is  an  algebraic  value,  p  is  positive  when  center  C 
is  on  the  positive  side  of  the  unit  normal.  ;  X  is  the  independent  variable  that  determines  the 
location  of  the  current  point  N  of  the  arc.  By  using  the  coordinate  transformation  from  5„  to  5f 
(Fig. 5. 2. 1(c)),  we  obtain  the  following  equations  of  the  surface  of  the  head-cutter: 


r> 


-I-  p  cos  A)  cos  9f 
( +  p  cos  A)  sin  dp 
z\f^  -I-  p  sin  A 
1 


(5.2.8) 


where  A  and  9y  are  the  surface  coordinates  (independent  variables). 

The  surface  unit  normal  is  represented  by  the  following  equations 


and 


(5.2.9) 


Then  we  obtain 
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(5.2.10) 


nf  =  -[cos  Acosflf-  cos  A  sin  sinA]^ 


The  variable  A  at  the  mean  contact  point  M  has  tlie  same  value  as  the  standardized  blade  angle 
Qf .  The  principal  directions  on  the  head-cutter  surface  are 


^  =  [-sin^f  cos  Of  0]^ 

orir 


(5.2.11) 


=  [sin  A  cos  Of  sinAsin^/r  -  cos  A 

!  — I 
'  d\  ' 


(5.2.12) 


The  principal  curvatures  are 


(/T)  cos  A 

~  .r{c)  .  T 

+  p cos  A 


,  (F)  t 

and  ~  ~ 


(5.2.13) 


The  radius  R„,  of  the  head-cutter  in  plane  (Fig.  5.2.1)  can  be  determined  from  the  equations 


(5.2.14) 


5.3  Observation  of  a  Common  Normal  at  the  Mean  Contact  Point  for  Surfaces  Ep, 


E2  ,  Sf  and  S] 


We  consider  that  at  the  mean  contact  point  M  four  surfaces-  Sp,S2,5^F  and  Si-  must  be  in 
tangency.  The  contact  of  Sp  and  S2  at  M  has  been  already  provided  due  to  the  satisfication  of 
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their  equation  of  meshing  (3.1.10).  Our  goal  is  to  determine  the  conditions  for  the  coincidence 
at  M  of  the  unit  normals  to  Ef,  Ep  and  E2.  The  tangency  of  Ej  with  the  three  above  mentioned 
surfaces  will  be  discussed  below. 

We  will  consider  the  coincidence  of  the  unit  normals  in  coordinate  system  5„,i .  To  determine  the 
orientation  of  coordinate  system  Si,  with  respect  to  5^1,  let  us  imagine  that  the  set  of  coordinate 
systems  Si,,  5i  and  S2  (Fig. 4. 2.1)  with  gears  1  and  2  is  installed  in  5,„i  with  observation  of  following 
conditions  (Fig. 5. 2. 2):  (i)  axis  xi,  of  Si,  coincides  with  axis  Xj,  of  5p;  (ii)  coordinate  system  5i 
coincides  with  5p  and  the  orientation  of  Si,  with  respect  to  Si  is  designated  with  angle  4>i,  =  (d>i)° 
where  </>/,  is  the  to  be  determined  instalment  angle.  Angle  <f>i,  will  be  determined  from  the  conditions 
of  coincidence  of  the  unit  normals  to  Ef',E2  Ep  and  Ei.  The  procedure  for  derivation  is  as  follows: 

Step  1:  Consider  that  the  coordinate  system  S/,  with  the  point  of  tangency  of  surfaces  E2  and 
Ep  is  installed  in  S„,i.  W'e  may  represent  the  surface  unit  normal  in  Smi  by  using  the  following 
matrix  equation  (Fig.  5.2.2). 


■  COS71 

0 

-sin7i 

■  1 

0 

0 

0 

1 

0 

0 

cos  01, 

-  sin  (t>ii 

sin7i 

0 

COS71 

0 

sin  <pi. 

cos  (t>l. 

. 

(5.3.1) 


—( 2 ) 

The  unit  vector  has  been  represented  by  equation  (4.2.3). 

Step  2:  The  unit  vector  to  the  surface  of  Ef  of  the  head-cutter  that  generates  the  pinion 
has  been  represented  in  Sf  by  equation  (5.2.3)  for  a  cone  and  equation  (5.2.10)  for  a  surface  of 
revolution.  Axes  of  coordinate  systems  5f  and  .S',,,!  have  the  same  orientation  and 


-if) 

«ml  =  n 


(f) 


(5.3.2) 
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Equations  (5.3.1),  (5.3.2),  (5.2.3)  and  (5.2.10)  yield  the  following  equations 


cos  6p  —  ~ 


(2)  ,  ■ 

‘^Xh  +  5in  7i 

cos  7i  cos  a  F 


(5.3.3) 


cos  0;,  = 


(-1,1’)'  +  ("'IV 


sin^/,  = 


(nil’)^  +  ("il’)^ 


(5.3.4) 


Here: 


Ui  =  -  cos  Of  sin  ^2  =  cosQ;fsin7i  cosff/r  ~  sin  of  cos 


(5.3.5) 


The  advantage  of  the  proposed  approach  is  that  the  coincidence  of  the  unit  normals  to  surfaces 
Sfi  £21  Sp  and  Ej  can  be  achieved  with  standard  blade  angles  and  without  a  tilt  of  the  head-cutter. 


5.4  Basic  Equations  for  Determination  of  Pinion  Machine-Tool  Settings 

At  this  stage  of  investigation  we  will  consider  as  known:  ’  "mi  ^i- 

necessary  to  determine:  iZcp.£mi,-^BC,  and  m^-..  Here:  ^ and 

are  the  principal  curvatures  and  unit  vectors  of  principal  directions  on  the  pinion  surface  that  are 

taken  at  meein  contact  point  M;  and  are  the  position  vector  of  M  and  the  contact  normeil 

(  F\ 

at  M.  The  subscript  “ml”  indicates  that  the  vectors  are  represented  in  Smi.  Designations  k) 
and  indicate  the  principal  curvatures  of  the  surface  of  the  pinion  head-cutter  that  are  taken 
at  M.  The  angle  is  formed  by  the  unit  vectors  and  of  principal  directions  on  Ej  and 
Ef,  Rep  is  the  cutter  ’’point  radius”  (Fig.5.5.1)  that  is  measured  in  plane  rf  =  0  and  is  dependent 
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on  K^j^\  Ern\  and  Xci  are  the  pinion  settings  for  its  generation  (Fig. 2. 1.1  and  Fig. 2. 1.2);  mpi 


which  is  equal  to  and  are  the  cutting  ratio  and  its  derivative. 

We  recall  that  the  pinion  surface  curvatures  and  k'jj  have  been  determined  in  the  process 
of  local  synthesis.  Vectors  i  *  have  been  determined  in  system  S^.  To  represent  these 

vectors  in  5t„i  we  have  to  apply  the  coordinate  transformation  from  5/,  to  5„ii  similar  to  equation 
(5.3.1). 


^Tn\  —  [Tnijp]  [Tp/i]u/j 


(5.4.1) 


where  represents  that  principal  directions  of  the  pinion  surface  and  the  position  vector 
of  mean  contact  point  Umi  represents  the  corresponding  vectors  and 

Now  our  goal,  as  it  was  mentioned  above,  is  to  determine  k^/  Emi^XcurriFi  and 

m^j.  We  recall  that  vectors  and  e'^  are  known  from  the  local  synthesis,  and  and 
become  known  from  equations  (5.2.4)  and  (5.2.5)  for  straight  blade,  and  from  equations  (5.2.11) 
and  (5.2.12)  for  curved  blade,  after  the  coincidence  of  the  contact  normal  to  surfaces  E2,  Ep  and 
Ef  is  provided.  Thus  parameter  can  be  determined  from  the  equations 


sin  (rdf )  =  X  ] 


cos^dF)  _ 

cos  a  - 


(5.4.2) 


According  to  Fig.  5.5.1,  since  the  Z„,i-axis  is  parallel  to  Z/r-axis,  surface  parameter  Sf  for  the 
cone  surface  at  mean  point  can  be  determined  as: 


Z 


ml 


COS  Qf 


(5.4.3) 
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/  p<  \ 

Parameter  Kjj  is  equal  to  zero  for  a  cone  surface  of  the  head-cutter  and  it  must  be  chosen 
for  a  head-cutter  with  a  surface  of  revolution.  Then,  the  number  of  remauning  pareuneters  to-be 
determined  becomes  equal  to  five  and  they  are;  Kj  ,  £mi,  and 

It  will  be  shown  below  that  we  can  derive  only  four  equations  for  determination  of  the  unknowns 
of  the  output  data.  Therefore  one  more  parameter  has  to  be  chosen,  and  this  is  m^j-the  modified 
roll.  Usually,  it  is  sufficient  to  choose  =  0,  but  the  more  general  case  with  0  is 

considered  in  this  report  as  well. 

The  to  be  derived  equations  are  as  follows. 


_(A/)  ^IF)  „ 


^Jn<^22  —  On 


an(l23  =  012013 


O12O33  =  013023 


(5.4.4) 


(5.4.5) 


(5.4.6) 


(5.4.7) 


Equation  (5.4.4)  is  the  equation  of  meshing  of  the  pinion  and  head-cutter  that  is  applied  at 
the  mean  contact  point.  Equations  from  (5.4.5)  to  (5.4.7)  come  from  the  conditions  of  existence  of 
instantaneous  line  contact  between  Ei  and  Ef .  The  coefficients  in  equation  (5.4.5)-(5.4.7)  are 
represented  as  follows, 


On  = 


{F} 


» M  2 
-  K  f  COS 


-  K 


(II 

// 


(5.4.8) 
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ai2  =  ^21  =  — — ^ — ^sin2£r(^^^ 


„  _  ^  _  ^(^)„(lf')  \-;{ini:JF) 

®13  —  ^31  —  ^  I  —  U}  n  Cj 


{^)  (1)  •  2  (IF)  (>)  2  (IF) 

“22  =  K/;  -  K/  Sm  a'  ’  -  K}/ COS  <t'  ' 


(F)  (IFI 

^23  =  032  =  -«//  I’// 


a33  = 

n„.  wl'l  •i'-']  +  .f'j  +  („-.  .t^') 


Vectors  in  c^iaation  of  meshing  (5.4.4)  can  be  represented  as  follows 


4ni  =  (cos7i  0  sin7i]^ 


(iuJ">l  =  1) 


Sf) 

«^rr.l 


;^!0  0  ij 

^ap 


where  Rap,  which  is  equal  to  - ,  is  the  ratio  of  roll. 

TTlFl 


-.(IF) 


-(1) 


“'ml 


‘ml 


(5.4.9) 

(5.4.10) 

(5.4.11) 

(5.4.12) 

(5.4.13) 

(5.4.14) 

(5.4.15) 

(5.4.16) 

(5.4.17) 
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41)  _  -.(1)  ^  JM) 

^tr  ~  “'rnl  ^  ^ml 


-  Vmi  sin  7l 
■^rnl  sin'll  ^ml 
^  ml  COS  'yi 


(5.4.18) 


’  =  -Z'  X  (■!"'  -  0707)  = 


^  ml^Fl  "7  £'mimf'l 

+  mfiA'ci  cos7i 
0 


(5.4’P) 


5.5  Determination  of  Cutter  Point  Radius 


Step.l;  Equations  (5.4.3),  (5.4.6)  and  (5.4.7)  yield  the  following  expression  for  Kj 


Jf)  K;^*(Ky  *  COS^  +  K;’,'  sin‘ 

^  -  Kj^  *  sin^  -  K^y  cos^ 


(5.5.1) 


Step  2:  According  to  Meusnier’s  theorem,  the  cutter  radius  Rm  at  the  mean  contact  point  is 
(Fig.5.5.1) 


Rm  - 


cos  a/r 


(5.5.2) 


As  shown  in  Fig.5.2.1.  the  cutter  point  radius  can  be  determined  for  a  straight  blade  cutter  as 


follows. 


Rirp  -  Rrn  “  -V  ^'n  Of 


(5.5.3) 


For  the  arc  blade,  the  location  of  the  center  of  the  arc  can  be  determined  in  5„  by  following 


equations. 
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A",*,'  '  =  R,„  -  p  cos  a  f 


(5.5.4) 


=  Rin  -  p  sin  ft  f 


(5.5.5) 


Knowing  Xo'^  and  Zlf'\  we  can  determine  the  point  radius  for  the  arc  blade  by  equation  (5.2.14). 
In  order  to  find  the  position  vector  of  the  center  of  the  head  cutter,  we  define  the  following  two 
vectors  in  as  shown  in  Fig.  5.5.1. 


xi'>)  _ 


p  =  n  cos  Of  -  Cjj  sin  ft  f 


(5.5.6) 


cos  df 

-  sin  df 

0 

0  ' 

■  A-r’  ■ 

A',','^*  cos  df 

sin  df 

cos  df 

0 

0 

0 

A',’,' '  sin  df 

0 

0 

1 

0 

z!''' 

0 

0 

0 

1 

1 

1 

(5.5.7) 


where,  is  a  unit  vector  directed  from  the  blade  tip  to  the  cutter  center  Of,  and  p^^'^  is 
a  position  vector  directed  from  Of  to  the  arc  center  C .  Referring  to  Fig. 5. 2. 2  and  Fig. 5. 5. 2,  the 
position  vector  of  the  cutter  center  Of  with  respect  to  O/,,  can  be  determined  in  system  5,„i 

as  follows, 

For  straight  blade: 


AM) 

^rn  1 


af^ 

^F*  II 


(5.5.8) 
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For  arc  blade  . 


r. 


AM) 

^rnl 


pnmi  ~ 


(5.5.9) 


It  can  be  verified  that  the  Z„a  component  of  is  zero,  since  equations  (5.4.3)  and  (5.5.5) 
are  observed.  It  is  worth  to  mention  that  6p  and  sjr  are  the  surface  coordinates  where  the  contact 
is  at  mean  point.  The  values  of  and  sjp  will  serve  as  the  initial  guess  in  tooth  contact  analysis. 

5.6  Determination  of  mfi  =  75^-.  E\fi  and  Xgi 

i  tap 

The  determination  of  cutting  ratio  i2„p,  settings  Emi  and  Xqi  is  based  on  application  of  equations 
(5.4.2),  (5.4.4)  and  (5.4.5). 

Initial  Derivations 

It  is  obvious  that  equation  of  meshing  (5.4.4)  is  satisfied  at  point  M  if  the  relative  velocity 
lies  in  plane  that  is  tangent  to  the  contacting  surfaces  at  M.  Thus,  if  velocity  satisfies  the 

equation. 


f(F\)  _  JF\)aF)  ^  (FDAF) 
I  _  ij  Cj  -r  Vjj  ejj 


(5.6.1) 


it  means  that  equation  of  meshing  (5.4.4)  is  also  satisfied.  Assuming  that  vectors  of  equation 
(5.6.1))  are  represented  in  coordinate  system  S,„i,  we  obtain 
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">1  ^ImlY  +  ^7/  ^IIm\Y 


(5.6.2) 


(FI)  (F) 

^7  ^ImlZ  +  ^//  ^llmlZ 


For  further  derivations  we  will  use  the  foUow'ing  expressions  for  ai3  and  033. 


Ul3  —  ^  +  M\\t\  +  A/12 

<123  —  '^7/^1'//  ^  A-  M2\t\  +  A/22 


(5.6.3) 


Here, 


Kf  (P)  (P) 

A/ll  —  ^mlX^ImlY 

A/12  =  -  C0S7i[n^iye^^\^  -  -nmi^e/^’iy] 

(f)  (P) 

A/21  -  llmlAe//„iy  - 


1  *  )  (  /'  ) 

A/22  =  -cos7i[n,„iyejj^i^  - 


(5.6.4) 


(1  =  rrifj  -  sin7i 


(5.6.5) 


Using  equations  (5.6.2),  (5.6.3)  and  (5.4.16),  we  obtain 


(Fl)  (F'  ,  (Fl)  (F)  ,  ,,  „ 

^ImlZ  +  ^I/rnlZ  +  ^  7l  =  0 


(5.6.6) 


Following  Derivations 
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/  JT!  \  ^  J  ) 

Step  1:  Expressions  for  Vj  and  Vjj  .  Equations  (5.6.6)  and  (5.4.4)  represent  a  system  of 


two  linear  equation  in  unknowns  v)  and  Vjj  .  The  solution  of  these  equations  for  the  unknowns 
yields; 


(Fl) 


—  1/21^1  +  L 


27 


(Fl) 


—  Liiti  +  L 


12 


(5.6.7) 


(5.6.8) 


Here, 


^21  -  (f)  (F)  ,  (F)  (F) 

ai2«/  ^llmlZ  +  ®11«//  ^ImlZ 


L22  = 


_  gf/,!iiz(Qn-^^22  -  o-nMu)  -  aiiK/;’Vmi  cos7i 


„  ,.(F)  (F)  JF)AF) 


(5.6.9) 


(5.6.10) 


ill  - 


il2  — 


~Ql2'^^n) 


ai2^/  ^IlmlZ  +  “11''//  ^hn\Z 

-  0'12Mi2)  -  flllK/^'V^l  COSqi 
(F)  (F)  ,  (F)  (F) 

“12«/  ^IlmlZ  +  ^U^IJ  ^ImlZ 


(5.6.11) 


(5.6.12) 


Step  2:  Expression  for 

Substituting  the  above  equation  in  equation  (5.6.2)  ,  we  obtain 
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^2lt\  4-  A’^22 
t-  X32 


Here; 


V  —  T  ^  T 

At,  -  i>2ie/„aA'  f  ■^•ne//,„iA- 

Y  _  r  T 

A 12  -  ^22^IrnlX  ^12^//„aX 

X21  =  ■^2ieiCiy  +  ■i'llCnilV 
A22  =  L22t\^niy  +  ^:2^;//,Li1' 
A31  =  +  Luc[^rnlZ 

^22  =  L22e\\,\z  +  i'nCulaz  . 

Step  3:  Expression  for 

Equations  (5.4.16)  and  (5.4.18)  yield 


A’litj  +  Ai3 

A2ltl  +  A23 

[  A'3]ti  +  A'^33 


(5.6.13) 


(5.6.14) 


(5.6.15) 


Here, 


A'ls  =  Ai2  -  y,„i  sin  7 


A’23  =  A''22  +  A'„,i  sin  7  -  Zrni  cos  7  > 


(5.6.16) 


A^33  =  A'^32  +  Ym\  cos  7 
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Step  4:  Expressions  for  triple  products  in  equation  (12. 2)  for  033 


[n  =  Exit\  ^  Ei2h  +  £13  (5,6.17) 

where 

£11  =  ^rniy-^'^ll  +  ’^rnl.Y^^'21 

£12  =  ^TniyA"i2  -  nmix''^22  “  '^ni.Z^21  COS  pj  +  n„,iy  A^3i  COS  pj  (5.6.18) 

£i3  =  ~il^:n\Z^22  ~  ^mlV' -^32)  COS  pj 

[n  =  121^1  A22  (5.6.19) 

where 


V21  =  -n,„iA-A’'2i  sinpi  +  nmiy(Au  sinpi  -  A%i  cos  71)  +  n,„i^A2i  C0S7] 

> 

A22  =  -n„nA'A'23sin7i  +  rimiy(Ai3 sin 71  -  A''33Cos7i) 


(5.6.20) 


In 


—  lutl  +  ^12 


(5.6,21) 


where 


^11  —  ^miy ( Afni  sin 7i  £miCOS7i)  sin 7] 

^12  =  sinpiYn 


(5.6.22) 


Step  5;  Expression  for  the  last  term  in  equation  (1-2.2)  for  033. 

We  havi  to  differentiate  between  two  derivatives:  mj]  and  m^j.  The  first  one,  mj],  is  applied 
to  provide  a  parabolic  function  of  transmissions  errors  for  the  case  of  meshing  of  the  generated 
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pinion  and  the  gear.  Such  a  function  is  very  useful  because  it  will  allow  to  absorb  linear  functions 
of  transmission  errors  caused  by  the  gear  misalignment.  The  other  derivative,  rn'p-^,  means  that 
the  cutting  ratio  in  the  process  for  pinion  generation  is  not  constant  and  it  is  just  an  additionad 
parameter  of  machine-tool  settings. 

In  the  approach  proposed  in  this  research  project  it  is  not  required  to  have  modified  roll. 
However  the  use  of  such  parameter  in  the  more  general  case  with  5^  0  is  also  included  to  offer 
an  extra  choice.  After  some  derivations,  we  obtain 


^tT  *)  —  ■^11^1  +  -^12^1  +  ^13 


(5.6.23) 


Here 

(2C)(n,„ix A'li  -f-  -I- 

(2C)[n,„ix A'ls  -t-  n,„iyA''23  +  +  sin7i(n,„i_x'An  i-  ATrii>'A2i  +  Ami2A'’3i)]  ■ 

(2C)sin7i(nT„iAAi3  +  ^miY^zz  +  ^TnizA'33) 

(5.6.24) 

where 


2^11  = 
Zi2  = 
•^13  = 


2C  = 


{mFiV 


(5.6.25) 


Step  6:  Final  expression  for  033. 

Using  the  expressions  received  in  steps  4  and  5,  we  obtain  the  following  expression  for  033 


<133 


Zjtj  -t-  Z2t  Z3 


(5.6.26) 
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where, 


Z\  —  L\i  —  Ell  +  ^11 


Z2  =  2K|  ^  1^21  L 22  +  L11L12  ~  Ei2  ~-  ^'21  +  ^  11  "t"  Zi2 


Zz  =  LI2  +  “  ^13  “  ^^22  +  ^12  +  ^13 


Step  7:  New  representations  of  coefficients  013  and  023- 
Equations  (5.6.3),  (5.6.7)  and  (5.6.8)  yield 

ai3  =  N2iti  +  N22 

0.23  —  Niiti  +  Ni2 

Here; 


_  Kjf^Lil  +  M21 

■N12  =  «//  Lu  +  -^22 
N21  =  L21  +  Mil 

N22  =  'f  /  ^  L22  +  Mi2 

Step  8;  Derivation  of  squared  equation  for  /i 
Equations  (5.6.30)  and  (5.4.5)  yield 


ajfj  4-  02^1  +  03  —  0 


(5.6.27) 


(5.6.28) 


(5.6.29) 


(5.6.30) 


(5.6.31) 


(5.6.32) 
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where, 
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,,  ^’^21^1  +  ^23  ~  ^rnlT^Fl 

Agi  =  - 

mpi  cos  7i 


(5.6.38) 


Knowing  and  Xqx,  we  may  represent  the  position  vector  of  the  center  of  head-cutter  w’ith 
respect  to  the  cradle  center  as  follows, 

A'gi  cos  7] 

=  (5.6.39) 

Xci  sin  7i  4-  Agi 

In  practice,  the  position  of  the  center  of  the  head  cutter  is  defined  by  radial  setting  5^]  and 
cradle  angle  qx,  which  may  be  determined  by  the  following  equations, 

_  (5.6.40) 

q,=sM-~) 

Orl 

Since  the  cutter  center  Of.  must  lie  in  the  machine  plane,  the  component  must  be  zero. 

Thus,  the  sliding  base  Xb\  niay  be  determined  as, 

Xbi  =  -Agi  sin7i  (5.6.41) 
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6  Tooth  Contact  Analysis 
6.1  Introduction 

The  tooth  contact  analysis  (TCA)  is  directed  at  simulation  of  meshing  and  contact  for  misaligned 
gears  and  enables  to  determine  the  influence  of  errors  of  manufacturing  , assembly  and  shaft  deflec¬ 
tion.  The  basic  equations  for  TCA  are  as  follows; 

=  ^u\^G,<t>p.02)  (6.1.1) 

n\]\eF,4>FJi)  =  n\^,\0G.4>p-<t^2)  (6-1.2) 

Equations  (6.1.1)  and  (6.1.2)  describe  the  continuous  tangency  of  pinion  and  gear  tooth  surfaces 
El  and  Ej.  The  subscript  h  indicates  that  the  vectors  are  represented  in  fixed  coordinate  system 
Sh-  The  superscripts  1  and  2  indicate  the  pinion  tooth  surface  Ei  and  gear  tooth  surface  E2, 
respectively.  Vector  equation  (6.1.1)  describes  that  the  position  vectors  of  a  point  on  Ei  and  a 
point  on  E2  coincide  at  the  instantaneous  point  of  contact  M;  vector  equation  (6.1.2)  describes 
that  the  surface  unit  normals  coincide  at  M . 

Parameters  Op  and  <t>F  represent  the  surface  coordinates  for  Ei;  Be  and  Op  are  the  surface 
coordinates  for  E2.  Parameters  d>j  and  d»2  represent  the  angles  of  rotation  of  the  pinion  and  gear 
being  in  mesh. 

Two  vector  equations  (6.1.1)  and  (6.1.2)  are  equivalent  to  five  independent  scaler  equations  in 
six  unknowns,  which  are  represented  as 

MOf,  <t>F,  4>[,  Og,  Bp,  4>'2)  =  0  (i  =  1,2,-  -,5)  (6.1.3) 
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The  continuous  solution  of  equations  (6.1.3)  means  determination  of  five  fimctions  of  a  param¬ 


eter  chosen  as  the  input  one,  say  (p[.  Such  fimctions  are: 

^Fi4>i)^  4>f{4>\),  <Ppi4>\),  02(<z>i).  (6.1.4) 

In  accordance  with  the  theorem  of  Implicit  Function  System  Existence  [4],  solution  (6.1.4)  exists 
if  at  any  iteration  the  following  requirements  are  observed; 

(i)  There  is  a  set  of  parameters 


P{0f,  4>F,  <t>p^  <t>2) 


(6.1.5) 


that  satisfies  equations  (6.1.1)  and  (6.1.2) 

(ii)  The  Jacobian  that  is  taken  with  the  above  mentioned  set  of  parameters  emd  with  dj  as  an 
independent  variable,  differs  from  zero,  i.e. 

^  0  (6.1.6) 

{Of,  0f,  dc,  <Pp,  02) 

The  solution  of  the  system  (6.1.3)  of  nonlinear  equations  is  based  on  apphcation  of  a  subroutine, 
such  as  DNEQNF  of  the  IMSL  software  package.  The  first  guess  f  >r  the  starting  the  iteration  process 
is  based  on  the  data  that  are  provided  by  the  local  synthesis. 

The  tooth  contact  analysis  output  data,  functions  (6.1.4),  enable  to  determine  the  contact  path 
on  the  tooth  surface,  the  so  called  line  of  action,  and  the  transmission  errors. 

The  contact  path  on  pinion  tooth  surface  is  determined  in  5i  by  the  following  functions 
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4>f,  4>i)  , 


(6.1.7) 


Similarly,  the  contact  path  on  gear  tooth  surface  is  represented  by  functions 

^2(^0,  <t>p.  02)  .  0g{o\)  .  Op(o\)  (6.1.8) 

Function  C>2{0\)  relates  the  angles  of  rotation  of  the  gear  and  the  pinion  being  in  mesh.  Devi¬ 
ations  of  <t>2[o\)  from  the  theoretical  linear  fimction  represent  the  transmission  errors  (see  section 
6.4).  TCA  is  accomplished  by  the  following  procedure;  (i)  derivation  of  gear  tooth  surface,  (ii) 
derivation  of  pinion  tooth  surface,  (iii)  determination  of  transmission  errors,  and  (iv)  determination 
of  bearing  contact  as  the  set  of  instantaneous  contact  ellipses. 

6.2  Gear  Tooth  Surface 

The  gear  tooth  surface  112  and  the  surface  unit  normal  have  been  represented  in  52  by  equations 
(3.1.5)  and  (3.2.10),  where  6^  is  the  parameter  of  generating  cone  and  o,,  is  the  rotational  angle 
of  the  cradle.  Coordinate  system  S2  is  rigidly  connected  to  the  gear  To  represent  the  gear  tooth 
surface  E2  and  its  unit  normal  in  fixed  coordinate  system  5/,  we  can  use  the  following  matrix 
equations: 


<t>p-  <t>2)  -  i^//.2(<2>2)j^2(^G.<5r) 


(6.2.1) 


fl'h\^G,  <t>p,  <t>2)  -  \Lh2{<t>2)\'^2{^G-Op) 


(6.2.2) 
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6.3  Pinion  Tooth  Surface 

We  will  consider  two  cases  for  generation  of  pinion  tooth  surface;  (i)  by  a  cone,  and  (ii)  by  a  surface 
of  revolution  that  is  formed  by  rotating  curved  blades. 

Generation  by  a  Cone  Surface 

Step  1:  We  recall  that  the  generating  cone  surface  and  the  surface  unit  normal  has  been  repre¬ 
sented  in  Sp  by  equations  (5.2.1)  and  (5.2.3). 

{R^j,  Sf- sin  Of )  cos 
(/?,,,,  -  sf  sin  Of)  sin  0p 

rp  =  (6.3.1) 

-sp  cos  a f 

1 

—  cos  Of  cos9p 

iff  =  -  cos  Of  sin  ^f  (6.3.2) 

-  Sf  cos  O f 

where  5f  and  6p  are  the  surface  coordinates. 

Step  2:  During  the  process  for  generation  the  cradle  with  the  mounted  cone  surface  performs  a 
rotational  motion  about  the  Z„,i-axis  and  a  family  of  cone  surfaces  with  parameter  Op  is  generated 
in  Smi-  This  family  is  represented  in  .9„,)  by  the  matrix  equation 

^rn\{SF.6p.Op)  =  Mrn\rl(Op)  fri{sp.9p)  (6.3.3) 
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where 


=  f>  +  [  5ri  cos9i  -  Sri  sin^i  0  (6.3.4) 

The  position  vector  rvi  represents  a  point  of  the  cone  surface  in  coordinate  system  Sd  ',  Sri  and 
qi  awe  the  settings  of  the  head-cutter  center  Of  in  Smi 
Matrix  [Mmici]  is  (Fig.2.1.1) 

cos<t>p  sin(pf  0  0 

-sin<^f  cos0f  0  0 

(6.3.5) 

0  0  10 

0  0  0  1 

The  unit  normal  at  a  point  of  the  generating  surface  is  represented  in  Smi  by 

n„n{^F-<t>F)  =  [^mici(<i>f  )]nci(0F)  (6.3.6) 

where  =  np. 

We  recall  that  the  generating  cone  surface  is  a  ruled  developed  surface  and  the  surface  unit 
normal  does  not  depend  on  s/r  (Parameter  Sf  determines  the  location  of  a  point  on  the  cone 
generatrix.)  Matrix  [Lmiri.  is  the  3x3  rotational  part  of  [imirii  and  is  represented  as  follows, 

cos  sin  0F  0 

Tmiri!  =  -sin^F  cosOf  0  (6.3.7) 

0  0  1 
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Step  3:  Equation  of  meshing  of  the  head-cutter  cone  with  the  pinion  tooth  surface.  The  equation 


of  meshing  is  considered  with  vectors  that  are  represented  in  Smi-  Thus: 


AF)  41F)  _ 


(6.3.8) 


Here:  is  the  sliding  (relative)  velocity  represented  as  follows 


‘ml 


ml  ‘‘■ml  )  ^ 


-d) 

ml 


(6.3.9) 


While  deriving  equation  (6.3.9),  we  have  taken  into  account  that  vector  of  angular  velocity 
of  pinion  rotation  does  not  pass  through  the  origin  Omi  of  5mi;  -Rmi  represents  the  position  vector 
that  is  drawn  from  Omi  to  a  point  of  line  of  action  of  can  be  represented  as  (Fig. 2. 1.2): 


^1  =  [A'ci  cos  7, 


•'ml 


A 


Gi  sm7i] 


(6.3.10) 


vectors  Jd)  and  are  represented  in  S,ni  as  follows 


‘^mi  =  [cos7i  0  sin7i]^  =  Ij) 


(6.3.11) 


“’ml  - 


R 


do  0  if  ( 


U) 


[F] 


ap 


R, 


ap 


^•d) 


(6.3.12) 


Equations  from  (6.3,8)  to  (6.3.12)  yield 
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(6.3.13) 


T\[6p,4>f) 

T2{0F,<i>F) 


Here: 


Ti  =  Xnmi{-E,ni  sin^i  -  ^i(sin7i  -  mf-i))  +  YmniiXs-i  cos'll  >l2(sin7i  -  tufi)) 

+  2nml(-Eml  COS71  +  yli  COS71)  (6.3.14) 

T2  =  A''„„u(sin7i  -  mf’i)sinaF  sin(6'F  +  (pp)  -  Vnmi[(sin7i  -  mfi )  sin  Of  cos(0t’  +  pp) 

-  cosaf’Cos7i]  -  Znm\  COS71  sina/:’sin(0f  +  pp)  (6.3.15) 


where 


Ai  =  RcpS'm{6F  +  Pf)  +  Sri  sin(-gi  4-  pp) 
A2  =  Rnp  COs{6f  +  0p)  +  Sri  COs{- Qi  +  pp) 


(6.3.16) 


Step  4:  Two-parametric  representation  of  surface  of  action 

The  surface  of  action  is  the  set  of  instantaneous  lines  of  contact  between  the  generating  cone 
surface  and  the  pinion  tooth  surface  that  are  represented  in  the  fixed  coordinate  system  .  The 
surface  of  action  is  represented  by  equations  (6.3.3)  and  (6.3.13)  being  considered  simxiltaneously. 
These  equations  represent  the  surface  of  action  by  three  related  parameters.  Taking  into  account 
equation  (6.3.13)  ,  we  can  eliminate  sp  and  represent  the  surface  of  action  in  two-parametric  form 

by 
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Knl  =  rml(ffF,4>F) 


(6.3.17) 


The  common  normal  to  contacting  surfaces  has  been  already  represented  in  two-parametric 
form  by  equations  (6.3.6). 

Generation  by  a  Surface  of  Revolution 

Step  1;  The  shape  of  the  blades  is  a  circular  arc  (Fig.5.5.1)  and  such  blades  generate  a  surface 
of  revolution  by  rotation  about  the  head-cutter  axis. 

The  position- vector  of  the  center  of  the  generating  arc  is  represented  in  5mi  by  the  equation 


&i^F,4>F)  =  COS  -  Sri  Sm  <?!  0]^} 


(6.3.18) 


where, 


cos<Af  sin<^F  0  0 

-  sin  (t>F  COS  4>f  0  0 
0  0  10 

0  0  0  1 


and  has  been  expressed  by  equation  (5.5.7). 

Step  2;  We  will  need  for  further  transformations  the  following  equations 


C/mi  — 


Tnlcl\*-FI 


-  sin(^F  -i-  (f>F) 
cos{9f  +  4>f) 


0 


(6.3.19) 


(6.3.20) 
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and 


^rnl 


-  cos{6r  1-  <pp) 

-  sin{6F  +  (pp) 

0 


(6.3.21) 


/  p\ 

Here;  is  the  unit  vector  of  principal  direction  I  on  the  head-  cutter  surface  and  i?  a  unit 
vector  that  is  perpendicular  to  and  the  axis  of  the  head-cutter  (Fig.  6.3.1) 

Step  3:  To  simplify  the  equation  of  meshing  we  will  represent  it  by  the  following  equatio'. 


•tr'  =  0 


(6.3.22) 


/  1  F  ^  \ 

where  '  is  the  relative  velocity  of  the  center  of  the  circular  arc  that  generates  the  head-cutter 
surface  of  revolution.  The  proof  that  (6.3.22)  is  indeed  the  equation  of  meshing  is  based  on  the 
following  considerations; 

(i)  The  relative  velocity  for  a  point  of  the  head-cutter  surface  is  represented  b}'  equation  (6.3.9), 
given  as 


J(\F) 

ml 


(-111 


ml  / 


X  T. 


ml 


+  Rjn\ 


^ml 


(6.3.23) 


We  can  represent  posi*ion  vector  r^i  for  a  point  M  as 


Fml  =  Kn\  +  P^rn\ 


(6.3.24) 
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where  p  is  the  radius  of  the  arc  blade. 

While  deriving  equation  (6.3.24),  we  have  taken  into  account  that  a  normal  to  the  head-cutter 
surface  passes  through  the  current  arc  center  C\  the  sii'n  of  p  depends  on  how  the  surface  unit 
normal  is  directed  with  respect  to  the  surface. 

Then,  we  may  represent  the  equation  of  meshing  as  follows 


AIF)  -(F) 


'■ml 


‘ml 


=  r 


{(‘^ml  -  4,fl’)  X  [(^ml  +  pnl^i)  +  {Rrnl  X 

[(‘^mi*  -  ‘^ml’)  X  [(^"ml  +  (^ml  X 
JXF.C)  -(F) 


ml 


”ml'  =  0 


(6.3.25) 


Thus,  equation  (6.3.22)  is  proven. 

Step  4:  It  follows  from  equation  (6.3.22)  that  vector  '  belongs  to  a  plane  that  is  parallel 

to  the  tangent  plane  T  to  the  head-  cutter  surface  (Fig.6.3.2).  Tliis  means  that  if  vector  is 

JT  1 

translated  from  point  C  to  M  it  will  lie  in  plane  T.  The  u.iit  vector  ^  lies  in  plane  T  already. 
Then,  we  may  represent  the  unit  normal  n„,]  by  the  equation 


^ml  (^Fi  *t*F) 


^  fmi 
l^/mi 


-(IF.C) 


X 


(6.3.26) 


where 


is  represented  as  follows, 


‘ml 


Pml 


(6.3.27) 
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X  4- [-Xgi  COS71  -XGisin7i]^j  (6.3.28) 


^IKC)  ^1)  -4F) 

^ml  Tfil 


(6.3.29) 


The  advantage  of  vector  equation  (6.3.26)  is  that  the  surface  unit  normal  at  the  point  of  contact 
is  represented  by  a  vector  function  of  two  parameters  only,  Of  and  (j>f  \  this  vector  function  does 
not  contain  the  surface  parameter  A. 

The  order  of  co-factors  in  vector  equation  must  provide  that  the  direction  of  nmi  is  toward  the 
axis  of  the  head-cutter.  The  direction  of  Umi  can  be  checked  with  the  dot  product 

A  =  n^i  •  fm\  (6.3.30) 

The  surface  unit  normal  has  the  desired  direction  if  A  >  0.  In  the  case  when  A  <  0,  the  desired 
direction  of  fim\  can  be  observed  just  by  changing  the  order  of  co-factors  in  equation  (6.3.26). 

To  determine  parameter  A  for  the  current  point  of  contact  we  can  use  the  equation, 


cos  A  —  '  *^ml 


(6.3.31) 


Step  5;  Our  final  goal  is  the  determination  in  5,„i  of  a  position  vector  of  a  current  point  of 


contact  of  surfaces  E/r  and  Ej.  This  can  be  done  by  using  the  equation, 


»Vtii(^F.<^F') 


Pin\ 


p^jnl 


(6.3.32) 
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where  p  is  the  radius  of  the  circular  arc. 

Finally,  the  pinion  tooth  surface  may  be  determined  in  5i  as  the  set  of  contact  points.  Thus: 

=  W\p][Mj,mxVTtii{^F,<i>F)  (6.3.33) 

The  unit  normal  to  surface  Ei  is  determined  in  Si  with  the  equation 

ni(^F,</)F)  =  [iij.][Tpm,]nm,  (^F,  </'f)  (6.3.34) 

Here:  f,ni(^Fi </’f)  and  n,n\{6F^<t>F)  have  been  represented  by  equations  (6.3.17)  and  (6.3.6)  for 
strEiight  blade  cutter  and  by  equations  (6.3.26)  and  (6.3.32)  for  curved  blade  cutter.  Here  (Fig. 2. 1.2): 

0  sin7i  -A'Gisin7i 

1  0  Fm, 

(6.3.35) 

0  cos7i  -(Xgi  sin7i  + 

0  0  1 

1  0  0  0  ’ 

0  cos<^i  sin  (^1  0 

[AUj,]  =  (6.3.36) 

0  -  sin  <pi  cos  <i>\  0 

0  0  0  1 


cos  7i 
0 

[Mpmx  ]  = 

-  sin7i 
0 
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where  (f>i  is  the  angle  of  the  pinion  rotation  in  the  process  for  generation.  Angles  0i  and  (t>F  (the 
angle  of  rotation  of  the  cradle)  are  related  as  follows: 

(i)  in  the  case  when  the  modified  roll  is  not  used  and  Rap  is  constant,  we  have 

01  =  Rap(l>F  (6.3.37) 

(ii)  when  the  modified  roll  is  used,  <t>i  is  represented  by  the  Taylor’s  series 

01  =  fi<i>F)  =  R,ip{<pF  -  C(pF  -  D0f  -  E(t>f-  -  F 0^f)  (6.3.38) 

where  C,  D,  E  and  F  are  the  coefficients  of  Taylor’s  series  of  generation  motion  (see  Appendix 
B). 

Step  7:  The  tooth  contact  analysis,  as  it  was  mentioned  above,  is  based  on  conditions  of  tan- 
gency  of  the  pinion  and  gear  surfaces  that  are  considered  in  the  fixed  coordinate  system  Su  (see 
section  6.1).  To  represent  the  pinion  tooth  surface  and  the  surface  unit  normal  in  we  use  the 
matrix  equations 

=  \MhM^\eF.<i>F)  (6.3.39) 

4'’  =  (6.3.40) 

Here: 
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[Mm]  = 


1  0 
0  cos 


0  0 

■  sin  0 


0  sm<^j  cosi^j  0 


0  0 


0  1 


(6.3.41) 


where  (i>\  is  the  angle  of  rotation  of  the  pinion  being  in  mesh  with  the  gear. 


6.4  Determination  of  Transmission  Errors 


The  function  of  transmission  errors  is  determined  by  the  equation 


^(<^l)  -  [4>2  -  (<^2)°]  “  ~  ('^1)°] 


(6.4.1) 


Here:  (<^')°  (i  =  1,2)  is  the  initial  angle  of  gear  rotation  with  which  the  contact  of  smfaces  Ei 
and  E2  at  the  mean  co’.tact  point  is  provided.  Linear  function 


Nx 


-(</>!)”] 
A  2 


(6,4.2) 


provides  the  theoretical  angle  of  gear  rotation  for  a  gear  drive  without  miscilignments.  The 
range  of  4>'2  is  determined  as  follows 


(‘^2)°  X  -  *^2  -  (*^2)^^  Y 
A  2  ^’2 


'  ^o 


(6,4.3) 
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The  function  of  transmission  errors  is  usually  a  piecewise  periodic  ftmction  with  period  equal 
2tt 

to  4>i  =  —  {i  =  1,2)  (Fig.6.4.1).  The  purpose  of  synthesis  for  spiral  bevel  gears  is  to  provide  that 
iV  j 

the  function  of  transmission  errors  will  be  of  a  parabolic  type  and  of  a  limited  value  S'  (Fig.6.4.1). 

The  tooth  contact  analysis  enables  to  simulate  the  influence  of  errors  of  assembly  of  various 
types,  particularly,  when  the  center  of  the  bearing  contact  is  shifted  in  two  orthogonal  directions 
(see  section  7). 

6.5  Simulation  of  Contact 

Mapping  of  Contact  Path  into  a  Two-Dimensional  Space 

It  was  mentioned  above  that  the  contact  path  on  the  pinion  and  gear  tooth  surfaces  is  determined 
with  functions  (6.1.7)  and  (6.1.8),  respectively.  For  the  purpose  of  visualization  ,  the  contact  path 
on  the  gear  tooth  surface  is  mapped  onto  plane  (^Yclr)  that  is  shown  in  Fig.6.5.1.  The  Xc-axis 
is  directed  along  the  root  cone  generatrix  and  is  perpendicular  to  the  root  cone  generatrix  and 
passes  through  the  mean  contact  point  (Fig.  6.5.1). 

Consider  that  a  current  contact  point  N*  is  represented  in  52  (Fig. 6. 5. 2)  by  coordinates: 
X2{<l>'2),  RL'{<t>2)  where  is  the  angle  of  rotation  of  the  gear  and  RL'  —  \EN\  =  [Y'2  +  2^1)’ ■ 
Axis  A"'2  belongs  to  plane  (X,;,  Ic)  (Fig. 6. 5. 2).  While  mapping  the  contact  path  onto  plane  (A'cl^), 
we  will  represent  its  current  point  N'  by  N  that  can  be  determined  by  coordinates  A’2  and  RL', 
where  RL'  =  (Fig.6.5.3).  The  coordinates  of  mean  contact  point  M,  XL  and  R, 

have  been  previously  determined  by  equations  (3.2.1)  and  (3.2.2).  Drawing  of  Fig.6.5.3  yield 


0,N  ^0r02  +  02E  +  EN  (6.5.1) 


Here: 
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(6.5.2) 


0c02  =  0,K  +  KO2 


KO2  =  -RLcosi'jk  -  72)ic 


(6.5.3) 


where  ')k  is  determined  by: 


Ik 


(6.5.4) 


Equations  from  (6.5.1)  to  (6.5.8)  yield 


OcK  =  -i0/i02|sin72ic 


(6.5.5) 


O2E  =  A''2COS72tc  -  Ar2sin72jc 


(6.5.6) 


EN  =  RL'ism-j/zic  +  cos72ic) 


(6.5.7) 


OcN  =  Xcic  +  ^cjc 


(6.5.8) 


A,  =  A2(</>2)cos72  +  RL' {02)  sm'f2  -  [(X Lf  +  {RL)^]'^  005(7^  -  72) 
y;  =  A2(<^2)  **>^72  +  RL' {02)  ^2  -  An  sin  72 


(6.5.9) 


Contact  Ellipse 
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Theoretically,  the  tooth  surfaces  of  the  pinion  and  the  gear  are  in  point  contact.  However,  due 
to  the  elastic  deformation  of  tooth  surfaces  their  contact  will  be  spread  over  an  eUiptic  area.  The 
dimensions  and  orientation  of  the  instantaneous  contact  ellipse  depend  on  the  elastic  approach  ft  of 
the  surfaces  and  the  principal  curvatures  and  the  angel  formed  between  principal  directions 

and  of  the  surfaces.  The  elastic  approach  depends  on  the  magnitude  of  the  applied  load. 
The  value  of  b  can  be  taken  from  experimental  results  and  this  will  enable  us  to  consider  the 
determination  of  the  instantaneous  contact  ellipse  as  a  geometric  problem.  Usually,  the  magnitude 
h  is  taken  as  6  —  0.00025  inch. 

In  our  approach  the  curvatures  and  principal  directions  of  the  pinion  and  the  gear  are  determined 
with  the  principal  curvatures  and  directions  of  the  generating  tools  and  parameters  of  relative 
motion  in  the  process  for  generation. 

Gear  Tooth  Principal  Curvatures  and  Directions 

The  procedure  for  determination  of  gear  tooth  principal  curvatures  and  directions  was  de¬ 
scribed  in  section  1.2.  Knowing  functions  0p(e!>2), £>,,(<>2)  from  the  TCA  procedure  of  computa¬ 
tion,  we  are  able  to  determine  the  position  vector  rm2(^p(d>2)- <^^(^2))  smface  unit  normal 

fiTn2{^pi<t>2)-<t>j,{02))  for  an  instantaneous  point  of  contact.  The  principal  directions  and  curvatures 
for  the  generating  surface  can  be  determined  from  equations  (5.2.4),  (5.2.5)  and  (5.2.6).  The  pa¬ 
rameters  of  relative  motions  in  the  process  for  generation  can  be  determined  with  equations  (3.1.12) 
and  (3.1.13). 

Pinion  Tooth  Principad  Curvatures  and  Directions 

As  it  was  mentioned  above,  the  pinion  tooth  surface  can  be  generated  by  a  cone  or  by  a  surface 
of  revolution.  The  derivation  of  principcd  curvatures  and  directions  on  the  pinion  tooth  surface 
is  based  on  relations  between  principal  curvatures  and  directions  between  mutually  enveloping 
surfaces  Ef  of  the  head-  cutter  and  Tj  of  the  pinion.  The  procedure  of  derivation  is  as  follows: 

Step  1:  We  represent  in  Smi  the  principal  directions  on  the  head-  cutter  surface  Ef  using  the 
following  equations 
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O' =  /,//) 


(6.5.10) 


Step  2:  Parameters  of  relative  motion  in  the  process  for  pinion  generation  have  been  represented 
by  equations  (5.4.14)  to  (5.4.19).  The  derivative  of  cutting  ratio,  is  equal  to  zero  for  the  case 
when  the  modified  roll  is  not  used,  and  can  be  determined  when  the  modified  roll  is  applied  as 
follows  (see  the  Appendix) 


m 


Fl 


d^0F 

d0\ 


f'M 

[fiorW 


(6.5.11) 


where, 


/  id>F)  =  diap(l  —  2C0F  —  -  4E0^  —  5F0f-) 


f"(0!r)  =  -Eari^C  +  6D0F  ~  12E0}r  +  20Fc>f ) 


(6.5.12) 


Step  3:  Now,  since  the  principal  curvatures  and  directions  on  E/r  are  known  and  the  relative 
motion  is  also  known,  we  can  determine  for  each  point  of  conta».t  path  the  principal  curvatures  k/ 
and  K//  of  the  pinion  tooth  surface  Tj,  the  angle  ^nd  the  principal  directions  on 

Ej.  We  use  for  this  purpose  equations  (1.2.6)  to  (1.2.10).  The  principal  directions  on  Ei  can  be 
represented  in  coordinate  system  5/,  by  the  matrix  equation  (Fig. 5. 2. 2), 
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(i  =  /.//) 


(6.5.13) 


Orientation  and  Dimensions  of  the  Instantaneous  Contact  Ellipse 

Knowing  the  principal  directions  and  principal  curvatures  for  the  contacting  surfaces  at  each 
point  of  contact  path,  we  can  determine  the  half-axes  a  and  b  of  the  contact  ellipse  and  angle  q*^  ' 
of  the  ellipse  orientation  (Fig. 6. 5. 4).  The  procedure  of  computation  is  as  follows  [4]: 

Step  1;  Determination  of  a  and  b 


1  _  r.'(2)  _  .rj 


A=  -  I\\  -  h\  -  ygl  -  2gig2  cos2«t  +  g] 


^  -  \  +  \/9\  -  2gig2Cos2a  -f  gl 


(6.5.14) 


(6.5.15) 


(6.5.16) 


(6.5.17) 


where. 


4"  =  A-;"  +  A-|"'  j,=  Ai''-Aly'  (.  =  1,2) 


(6.5-18) 


Step  2:  Determination  of  (Fig. 6. 5. 4) 


66 


sincr'^^*  =  [  n,,  4*  ^n! 


cos  a 


(12)  _  ^1)  .  ;:^2) 


-Ih 


-Ih 


(6.5.19) 


Step  3:  Determination  of 

Angle  determines  the  orientation  of  the  long  axis  of  the  contact  ellipse  with  respect  to  4jV 
(Fig. 6. 5. 4)  and  is  one  of  the  angles  determined  by  the  following  equations, 


tan2o*’  *  = 


g2  sin 

9i  -  92  cos  2<r(^^l 


(6.5.20) 


Step  4:  The  orientation  of  unit  vectors  ij  and  (J  of  long  and  short  axes  of  the  contact  ellipse 


(Fig, 6. 5. 4)  with  respect  to  the  pinion  principal  directions  is  determined  with  the  equations 


-ti) 

rih  =  C/,/  cos  o 


(1)  _  rll) 


sin  Q 


(1) 


(6.5.21) 


-<i) 
*  ui 


sin  a 


1 ) 


,41) 


v(') 


(6.5.22) 


Step  5:  In  order  to  visualize  the  contact  ellipse  we  represent  its  axes  of  contact  ellipse  in  plane 


(AcTr)  (Fig. 6. 5.1),  using  the  following  equations 


%  =  [^2hiWi  Cz  =  li-2lil(,h 


(6.5.23) 


where 
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The  unit  vectors  of  axes  of  contact  ellipse  form  in  plane  (Xc,Tc)  the  following  angles  with  the 
X^-axis  (the  generatrix  of  the  root  cone)  : 

Ti  =  A„  -  72  r2  =  A|^  -  72  (6.5.26) 
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7  V  and  H  check 


The  purpose  of  the  so  called  V  and  H  check  is  the  computer  aided  simulation  of  the  shift  of 
the  bearing  contact  to  the  toe  and  to  the  hill  of  the  gear.  The  gear  quahty  is  judged  with  the 
sensitivity  of  the  shape  of  the  contact  pattern  and  the  change  in  the  level  of  transmission  errors  to 
the  above-mentioned  shift  of  contact. 

7.1  Determination  of  V  and  H  values 

Fig. 7. 1.1  shows  the  initial  position  M  of  contact  point  (it  is  the  mean  contact  point),  and  the  new 
position  M*  of  the  contact  point).  The  shift  of  the  contact  pattern  was  caused  by  the  deformation 
under  the  load.  Coordinates  X  L  and  RL  determines  the  location  of  M.  For  the  following  derivations 
we  will  use  the  following  notations. 

(i)  PF  =  A  -  A*  is  the  shift  of  the  center  of  bearing  contact,  where  F  is  the  tooth  length 
measuring  along  the  pitch  line;  p  is  an  algebraic  value,  that  is  positive  when  .4’  <  .A  and  the  shift 
is  performed  to  the  toe  as  shown  in  Fig. 7. 1.1.  Usually,  p  is  equal  to  0.25. 

(ii)  Sg  and  qq  are  the  gear  dedendum  and  addendum  angle. 

(iii)  PD  =  bo  and  P* D*  =  are  the  gear  dedendums  that  are  measured  in  sections  I  and  I*. 

(iv)  hm  =  BD  and  h*  =  B* are  the  gear  tooth  heights. 

(v)  r2  is  the  pitch  cone  angle 

The  determination  of  V  and  H  for  point  contact  M*  is  based  on  the  following  procedure. 

Step  1;  Determination  of  A"Z*  and  RL* . 

Fig. 7. 1.1  results  in  . 


h*  -  hjn  ~  pF(tan(5G  +  tanoc) 


(7.1.1) 
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-  pF  tan  6g 


(7.1.2) 


where  =  P*D*  and  be  =  PD 

h*  c  h  c 

We  assume  that  M'D*  =  -  and  MD  = - ,  where  c  is  the  clearance. 

2  2 

Taking  into  account  that 


O2A/’  =  O2P’  +  P'M' 


(7.1.3) 


we  obtain 


L »  I  ^ 

X L*  =  A*  cos  r2  4-  P” M*  sin  r2  =  >1*  cos  r2  -t-  {b^ - - — )  sin  r2 


(7.1.4) 


RL*  =  A*  sinr2  +  P' M'  cos  r2  =  A*  sin  r2  -  [be  -  - — )  cos  r2 


(7.1.5) 


The  surface  coordinates  (^e,<pp)  can  be  determined  by  solving  the  following  two  equations, 


A^2(^g,C)  = 


(7.1.6) 


iy2(&b, </>;)!'  +  =  (Frf 


(7.1.7) 
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Step  2:  Determination  of  V  iid  H 


We  introduce  the  shift  of  the  bearing  contact  in  coordinate  system  Sh  by  V  and  H  that  are 
directed  along  the  shortest  distance  between  the  pinion  and  gear  axes,  and  the  pinion  axis,  respec¬ 
tively  (Fig. 7. 1.2).  V  is  positive  when  the  gear  is  shifted  apart  from  the  pinion  in  1),  direction,  H  is 
positive  when  the  pinion  is  withdrawn.  It  is  obvious  that 


+  Vji. 


(7.1.8) 


Hn 


(7.1.9) 


Here:  (i  -  1,2)  is  the  position  vector  for  the  initial  point  of  contact,  (i  =  1.2)  is  the 

position  vector  for  the  shifted  contact  point;  and  Ivi  are  the  unit  vectors  of  coordinate  axes 

5,. 

Equations  of  tangency  at  the  new  contact  point  provide 


f^2) 


(7.1.10) 


(7.1,11) 


Gear  surface  ~>ordinates  0^  and  C>*  can  be  determined  from  equations  (7.1.6)  and  (7.1.7) 
Equations  (7.1.10)  and  (7,1.11)  yield 


(7,1.12) 
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(7.1.13) 


(1)^ 


(7.1.14) 


sin  (^)■^  = 


(2)  (1) 
+  ^hZ^iY 


(AyV  +  (n'/z) 


(1)\2 


COS  <t)i  = 


'  _  ^hY^lY  +  ^ 


(2),(1) 


("'.V)’  +  ("Vz)' 


(1)n 


(7.1.15) 


4x(^g^<^;><^2)]*  -  =  0 


(1), 


(7.1.16) 


Equations  from  (7.1.12)  to  (7.1.16)  represent  a  system  of  five  independent  equations  in  six 
unknowns;  V,  H,  <p2,<f)\,0f  and  <pf.  The  sixth  independent  equation,  that  is  required  for  the 
solution  of  unknowns,  can  be  derived  based  on  the  condition  that  the  equation  of  meshing  must  be 
satisfied  with  the  designed  gear  ratio,  i.e., 

=  /(^G-<^*G,^2^<i>'l-^F.<i>p.r,i/)  =  0  (7.1.17) 


In  solving  the  above  system,  we  first  solve  a  sub-system  composed  of  equations  (7.1.14),  (7  1.16) 
and  (7.1.17)  for  (1>2,<^f  and  Of,  and  then  calculate  the  values  of  V  and  H  directly,  by  equations 
(7.1.12),  (7.1.13)  and  (7.1.15). 
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7.2  Tooth  Contact  Analysis  for  Gears  with  Shifted  Center  of  Bearing  Contact 

After  the  determination  of  pareuneters  V  and  H,  the  tooth  contact  analysis  for  gears  with  shifted 
center  of  bearing  contact  can  be  performed  similarly  to  the  analysis  described  in  sections  6.4.  and 
6.5.  The  initial  guess  for  the  first  iteration  in  the  procedure  of  computations  is  provided  by  the  set 
of  six  unknowns  obtained  in  section  7.1. 
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Appendix  A 

Generation  witli  Modified  Roll 

1  Introduction 

Modification  of  roll  or  sometimes  called  modified  roll  means  that  the  cutting  ratio  is  not  constant 
but  varied  in  the  process  for  generation.  The  variable  cutting  ratio-the  variable  ratio  of  roll-  can 
be  provided  by  a  cam  mechanism  of  the  transmission  of  the  cutting  machine  or  by  the  servo-motors 
of  a  computer  controlled  cutting  machine.  According  to  the  developments  of  Gleason,  the  TCA 
program  can  analyze  the  process  for  generation  up  to  members  of  the  fifth  order.  However,  due  to 
the  limitations  caused  by  application  of  cam  mechanisms  only  the  pareuneters  up  to  the  third  order 
are  controllable  in  the  process  for  generation. 

The  modified  roll  is  an  additional  parameter  for  the  synthesis  of  spiral  bevel  gears.  In  our 
approach  the  synthesis  of  spiral  bevel  gears  can  be  performed,  as  it  was  mentioned  above,  with  a 
constant  cutting  ratio.  However,  we  consider  in  this  section  the  apphcation  of  modified  roll  as  weU 
to  provide  a  broader  point  of  view  on  synthesis  of  spiral  bevel  gears. 

2  Taylor  Series  for  the  Function  of  Generation  Motion 

According  to  the  practice  of  Gleason,  the  kinematic  relation  between  the  angles  of  rotation  of  the 
workpiece  and  the  cradle  is  represented  by  a  Taylor’s  series  up  to  fifth  order.  To  the  knowledge  of 
the  authors,  Gleason  has  never  published  any  materials  related  to  the  kinematics  of  the  modified 


roll.  However,  Professor  Zheng  had  done  a  good  job  in  deciphering  Gleason’s  mechanisms  for 
modified  roll  and  represented  the  kinematic  relations  in  his  valuable  book  [5]. 


Consider  that  the  angles  of  rotation  of  the  pinion  and  the  cradle  are  related  by 
function 


<J>1  =  f{<i>F)  fec^  {K  >  3) 

We  assume  that  (f>i  =  0  at  (pf  —  0  and  represent  f{<pF)  in  the  neighborhood  of  (pp 
Taylor  series  as  follows, 


<t>\  —  f  {0)<Pf  +  •  •• 


Taking  into  account  that 


d<Pi 

d<pF 


f 


We  obtciin 


/'(O) 


u> 


(1) 


(F) 


u' 


|<^p=:0 


where  R„p  is  the  ratio  of  roll. 


Without  loosing  generality  of  the  solution,  we  can  take  =  1  and  then  obtain 


nonhnear 


(A.l) 


=  0  bv  the 


(A.2) 


(A.3) 


(A.4) 
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/Vf) 


d<t>F 

dt 


^  1 


(A.5) 


Differentiation  of  equation  (A.5)  yields 


d<i>F 

dt  ’ 


=  /'(0f) 


df^<i>F 

dt^ 


(A.6) 


Equation  (A.6)  yields 


^  _  /"(d>F) 
4  ~  /'(<^f) 


(A.7' 


where  oj  - 


(i^<i>F 


is  the 


angular  acceleration  of  the  cradle. 


Equation  (A.7)  with  new  designations  can  be  represented  as  follows 


R 


7  (0) 


ap 


(A.8 


Similar  differentiation  of  higher  order  of  equation  (A. 3)  yields: 


<**i  —  Rap{*PF  ~  ^  “■  Dd’f  ~  E<P^f-  — 


2C  =  -^/"(0) 

^ap 


QD^-~r{Q) 

/lap 


24£=  -^^(0) 

/lap 

120F=  -^/‘'(O) 

/lap 


Unfortunately,  function  f{<t>F)  caimot  Le  represented  in  explicit  form  for  certain  cutting  ma¬ 
chines,  for  instance,  for  the  Gleason  spiral  bevel  grinder.  For  such  a  case  we  will  consider  the 
following  cuxiliary  expressions 


d^<t>F 


QCX  = 

^F 


(A.IO) 


04  = 


d^0F 

IF 


2iDX  = 


04 
w  p 


(A.ll) 


d^<pF 

dt^ 


\2^)EX 


(A. 12) 


Then,  differentiating  equation  (A. 6)  and  taking  Of  =  0,  we  may  obtain  the  following  equations 


6D  =  6CX  -  3(2Cf 


(A. 13) 
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2iE  =  24DX  +  (2C)[15(2C)^  -  10(6CA'')] 


(A. 14) 


120F  =  120EX  -  l5(2C)i24DX)  +  105(2r)2[6rA  -  (2Cf\  -  10(6rA)^  (A. 15) 


The  procedure  for  determination  of  coefficients  C,  D,  E  and  F  for  the  Taylor’s  series  (A. 9)  when 
function  /i (</>/’)  cannot  be  represented  in  explicit  form  is  as  follows; 

Step  1:  Differentiate  the  implicit  equation  that  relates  4>p  and  4>i  up  to  five  times  and  then  find 

02, 03,  04,  05  in  terms  of  <pi  and  </>2  at  </>!  =  4>f  =  0. 

Step  2:  Considering  ((>1  =  =  0,  find  6CX,24DX,120EX  by  equations  (A. 10)  -  (A. 12). 

Step  3:  Find  2^,  6D,24E,  120F  by  equation  (A. 8), (A. 13)-  (A. 15). 

3  Synthesis  of  Gleason’s  Cam 
Introduction 

Glecison’s  cam  mechanism, as  shown  schematically  in  Fig.  A. 3.1,  is  an  ingenious  invention  that 
has  been  proposed  and  developed  by  the  engineers  of  the  Gleason  Works.  The  mechanism  trans¬ 
forms  rotation  of  the  cam  about  O,  into  rotation  of  the  cradle  about  Or  The  rotation  of  the  cam 
about  0,j  is  related  with  the  rotation  of  the  pinion  being  generated,  but  the  angles  of  cam  rotation 
and  pinion  rotation,  <pq  and  <f>i,  are  related  by  a  linear  function  when  there  is  no  cam  settings. 

To  authors’  knowledge,  the  engineers  of  the  Gleason  Works  have  not  published  the  principles 
of  synthesis  and  analysis  of  this  mechanism.  However,  H. Cheng  [6],  Zheng  5  have  made  good 
contributions  to  the  deciphering  of  this  mechanism.  The  following  is  a  systematic  representation 
of  synthesis  and  analysis  of  Gleason’s  mechanism 

The  purpose  of  cam  synthesis  is  to  obtain  the  shape  of  the  cam,  considering  that  the  angles  of 
rotation  of  the  cam  and  the  cradle  are  related  by  a  linear  function,  However,  this  function 
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can  be  modified  into  a  nonlinear  function  by  changing  the  location  of  the  designed  cam  with  respect 
to  Og  and  the  orientation  of  the  cam  guides  that  are  installed  on  the  cradle.  Fig.  A. 4.1  shows  the 
settings  of  the  cam  mechanism  with  the  designed  shape  :  (i)  the  cam  is  translated  along  the  line 
OcO(j  an  amount  AT  ;  (ii)  and  then,  the  cam  guides  are  rotated  about  the  cam  rotation  center 
and  formed  angle  q  with  Or.Oq-  It  is  obvious  that  the  cam  mechanism  with  the  settings  AF  and 
a  will  transform  rotation  about  0^  to  Or.  with  a  nonlinear  function  between  the  angles  of  rotation 
of  the  cam  and  the  cradle.  The  deviation  of  this  function  from  a  linear  one  depends  on  settings  of 
the  cam  mechanism  and  will  be  discussed  in  section  A. 4. 

Coordinate  Systems 

While  considering  the  synthesis  of  the  cam  mechanism,  we  will  use  three  coordinate  systems: 
the  movable  coordinate  systems  Sr  and  S,  that  are  rigidly  connected  to  the  cradle  and  the  cam, 
and  Sf  that  is  the  fixed  coordinate  system  (Fig.  A. 3.2). 

Equation  of  Meshing,  Contact  Point  in  Sr 

Assuming  that  the  transformation  of  motion  is  performed  with  constant  ratio  of  angular  veloc¬ 
ities  and  in  the  same  direction,  we  can  determine  the  location  of  instantaneous  center  of  rotation, 
7,  in  coordinate  system  Sj  by  using  the  equation  (Fig.  A. 3.2) 


(A. 16) 


Where,  E  in  the  distance  between  the  cradle  center  O,  and  the  cam  rotation  center  0,^.  is  the 
so-called  pitch  radius  of  the  cam. 

The  location  of  instantaneous  point  of  contact  M  on  the  guides  can  be  deterrruned  by  using 
the  theorem  of  planar  gearing  4  .  According  to  this  theorem  the  common  normal  to  the  guides 
and  the  cam  at  the  point  of  their  contact  must  pass  through  the  instantaneous  center  of  ri  tation 
7.  Thus,  contact  point  A7  and  the  unit  normal  at  M  are  represented  in  >V  as  fo'lows 
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r,  =  [  -6  -  u  0  1  f 


(A.17) 


=  [  1  0  0  f 


(A.18) 

(A.19) 


Here;  6  is  an  algebraic  value  (6  is  positive  if  the  left  side  of  guides  is  considered  and  b  is  negative 


if  the  right  side  of  guides  is  considered);  u  is  a  variable  parameter  that  is  determined  with  the 
equation 


u  =  (£  +  r„)  cos  Or,  -  E 


(A.20) 


Equations  (A.17)  and  (A.20)  yield 


fr[6r)  =  [-h  f{9r)  0  1 


(A.21) 


where 


f{0r)  =  E  -  {E  -r  r^)  cos  Or 


(A.22) 


Shape  of  the  Cam 

The  shape  of  the  cam  is  a  planar  curve  that  is  represented  in  5,  the  matrix  eqr.ation 
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(A.23) 


Here:  coordinate  system  5p  is  an  anxiliarj’  fixed  coordinate  system  (Fig.  A. 3. 2).  Matrices  in 
equation  (A.23)  are  represented  as  follows 


= 


cos  6q 

-  sin  Oq 

0 

0 

sin  9q 

cos  6q 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

[Mp/l  = 


10  0  0 

0  1  0  i: 

0  0  10 
0  0  0  1 


W]r\  - 


cos  Or  sin  Or  0  0' 

-  sin  Or  cos  ^^0  0 

0  0  10 

0  0  0  1. 


The  normal  to  the  cam  shape  is  represented  by  the  matrix  equation 


^q{0r')  —  Lqjr  LpJ  Ljrj'^r{0r'l 


(A.24] 


(A.25;: 


(A. 26 


Here:  [ipy]  is  the  identity  matrix  and  is  the  (3  x  3)  submatrix  of  the  respective  matrix  [M],  We 
consider  that  the  shape  of  the  cam  and  its  normal  depend  on  the  generaUzed  parameter  6^  only 
since 


Br  =  ( 


(A.28) 


The  final  equations  of  the  cam  and  its  normal  are  represented  as  follows 


+  (r„  +  f;)cos(^^,,)sin(^— -  £  sin 

-  i^u  +  E)  cosine,,)  cosi^J,,)  -  EcosB, 

0 

1 


(A.29) 


cos(^^,) 

0 


(A.30) 


4  Cam  Analysis 

The  cam  analysis  is  directed  at  the  determination  of  function  0'{9'^)  for  a  cam  and  guides  with 
modified  settings.  The  analysis  is  based  on  simulation  of  tangency  of  the  designed  cam  with  the 
cradle  guides  taking  into  account  the  settings  of  the  cam  and  the  guides. 

Cofirdinate  Systems  and  Coordinate  Transformation 

(^)ordinate  systems  S,i,  S,  and  S,-.  are  rigidly  connected  to  the  guides  and  the  cradle  (Fig. 
A. 4. 1(a)).  The  guides  after  rotation  about  O,:  form  angle  o  with  the  i/,.-axis 
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Coordinate  systems  5q,  Sn  and  Sm  are  rigidly  connected  to  the  cam.  The  settings  of  the  cam 
with  respect  to  Sm  are  determined  by  AT  and  angle  a. 

The  cam  and  the  cradle  perform  rotations  about  Om  and  Of,  respectively  (Fig.  A. 4. 2).  The 
conditions  of  continuous  tangency  mean  that  the  designed  cam  and  the  guides  have  a  common 
normal  and  a  common  position  vector  at  every  instcint  in  5y. 

A  current  point  N  of  the  guide  is  determined  in  5y  with  the  equation  (Fig.  A. 4.1  and  Fig. 
A.4.4): 


(A.31) 


where 


=  [  -6  -  A  0  1  I’’  (A.32) 


The  unit  normal  is  determined  in  5/  as  follows 


'  =  [L*fc][^cr.][L,,d]nd 


(A.33) 


where 


n,i  =  [  1  0  0  ]’' 


(A.34) 


Here: 
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m.]  = 


cos  -  sin  6*  0  0 

sin  O'  cos  O'  [  0 

0  0  10 


0 


0  0  1 


(A.35) 


^MrJ 


[M.d\  = 


10  0  0 
0  1  0  -£ 

0  0  1  0 

0  0  0  1 


cos  Q  sin  Q  0  0  1 


sin  a  cos  q  0  0 


{A.36) 


0  1  0 


0  0  1 


(A.37) 


A  current  point  of  the  cam  and  the  unit  normal  at  this  point  are  represented  in  Sj  by  the 
equations  (Fig.  A. 4. 1(b),  Fig.  A. 4. 2). 


(A.38) 


(A.39) 
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Equations  (A. 38)  and  (A. 23)  yield 


(A.40) 

where 

=  [  -i>  -  u  0  1  (A. 41) 

Matrices  [Mqp],[Mpf]  and  [M/c]  have  been  represented  by  equations  (A. 24),  (A. 25)  and  (A. 26), 
respectively.  Matrices  and  [M*,]  are  repr<*sented  as  follows  (^ig.  a. 4, 1(b), 

Fig.  A. 4.2): 

■  1  0  0  o' 

0  1  0  -E 

[Mfp]  =  (A.42) 

0  0  1  0 

0  0  0  1 

cos  Og  -  sin  0  0 

sin  9*  cos  ^*0  0 

[m;„]  =  (A.43) 

0  0  10 

0  0  0  1 
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cos  Q 

sin  a 

0 

0 

—  sin  a 

cos  a 

0 

0 

[M,.d]  = 


0  0  10 
0  0  0  1 


(A.44) 


[KJ  = 


10  0  0 
0  1  0  -AT 
0  0  1  0 
0  0  0  1 


(A.45) 


Equations  (A. 39)  and  (A. 27)  yield 


(A.46) 


Here; 


n.  ]  =  [  1  0  0  ] 


(A.47) 


Matrices  [L*]  and  [L]  are  3x3  submatrices  of  matrices  [A/*]  and  [A/]. 

Equations  of  Tangency 

The  tangency  of  cam  and  guides  with  modified  settings  is  represented  by  equations 


(A.48) 
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(A.49) 


n^Pi6:,e,,a)  =  n^p{e;,0„,a) 

We  recall  that  vector  equations  (A.48)  and  (A.49)  yield  a  system  of  only  three  independent 
equations  in  four  unknowns:  and  A;  setting  parameters  AT  and  q  are  considered  as  given; 

and  6c  are  related  with  equation  (A. 28)  and  6g  is  considered  as  a  generalized  parameter.  Our 
goal  is  to  determine  the  function  that  relates  angles  of  rotation  of  the  cam  and  the  cradle,  (p^  and 
</>*,  and  the  parameters  of  settings  a  zind  AT,  i.e.  the  function 

F{e'^,6;,AT,a)  =  0  (A.50) 

Equality  of  Contact  Normal-Satisfaction  of  Equation  (A.49) 

Equations  (A.49),  (A. 33)  and  (A.46)  yield 

[L%][A]  =  [r^][A][L,,][Lj^i  (A.51) 

Then  we  obtain 

\B]\A]  =  [A]}c]  (A.52) 

Here: 
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■  cos{e;~ei)  sm{e;-e:)  o' 

[B] ^[L%][LUn]-  cos{e;-e;)  o  (A.53) 

0  01 

cos(0q  -  0c)  -  sin(0q  -  0c.)  0 

[C]  =  \L^j]\L]c]  -  sin(0q  -  6^)  cos(0g  -  0c)  0  (A. 54) 

0  0  1 

COSO  sin  Q  0 

[A]  =  [Led]  =  [Lng]  =  -sino  cosq  0  (A. 55) 

0  0  1. 

Matrices  (A. 52)  are  rotational  matrices  that  describe  rotation  about  axes  of  the  same  orienta¬ 
tion.  This  means  that  we  can  change  the  order  of  co-factor  matrices  and 

{B]{A]  =  fr][A]  (A.56) 

This  yields  that 

[r]-^:5][A]  =  [A]  (A,57) 

=  (A. 58) 

where  [7]  is  a  unitary  matrix  and 
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[B]  =  [C] 


(A.59) 


Equation  (A, 59)  yields 


e:  =  -{e,  -6,)  =  - 


since 


Or. 


E 


6q 


Equality  of  Position  Vectors-Satisfaction  of  Equation  (A. 48) 


Equations  (A. 48),  (A. 31)  and  (A. 38)  yield 


After  transformations  we  obtain 


[Q]fd  =  [5]fc 


Here: 


(A.60) 


(A.61) 


(A.62) 


(A.63) 
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[Q]  = 


Oil 

012 

0 

014 

021 

022 

0 

024 

0 

0 

1 

0 

0 

0 

0 

1 

On 

O12 

0 

aU 

021 

022 

0 

“24 

0 

0 

1 

0 

0 

0 

0 

1 

(A.64) 


(A.65) 


The  rotational  3x3  submatrices  of  [(?]  and  [5]  are  equal  due  to  the  equality  of  contact  normals 
(see  equation  (A. 52).  The  elements  of  [(?]  and  [S]  are  represented  by 


Oil  =  COST)  ai2  =  sin  t;  021  =  — 012  022  = 


(A,66) 


where  t)  =  6*  -  6*  +  a  ,  and 

au  =  E[-sm{0*  -  6')  +  s'mOg]  ,  024  =  f’i- cos(e*  -  P*)  ^  cos 0,:  (A. 67) 


=  -£'[sin(^q  -  q)  -  AT  sin  a]  ,  =  £';cos(0q  -  a)  -  AT  cos  a  (A. 68) 

Matrix  equation  (A. 62)  yields  the  following  system  of  two  linear  equations 


ai2(u  -  A)  +  ai4  -  aj4  =  0  ,  022(11  -  A)  +  024  -  024  =  0  (A. 69) 
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Eliminating  (u  -  A),  we  obtain 


(oi4  -  014)022  -  (024  -  054)011  -  0  (A. 70) 

Equations  (A. 70),  (A. 66),  (A.67)  and  (A. 68)  results  in 

r  AT 

FiO;j;,AT,a)  =  sin(^(^,*  -  ^*)]  +  —  sin(e;  -  91)  -  Bm{9;  -  a)  -  sina  =  0  (A.71) 

Equation  (A.71)  represents  in  implicit  form  the  displacement  function  for  the  cam  mechanism 
with  settings  a  and  AT.  It  is  easy  to  be  verified,  that  equation  (A.71)  with  AT  =  0,q  =  0 
represents  the  linear  function  , 

(A.72) 

For  Gleason’s  grinder,  E  is  equal  to  15  inch.  According  to  Gleason’s  practice,  the  sense  of 
rotation  of  the  cradle  and  the  cam  is  opposite  to  the  assumption  in  the  derivation  in  this  report. 
Without  loss  of  generality,  by  substituting  9*  =  -9*  and  8*  =  -9*  with  E  =  15  in  equation  (A.71) 
we  obtain  the  final  expression  of  the  relation  between  9*  and  9*  as  follows, 

AT  r 

sm{9:  +  a)  -  sina  +  —  sin(^*  -  9;)  +  sin  ^(^,*  -  9^)  =  0  (A.73) 
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5  Determination  of  Coefficients  of  the  Taylor’s  Series 


The  determination  of  the  coefficients  of  the  Taylor’s  Series  for  generation  motion  with  modified 

roll  is  a  lengthy  process.  Gleason  provides  its  customers  with  computer  program  which  can  select 

the  cams  with  settings  and  analyze  the  effects  of  the  modified  roll.  However  Gleason’s  program 

is  a  black  box  with  no  explanation  for  the  determination  of  the  coefficients  of  the  Taylor’s  Series. 

Valuable  contribution  to  the  understanding  of  Gleason’s  program  has  been  made  by  C.Q.  Zheng 

[5].  For  reader’s  convenience,  a  series  of  derivations  are  represented  in  this  section,  which  coincide 

with  the  equations  in  [5]  except  some  printing  errors. 

In  the  process  of  generation,  the  cam  rotates  at  a  constant  angular  velocity  .  Without  loss  of 

generalitv,  we  assume  that  the  cam  rotates  with  unitarv  velocitv,  i.e.  =  1.  Using  the  procedure 

at 

discussed  in  section  A. 2,  we  differentiate  equation  (A. 73)  five  time  as  follows, 

T  T 

uj;  cos(6>*  +  q)  +  (u;;  -  1)[-^  cos(6';  "  ^  cos  sin(6»;  -  ^i*)]  =  0  (A.74) 

02  cos(0*  +  q)  +  sin(0’  +  a) 

T  T  /\'T* 

T*  7*  /\'T^ 

-K*  -  sin  -  e;)]  +  —  sin(e;  -  0;)]  =  0  (A.75) 

03  cos(0*  +  a)  -  3o2(u.’*)^  sin(0*  +  q)  -  (a;*)  cos(0*  +  a) 

7*  V  /\T* 

+a3[^  cos  -  e;)  +  —  cos(0;  -  e;)] 

V  V  /\.T 

-302K  -  ~ 

V  V  /\  T* 

-K  -  i)"i(^)^cos  f^(e:  -9;)^—  cos{e;  -  =  o  (a.76) 
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a4  cos(9*  +  a)  -  4a3u;*  sin(0’  +  q)  -  Soj  sin(0*  +  a) 
-6(u)*)^a2  cos(0*  +  a)  +  (0;*)“*  sin(0*  +  q) 

7*  7*  7^ 

-4a3K  -  -0;)+  —  cos(0:  -0*)] 

7*  7*  /\  7^ 

-3a^[(^)^  sin  ^(0,*  sin(0:  -  0;)] 

T*  r  /\T 

+  04:  cos  ^(|9;  -  6»;)  +  —  cos(6>;  -  6'^)] 

V  7*  /\  T' 

-ea^K  -  l)^[(^f  cos^(^,*  -  ^ 

V  7*  T' 

«  -  1)^[(^)^  sm  f^{9:  -0;)^^  sin(e;  -  S,*)]  =  0 


as  cos(0*  +  q)  -  [10a3(w’)^  +  -  (u;*)^j  cos(^*  +  q) 

-[504^*  +  100303  -  10a2(u»')^]  sin(0*  +  a) 

V  7*  /K'T' 

+  05(^1  cos  -  ^g)  +  -^  -  ^g)l 

7*  j* 

-504(0-^*  -  l)[(j|)^sin^(0;  -^q)+  -  ^q)] 

V  T*  /\  'T' 

-10a3(a';  -  1)^[(y|)^cos  :^{e;  ~^q)^  '^os(^r  “  ^q)] 

T  V  /\7^ 

-10a2a3[(y^)^sin  -0*)+  —  sin(0;  -  P*)l 

7»  j» 

-i5a^K  -  I)[(^f  cos^(^;  -  0*)  ^  —  cos(0;  -  e;)] 

7*  7*  7^ 

+  1002^  -  if  i(jT)^sin  ^(6»;  -  sin(6';  -  6»;)i 

7*  ?'  /^T' 

4-(^;  -  i)^[(y|f  cos  ^(0;  -0;)^—  cos{0:  -  e;)]  =  o 


At  6*  =  6*  -  0,  we  can  determine  u;*, 03, 03, 04, 05  from  above  expressions  as  follows, 


Tu  +  ‘^T 

15  cos  a  4-  r^i  +  AT 


(A.77) 


(A.78) 


(A.79) 
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02  = 


15  sin  a 


15  cos  a  +  r„  +  AT 


.^2 


K*) 


03  = 


r®  AT 

3020;*  sina  +  (u);)3cosa  +  +  -^) 


cos  a  + 


Vu  -f-  AT 
15 


04  = 


05  = 


7.3 

6a2(u>;;)^  COSO  +  [4030;’  +  3o^  -  (u)*)^]sinQ  +  6o2(u)’  -  1)^(y^  +  -^) 


cos  a  -f 


+ AT 
15 


(A.80) 


(A.81) 


(A.82) 


COS  a  + 


;rq:^{[10o3(u>*f  +  I502U)’  -  (w^)  ]cosa 


15 


+  [5a4u;*  +  100203  -  1002(0;*)^]  sin  a  +  [1003(0;*  -  1)^ 


AT . .  r3  AT. 


V  /  \  /  V 

+  l5a’K  -  1)1( -  K  -  lf( +  ^)) 


(A.83) 


Using  equation  (A. 8)  and  (A. 10)  -  (A. 12),  we  obtain 


Rac  —  —  1  + 


15  cos  a 
r„  +  AT 


2C 


Rac  -  1 


tana 


6CA  = 


24DX  = 


1  +  3(2C)tana  +  +  AT) 


1  + 


15 COSO  ^15^ 

+  AT 
15  cos  a 


(A.84) 


(A.85) 


(A. 86) 


cos  Q  + 


Tu  + 


^{6(2C)co5a  +  [4(6CX)  +  3(2C)"  -  Ijsina 


15 


AT 

+6(2C)(1  -  Rac)\:^,  +  ^)} 


(A.87) 


94 


\2QEX 


- ^^j^-^{[10(6CX)  +  15(2C)*  -  l]cosa 

cos  a  +  “  -— - 

15 

+  [5(24DX)  +  10(2C)(6CJ>ir)  -  10(2C)]  sin  a 

_3  AT 

+  [10(6CX)(1  -  Racf  +  15(2C)2(1  -  Rac]{:^  +  — ) 
r®  AT 

-{l-Ra.njl-,+^)}  (A.88) 


Knowing  Rar,^C,6CX,2ADX  and  120£'X,  we  can  determine  6D,2AE  and  120F  by  equations 
(A.13)  -  (A. 15). 

6  Selection  of  Cams  and  Cam  Settings 

In  order  to  provide  the  desired  low  transmission  errors  and  bearing  contact,  the  ratio  of  roll  Rap 
and  second  ratio  of  roll  (2c),  which  are  determined  by  the  local  synthesis,  must  be  applied  for  the 
grinder.  Due  to  the  structure  of  Gleason’s  grinder,  the  ratio  of  roll.  Rap  is  related  to  Rac  as  follows, 

Rac  =  rnirUcRap  (A. 89) 


Here,  Ran  is  the  transmission  ratio  between  the  cam  and  the  cradle,  as  determined  by  equation 
(A. 84).  nic  is  a  fixed  gear  ratio  and  is  equal  to  1  in  Gleason’s  grinder;  is  the  gear  ratio  from  the 
workpiece  to  the  cam  and  is  determined  as. 


n 

TTlj  —  '' 

ni 


(A.90) 
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here,  n  is  the  number  of  teeth  of  the  workpiece  and  is  the  index  internal,  i.e.  the  gear  tooth 
number  skipped  over  in  indexing. 


Prom  equation  (A. 80)  and  (A. 89),  we  obtain 


a 


lr((2c)fi 


)1 


(A.91) 


+  AT  - 


15  cos  Q 
Rac-l 


(A.92) 


The  cams  and  their  pitch  radii  are  tabulized.  A  cam  with  pitch  radius  closest  to  (r^  +  AT) 
calculated  by  equation  (A.92)  should  be  selected.  After  the  cam  with  pitch  radius  is  selected 
the  corresponding  setting,  AT,  can  then  be  determined  as: 


AT  = 


15  cos  a 
Rac  -  1  " 


(A.93) 


In  some  cases,  it  is  also  necessary  to  control  6CX.  In  order  to  satisfy  Rac^2C  and  the 

value  of  Tii  can  be  used  together  with  AT  and  a.  Since  must  be  an  integral  number  it  is  difficult 
to  obtain  an  accurate  solution.  But  by  careful  selection  of  cams  and  index  interval  n^,  a  practical 
engineering  solution  is  often  achievable. 

When  the  cam  and  its  s<  itings  are  selected,  it  is  then  necessary  to  determine  the  coefficients  of 
the  Taylor’s  Series  of  the  generation  motion  and  carry  out  the  TCA  to  see  how  the  higher  order 
coefficients  (i.e.,  6D,24E  and  120P’)  affect  the  transmission  errors  and  bearing  contact.  If  the  res\ilt 
of  TCA  are  satisfactory,  then  the  gears  can  be  groimd  by  the  selected  cam  and  cam  settings. 
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Appendix  B 


Description  of  Program  and  Numerical  Example 


Input  and  Output  of  Program 

The  research  project  is  complemented  by  a  computer  program,  which  can  be  used  for  the 
determination  of  machine  tool  settings  through  the  method  of  local  synthesis  and  simulate  the 
transmission  errors  and  bearing  contact  through  TCA.  The  input  data  to  the  program  include  four 
parts. 

Part  1 .  Blan't  Data 
TNI  ;  pinion  number  of  teeth 
TN2  :  gear  number  of  teeth 
C  :  shaft  offset  (  zero  for  spiral  bevel  gear  ) 

FW  :  width  of  gear 
GAMMA  :  shaft  angle 
MCD  :  mean  Cone  distance 
RGMAl  ;  pinion  root  cone  angle 
B1  :  pinion  spiral  angle 
B2  :  gear  spiral  angle 
RGMA2  :  gear  root  cone  angle 
FGMA2  :  gear  face  cone  angle 
PGMA2  :  gear  pitch  cone  angle 
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D2R  :  gear  root  cone  apex  beyond  pitch  apex 
D2F  :  gear  face  cone  apex  beyond  pitch  apex 
ADD2  :  gear  mean  addendum 
DED2  :  gear  mean  dedendum 
WD  ;  whole  depth 
CC  ;  clearance 

DEL  :  elastic  approach  (experiment  datum) 

Part  2.  Cutter  Specifications 
RU2  :  gear  nominal  cutter  radius 
PW2  :  point  width  of  gear  cutter 
ALP2  :  blade  angle  of  gear  cutter 

Part  3.  Parameters  of  Synthesis  Condition 

FI21  :  derivatj'i'o  nf  transmission  ratio,  negative  for  gear  convex  side  and  positive  for  gear  concave 
‘'de  The  ritiigo  is  -l).008  <  F721  <  0.008  . 

KD  ;  percentage  of  the  half  long  axes  of  contact  ellipse  over  face  width.  K D  =  0  15  -  0.20  . 
ETAG  :  direction  angle  of  contact  path.  For  right  hand  gear,  -80°  <  ET AG  <  0°  for  gear  convex 
side  and  -80°  <  ETAG  <  0°  for  gear  concave  side;  For  left  hand  gear,  0°  <  ETAG  <  80°  for 
gear  convex  side  and  -80°  <  ETAG  <  0  for  gear  concave  side.  When  ETAG  is  close  to  zero,  the 
contact  path  is  along  the  tooth  height,  when  the  magnitude  is  increased,  the  contact  path  will  have 
bias  in  and  reach  almost  longitudinal  direction  if  ETAG  is  close  to  90  degrees. 

GAMAl  :  pinion  machine  root  angle,  which  is  the  same  as  the  pinion  root  angle  if  no  tilt  is  used. 
RHO  :  radius  of  the  arc  blade  if  curved  blade  is  used,  which  can  be  any  values  when  curved  blade 
is  not  used. 

C2  ;  second  order  ratio  of  roll  if  modified  roll  is  used.  C2  is  zero  without  modified  roll. 
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ALPl  ;  pinion  cutter  blade  angle.  ALPl  is  positive  for  gear  convex  side  and  negative  for  gear 
concave  side.  ALPl  can  be  the  same  as  ALP2.  For  better  result,  it  is  suggested  for  pinion  concave 
side  the  magnitude  of  ALPl  is  smaller  than  ALP2  and  for  pinion  convex  side,  the  magnitude  of 
ALPl  is  larger  than  ALP2.(As  shown  in  the  example). 

TNlI  :  number  of  teeth  skipped  over  indexing.  TNII  is  only  used  in  modified  roll,  the  ratio  between 
TNII  and  TNl  must  not  be  an  integer. 

Part  2.  Control  Codes 

JCL  ;  JCL  control  V  and  H  check,  JCL  -  1  means  no  V-H  check. 

JCH  ;  For  right  hand  gear,  set  JC H  =■  1.  for  left  hand  gear  set  JCH  =  2. 

JCC  ;  For  straight  blade,  set  JC'C  =  1,  for  cur%'ed  blade  set  JCC  =  2. 

TLl,  TL2  :  Extra  points  on  contact  path,  both  should  be  less  or  equal  than  2. 

Tne  program  output  includes:  (1)  the  machine-tool  settings  for  gear  and  pinion;  (2)  the  trans¬ 
mission  error;  (3)  the  contact  path;  (4)  the  length  and  orientation  of  the  long  axes  of  the  contact 
ellipse;  and  bearing  contact  at  toe  and  heel  position. 

Numerical  Example 

The  model  used  in  this  report  is  the  spiral  bevel  drive  with  the  shaft  angle  of  90  degrees.  In 
the  numerical  example,  modified  roll  and  curved  blad“  for  generation  of  gears  were  not  used  since 
favorable  results  were  attained  without  them.  The  list  of  the  blank  data  and  machine  tool  settings 
are  tatjulized  in  the  attached  tables. 

The  TCA  results  with  \'-H  check  are  shown  through  Fig.  B.l  through  Fig.  B.6.  The  V  and  H 
values  shown  in  the  figures  are  of  “(jj,  inch.  It  is  shown  also  that  ,he  transmission  errors  are  very 
small  and  the  bearing  contact  is  stable  for  both  side  at  the  three  positions,  toe,  mean  and  heel. 
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BLANK  DATA 


PINION 

GEAR 

NUMBER  OF  TEETH; 

11 

41 

PRESSURE  ANGLE: 

20" 

SHAFT  ANGLE; 

o 

e 

MEAN  SPIRAL  ANGLE: 

35.0" 

HAND  OF  SPIRAL; 

LF 

RH 

OUTER  CONE  DISTANCE; 

90.07 

FACE  WIDTH; 

27.03 

WHOLE  DEPTH; 

8.11 

8.11 

CLEARANCE; 

0.81 

0.81 

ADDENDUM; 

5.24 

2.061 

DEDENDUM; 

2.87 

6.05 

PITCH  ANGLE; 

i5«r 

74*59’ 

ROOT  ANGLE: 

13'20’ 

70"39' 

FACE  ANGLE: 

19021 ' 

76»40’ 

GEAR  CUTTER 

SPECIFICATIONS 

BLADE  ANGLE: 

0 

0 

CM 

CUTTER  DIAMETER; 

152.40 

POINT  WIDTH; 

2.79 

100 


GEAR  MACHINE  TOOL  SETTINGS 


RADIAL  SETTING (*) :  70.43577 

CRADLE  ANGLE (q) :  62.3981* 

MACHINE  CENTER  TO  BACK(X^;):  0.00 

SLIDING  BASECXg);  0.00 

RATIO  OF  ROLL(R^):  1.032397 

BLANK  OFFSET (E^) :  0.0 

MACHINE  ROOT  ANGLE (7„) :  70.65“ 


PINION  MACHINE  TOOL  SETTINGS 


CONVEX 

CONCAVE 

CUTTER  BLADE  ANGLE: 

21.5® 

18.5* 

CUTTER  POINT  RADIUS: 

80.4876 

71.7222 

RADIAL  SETTING (s): 

71.55166 

69.04316 

CRADLE  ANGLE (q) : 

59.4638* 

64.0624* 

MACHINE  CENTER  TO  BACK(Xg) : 

1.08497 

-1.58960 

SLIDING  BASE (Xg) ;  -0.25021  0.36659 

RATIO  OF  ROLL (R^) ;  3.898097  3.788604 

BLANK  OFFSET (E_):  -2. 56862 (Up)  2 . 19033 (Down) 

n 

MACHINE  ROOT  ANGLE (7^) : 


13.3333* 


13.3333* 
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Cycle  of  Meshing 


Fig.  1.2.3  Piecewise  Linear  Function  of  Transmission  Errors 
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Fig.  2.3.2  Gear  Generation:  Instantaneous  Axis  of  Rotation 
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^p2 


Fig.  3.1.1  Gear  Head-Cutter 
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Fig.6.3.2  Visualization  of  Orientation  of  Vectors  in  Plane 
Tangent  to 
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AT 


f 


Fig.  6.4.1  Parabolic  Function  of  Transmission  Errors 
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Fig.  6.5.1  Coordinate  Systems  Used  for  Visualization  of  Contact  Path 
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Vc 


Fig.  6.5.3  Coordinates  of  Points  of  Contact  M  and  A 
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Fig.  7.1 .1  Location  of  Mean  Contact  Point  of  Shifted  Bearing  Contact 


CAM  PROFILE  AND  ITS  MOTION 
b  ■  5.5,  E  *15,  »  3.6230 


Cam 


I 


0 


-4.® 


X 


4.0 


13.0 


Fig.  A.3.1  Cam  Profiles  and  Guides 


Transmission  Error  in  Meshing  Pesiod 

G«w  CoMTM  Si4«  -  !!••• 


TOOTH  COITTACT  PATTERM 
G«ar  Caavax  Si^ 


Root  Coaa  Gaaaratriz  (aa> 


Fig.  B.l  Contact  Pattern  and  Transmission  Errors  at  Gear  Mean 
Position  of  Convex  Side 
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Transmission  Error  in  Meshing  Period 

G««r  Coavra  SM«  To* 

V  -  H  •  0.0 


TOOTH  CONTACT  PATTERN 
Goor  Coaoos  SM* 

B  .  Too  V/H  -  »  /  0 


|S  -It.O  -l*.0  -U.t  >4.*  *.*  4.0  $.$  li.$  I*.*  It.O 


Root  Coo*  Goaorotrix  (aa) 


B.2  Contact  Pattern  and  Transmission  Errors  at  Gear  Toe  Position 
of  Convex  Side 
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Transmission  Error  in  Meshing  Period 

G««r  CoaTM  Sida  ~  Haal 
V  — §  ;  H  -  } 


TOOTH  CONTACT  PATTERN 
Gaar  Caaaax  Sida 
Haal  V/H  — S  /  3 


(2  -M.«  -u.a  -i.a  -«#  a.a  4.a  •#  ii.a  i».a  la.t 

Raot  Ceaa  Gaaaratrix  (■■> 


Fig.  B.3  Contact  Pattern  and  Transmission  Errors  at  Gear  Heel  Position 
of  Convex  Side 
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Transmission  Error  in  Meshing  Period 

G«ar  Coac«'^«  Si4a  MEAN 
V  -  0:  H  -  0 


TOOTH  CONTACT  PATTERN 
Gear  Coaeaee  Side 
Meaa  V/H  -0/0 


B  •' 


Fig.  B.4  Contact  Pattern  and  Transmission  Errors  at  Gear  Mean 
Position  of  Concave  Side 
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Transmission  Error  in  Meshing  Period 

G««r  CoaMT*  SM«  ^  TOB 
V  —11;  H  -  3 


TOOTH  CONTACT  PATTBKN 


Root  Coa«  Coaaratru  (aial 


Fig.  B.5  Contact  Pattern  and  Transmission  Errors  at  Gear  Toe  Position 
of  Concave  Side 
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Transmission  Error  in  Meshing  Period 

G««r  Coaeav*  Si4«  -  HmI 
V  -  9  ;  H  — 3 


TOOTH  COKTACT  PATTEXM 
Gaar  Coaeava  SMa 
Haal  V/H  -  9  /-3 


B.6  Contact  Pattern  and  Transmission  Errors  at  Gear  Heel 
Position  of  Concave  Side 
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c. . . 

C....  THIS  PROGRAM  IS  TO  DERIVE  THE  MACHINE  TOOL  SETTINGS 
C....  FOR  PINION  GENERATION  &  TEST  THE  RESULTS 

C.  .  . 

IMPLICIT  REALMS (A-H,0-Z) 

REALMS  KS,KQ,KlI,KIII,K2l,K2lI,KFI,KFII,KD,KF,KH,mcd 
REALMS  M11,M12,M13,L11,L12,L13,L14,M21,M22,M23,L21,L22,L23,L2A, 
5<N11,N12,N2I,N22 
real *8  xi (5) , x  (5) , f  (5) 

EXTERNAL  FCN I , FCN2 , FCN , FCNM , FCNR , FCNMR 

DIMENSION  CH(3) , P (3) , ElEF (3) ,ESN(3) ,EQN(3) ,W1VT2(3) ,WV12(3) , 
SW2VT1  (3)  ,EFIH(3)  ,EFIIH(3)  ,RH(3)  .  GNH  (3)  ,  E2IH  (3)  ,E2IIH''3)  ,PI2P(20)  , 

5. EIIH(3)  ,E1IIH(3)  ,EFI(3)  .EFII(3)  .Eli  (3)  ,  El  1 1  (3)  ,GN(3)  ,EFE1(3)  , 

6. ERR(20)  ,xcp(20)  ,ycp(20)  ,AX1  (20)  .AX2(20)  ,ANG1  (20)  ,ANG2(20) 
COMMON/Al/CNST, TNI, TN2,C,FW, GAMMA, xl.rl .mcd 
COMMON/A2/B 1 ,  RGMAl ,  FGMAl ,  PGMAl ,  DIR  ,  DIF , .At)Dl ,  DEDl 
COMMON/ a3/B2 , RGMA2 , FGMA2 , PGMA2 , D2R , D2F , ADD2 , DED2 , WD, CC , D2P 
COMMON/A4/SR2 , Q2 , RC2 , PW2 . XB2 , XG2 , EM2 , GaMA2 , CR2 , ALP2 . PHI2 , PHI2P 
C0MM0N/A5 / SG , XM , YM , ZM , XNM , YNM , ZNM , X2M . Y2M . Z2M . XN2M , YN2M , ZN2M , 

i<XNH2 ,  YNH2 ,  ZNH2 ,  XH2 ,  YH2  ,  ZH2 

C0MM0N/a6/ES(3) ,EQ(3) ,CN(3) .Wl(3) ,W2(3) ,W12(3) .VT1(3) ,VT2(3) , 
SV12(3) ,KS,KQ,KF.KH,EF(3) ,EH(3) ,SIGSF,PI21 
COMMON/A7/SR1 ,Q1 ,Rc£ , PWl , XBl ,XG1 , EMI , GaMAl , CR 1 , ALPl , PHI  1 , PHI  1 P 
C0MM0N/a8/S£,XM1,YM1 ,ZM1 ,XNM1 . YNMl , ZNMl , XIM , Y 1M . ZIM , 

&,XN1M.  YN1M,ZN1M,XNH1  ,  YNHl.ZNHl  .XHl,  YHl  ,ZH1 
COMMON/a9/PHI2PO,OX,OZ.XO,20,RHO.ALP, V.H.CRIT.PCR! I 
COMMON/A10/K1I,K1II.K2I.K2II,DEL,E1IH,E1IIH.E2IH,E2IIH,GNH, 

&<a2P  ,  B2P  ,  TAU IR  .  TAU2R  ,  a2L.  B2L 

COMMON/Al 1/RAM, PSI 1 .C2,D6,E24,F120,CX6.DX24,EX120.RU1 .uELT.RUP, 
SRAl ,CPF,DPF,EPF,FPF 
CNST=DARCOS (- 1 . ODOO) / 1 80 . ODOO 
C 

c. .  . 

C...  INPUT  THE  CONTROL  CODES 

C.  .  . 

c 

C  IF  V  AND  H  CHECK  IS  NOT  DESIRED,  SET  JCN  =  1 

C  DO  NOT  SET  JCN  TO  BE  3 

C 

JCL=2 

C 

C...  FOR  RIGHT  HAND  GEAR  JCH=1,  FOR  LEFT  HAND  GE.AR  JCH  =2 
C 

JCH=1 

C 

C...  FOR  STRAIGHT  BLADE  JCC=1,  FOR  CURVED  BLADE  JCC=2 
C 

JCC=1 

c 

C  TLl  AND  TL2  .ARE  NUMBER  OF  EXTRA  POINT  ON  CONTACT  PATH 
C  WHICH  SHOULD  NOT  BE  LARGER  THAN  2 

C 

TL1=1.0 

TL2-1.0 

C.  .  , 

C...  INPUT  BLANK  DAT/.  OF  GEAR  AND  PINION 

C.  .  . 

TN1=11.0 

TN2=41.0 

C=0.0 
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FW=27.03 

GAMMA=90.0’'‘CNST 

MCD-76.56 

RGMA1=13.3333*CNST 

B1=35.0*CNST 

B2=35.0*CNST 

RGMA2=70.6500*CNST 

FGMA2=76.6667*CNST 

PGMA2=74.9833*CNST 

D2R=  0.0 

D2F=0.0 

ADD2=2.06 

DED2=6.05 

WD=  8.11 

CC=0.81 

DEL=0.00025*25.4 

INPUT  NORMINAL  RADIUS  OF  GEAR  CUTTER  AND  POINT  WIDTH.  BL.^DE  .ANGLE 

RU2=152. 4000/2.0 
PW2=2.79 
ALP2=20.0*CNST 
DC2=2.0-'RU2 


INPUT  THE  SYNTHESIS  CONDITION  PARAMETERS  .AND  PINION  CUTTER  BLADE 
ANGLE,  ALPKFOR  GEAR  CONVEX.  .ALP1>0,  FOR  GEAR  CONCAVE  ALP1<0) 

GE.AR  CONVEX  SIDE 

FI21=-0.0008 
KD=0. 180 
ETAG=-65.0*CNST 
GaMA1=13.3333”CNST 
RHO=  250.0 
C2=  0.00 
ALP1=18.500*CNST 
TN1I=8.0 

SGN=DSIN(aLP1) /DaBS(DSIN(aLP1)) 

KSIDE=0 


GOTO  1989 

**••*  GEAR  CONCAVE  SIDE 

1990  CONTINUE 

FI21=  0.0008 

KD=0. 180 

ETAG=  65.0’'CNST 

GAMA1=13.3333''CNST 

RHO-  200.0 

C2-  0.00 

ALPl  — 21.50’'‘CNST 
TNlI-8.0 

C.  .  . 

SGN-DS IN (ALP  1) /DABS (DS IN (ALP  1 ) ) 

KSIDE=1 

jcl»2 
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c 

c 

C...  INPUT  GEAR  .MACHINE  TOOL  SETTINGS 
C 

cl989  Q2=  52.6589-'CNST 
c  SR2=3.8872"25.^ 

c  XG2=0.0 

c  XB2=-0. 0333*25. 4 

c  CR2=. 9772974 

C  RaG=1.0/CR2 

c  GAMA2=RGMa2 

c  E.M2=0.0 

c  RC2=RU2-SGN*PW2/2.0 

c  ALP2=  SGN*ALP2 

C 
C 

C...  CALCULATE  GEAR  .HaCHINE  TOOL  SETTINGS 
C 

1989  hg=nicd''dcos  (pgma2-rgma2)  -ru2''ds  in  ( b2  ) 
vg=ru2“dcos (b2) 
q2=datan (vg/hg) 
sr2=dsqrt (hg**2*vg**2 ) 
xg2=0.0 
GaKa2=RGMA2 
xb2=d2r "ds in (gama2 ) 

EM2=0.0 

rag=dcos (pgma2-rgma2) / ds in {pgma2) 
cr2=l  .0/rag 
RC2=RU2-SGN*PW2/2.0 
ALP2=  SGN*ALP2 
C 

C...  DELT  IS  THE  CA.M  SETING 
C 

DELT=0.0 

C 

C...  DEFINE  THE  MEAN  CONTACT  POINT 
C 

V=0.00G 

H=0.000 

FA=FGMA2-PG?1A2 

RA=PGMA2-RGMa2 

HH=CC*WD-0 . 5*FW* (DTAN  < FA) -DTAN (RA) ) 

DED2R  =  DED2-0 . 5*FW*DT.AN  (RA) 

XL=HCD*DCOS  (PGMA2)  -  (DED2R-H.M/2 . 0)  *DSIN  (PGMA2) 
RL=MCD*DSIN  (PG1U2)  -  (DED2R-H.M/2 . 0)  *DCOS  (PGflA2) 

C.  .  . 

AGL=DATAN(RL/XL) 

0X=-DSQRT  (XL**2*RL**2)  *DCOS  (AGL-RGM.A2) 
0Y=-D2R'-DSIN(RGMA2) 

C  WRITE(9.11)  OX.OY.XL.RL 
C 
C 
C 

C...  FIND  SURFACE  COORDINATES  OF  THE  MEAN  CONTACT  POINT 
C 

ERRREL=0. lD-10 
N=2 

ITMAX=200 

IF  (JCH.EQ.l)  THEN 
Q2='Q2 
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XI(1)=270.0*CNST+B2 

ELSE 

XI  (1)='B2 
END  IF 
XI (2)=0.0 

CALL  DNEQNF (FCNl , ERRREL , N . UMAX . XI , X . FNORM) 

TH=X(1) 

PH=X(2) 

ST=DSIN(TH) 

CT=DCOS(TH) 

SH=DSIN(PH) 

CS=DCOS(PH) 

SP=DSIN(ALP2) 

CP=DCOS (ALP2) 

SM=DSIN (GAMA2) 

CM=DCOS (GAMA2) 

THIG=TH 

C  WRITE(9,11)  Xn2M,Yn2M.Zn2M.ZNM.YNM,2NM 

C  WRITE(9,11)  XM,  YM.ZM.sg.hrn 

C 

C...  DEFINE  VECTORS  TO  COMPUTER  THE  SECOND  ORDER  PROPERTY  OF  GEAR 
C 

ES(1)=-DSIN(TH-PH) 

ES(2)=  DCOS(TH-PH) 

ES(3)=  0.0 

EQ  (1 )  =-SP'"DCOS  (TH-PH) 

EQ(2)=-SP''DSIN(TH-PH) 

Eg(3)=-CP 

CN(1)=XNM 

CN(2)=YNM 

CN(3)=ZNM 

KS=CP/(RC2-SG*SP) 

KQ=0.0 
W1  (1)=-CM 
Wl(2)=  0.0 
W1  (3)=-SM 
W2(l)=  0.0 
W2(2)=  0.0 
W2(3)=-CR2 

VT1(1)=  YM*SM*EM2*SM 
VTl (2) =-XM*SM* (ZM-XB2) "CM 
VTl  (3)=-YM--'CM-EM2"C.M 
VT2(1)=  YM--'CR2 
VT2(2)=-XM*CR2 
VT2(3)=  0.0 
DO  10  I»l,3 
W12(I)-W1(I)-W2(I) 

V12(I)-VT1 (I)-VT2(I) 

10  CONTINUE 
C 

C...  FIND  THE  PRINCIPAL  DIRECTION  A.ND  CURVATURES  AT  MEAN  POINT 
C 

PI21=0.0 
CALL  CURVAI 

C  WRITE(9,12)  KF.KH.SIGSF 
C  12  F0RMAT(3X,3(G1A.7,2X)) 

K2I=KF 
K2II=KH 
PHI2=PH/CR2 
sh2=ds in (phi  2) 
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ch2*dcos (phi2) 

xX-  CM’'ef  (l)+SM*ef  (3) 

yY-  ef(2) 

zZ=-SM*ef (l)+CM*ef (3) 
ef (1) =xx 

ef(2)-  CH2*yY-SH2*zZ 
ef(3)=  SH2*yY+CH2’'zZ 
C  WRITE(9,11)  xx.yy.zz 
c 

xX=  CM*eh(l)+SM=-eh(3) 
yY=  eh(2) 

zZ=-SM*eh(l)+CM*eh(3) 

eh (1) »xx 

eh(2)«  CH2’-yY-SH2*z2 
eh  (3)*  SH2*yY+CH2*zZ 
C  WRITE(9,11)  Ef  (1) ,Ef  (2) ,Ef (3) 

C  WRITE(9,11)  Eh(l) ,Eh(2) ,Eh(3) 

C 

C 

C 

ERRREL”0. lD-10 
N=1 

ITMAX=200 

XI  (i)=o.o 

CALL  DNEQNF(FCN2,ERRREL.N. IT^UX.XI.X.FNORM) 
PHI2P0=X(1) 

C  WRITE(9.11)  X(1) 

C  WRITE (9, 11)  XH2.YH2,ZK2 

C  WRITE(9.11)  XNH2, YNH2.ZNH2 

C 

CHP-DCOS(X(l)) 

SHP=DSIN(X(1)) 

CMiM=DCOS  (gamma) 

SMM=DS IN (GAMMA) 

XX=  ef(l) 

YY=-e£(2)-'CHP+ef  (3)*shp 
ZZ— ef  (2) ’'SHP-ef  (3)*chp 
EF(1)“  XX''CMM*ZZ’'SMM 
ef(2)-  YY 

EF  (3)  — XX-‘SMM*ZZ*CMM 

C.  .  . 

XX=  eh(l) 

YY— eh(2)*CHP‘eh(3)*shp 
ZZ— eh  (2)  *SHP-eh  (3)  *chp 
EH(1)-  XX*CMM+ZZ*SMM 
eh(2)-  YY 

EH  (3) --XX*SMM+2Z*CMM 
C  WRITE(9,11)  EF(1) ,EF(2) .EF(3) 

C  WRITE(9,11)  EH(1) ,EH(2) ,EH(3) 

ETAG“90 . 0*CNST>S IGSF+ETAG 
C 

C  LOCAL  SYNTHESIS  AT  MEAN  CONTACT  POINT 

C 

C.  .  . 

RH(1)-XH2 

RH(2)-YH2 

RH(3)-2H2 

GNH(1)-XNH2 

GNH(2)-YNH2 

GNH(3)-ZNH2 
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E2IH(1)=  EF(1) 

E2IH(2)=  EF(2) 

E2IH(3)=  EF(3) 

E2IIH(1)=  EH(1) 

E2IIH(2)=  EH(2) 

E2IIH(3)=  EH(3) 

K2I=KF 
K2II=KH 

...  RELATIVE  MOTION  P.ARAMETERS  IN  GEAR  i.  PINION  MESHING  PROCESS 

R12=TN1/TN2 
W1(1)=~1.0D00 
WI (2)=0.0D00 
W1(3)=0.0D00 
W2(l)=  R12’'CHM 
W2(2)=0.0D00 
W2(3)=-RI2"SMM 
W12(U=W1  (1)-W2(1) 

W12(2)=W1 (2)-W2(2) 

W12(3)>W1 (3)-W2(3) 

VT1(1)=  O.ODOO 
VTl (2)=  ZH2 
VTl (3)=-YH2 

VT2(1)=  R12*(YH2-C) -SMM 
VT2(2)=-R12*  (XH2''SMM-ZH2*CMM) 

VT2(3I=  R12-'(YH2-C)  •■•CMM 
V12(l)=  VT1(D-VT2(1) 

V12(2}=  VTl  (2)-VT2(2) 

V12(3)=  VTl (3)-VT2(3) 

WRITE(9,3)  V12(1) ,V12(2^ . V12(3) 

F0RM.AT(5X,3G1^.7,/) 

CALCULATE  THE  COEFFICIENT  a13.a23.A33 

ESN(1)=  GNH(2)  "E2IH(3)-GNH(?I  •'E2IH(2) 
ESN(2)=-(GNH(1)*E2IH(3)-GNH(3)  ■''E2IH(1)) 

ESN  (3)  =  GNH  ( 1 )  ■••E2IH  (.2)  -GNH  (2)  ••■•E2IH  ( 1 ) 

C.  .  . 

EQN(1)=  GNH (2) -E2 I IH (3) -GNH (3) *E2I IH (2) 
EQN(2)=-(GNH(1)*E2IIH(3)-GNH(3)*E2IIH(1)) 

EQN(3)=  GNH(l)  -E2IIH(2)-GNH(2)  •■•E2IIH(1) 

C.  .  . 

W1VT2(1)=  Wl  (2)*VT2(3)-W1  (3')*VT2(2) 

W1VT2  (2)  =-  (Wl  (1)  *VT2  (3)-Wl  (3)  ■•■'VT2  (1) ) 

W1VT2(3)=  Wl (1)*VT2(2)-W1 (2)*VT2(1) 

C. .  . 

W2VT1  (1)=  W2(2)*VT1 (3)-W2(3)*VTl (2) 

W2VT1  (2)=-(W2(l)*VTl  (3)-W2(3)*VTl  (D) 

W2VT1(3)-  W2(1)*VT1  f2')-W2(2)*VTl(l) 

C. .  . 

WV12(1)=  W12(2)*V12(3)-W12(3)*V12(2) 
WV12(2)=-(W12(1)*V12(3)-W12(3)*V12(1)) 

WV12(3)=  W12(1)*V12(2)-W12(2)*V12(1) 

C. . . 

VI2S=0.0D00 
V12Q=O.ODOO 
WNES=0.0D00 
WNEQ=O.ODOO 
VWN=  O.ODOO 
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WITN-O.ODOO 

W2TN-O.ODOO 

VT2N=0.0D00 

C.  .  . 

DO  1  I-l,3 

V12S=  V12(I)*E2IH(I)*V12S 
V12Q=  V12(I)*E2IIH(I)+V12Q 
WNES=  W12(I) -'ESNCD-^WNES 
WNEQ-  W12(I) -'EQN(I)-WNEO 
VWN  -GNH(I)*WV12(I)-^VWN 
W1TN=GNH(I)*W1VT2(I)+W1TN 
W2TN=GNH(I) *W2VT1 (I)+W2TN 
VT2N=  GNH(I)*VT2(I)+VT2N 
1  CONTINUE 

WRITE(9,6)  V12S.V12Q 
6  F0RMAT(5X.2G14.7,/) 

COMPUTER  THE  COEFFICIENTS  a13,a23,a33 

B13=-K2r-V12S-WNES 
B23=-K2II*V12Q-WNEQ 

B33=K2I*V12S-"-2-^K2II*V12Q*"2-VWN-W1TN-W2TN*VT2N’TI21*TN2/TM 

..  LOCAL  SYNTHESIS  OF  MESHING  AT  MEAN  CONTACT  POINT 

C.  .  . 

DL  =  KD-'FW 
SIGK2=  K2I+K2II 
SIGG2=  K2I-K2II 
A-DEL/DL-'*2 

T 1 = -  B  1 3  ■' V 1 2Q+  (  B  3  3  -  B  1 3  •••  V 1 2  S )  *DT AN  (  E  T AG ) 

T2=B33+B23-‘(V12Q-V12S -■DT.AN  (  ETaG ) 

ETAP-DATAN(T1/T2) 

VSl  ’B33/  (B13-B23’-DTAN<ETAP)  ) 

AM1=DSIN(ETAP)"--2 

AM2=-DS IN (2 . ODOO*ETAP) II.  ODOO 

AN1=  (B  13-B23"DTAN  (ETAP'  'i  /  ( (1 .  ODOO+DT.AN  ( ETaP)  •”•2)  ''VS  1 ) 

AN2=  (B  13*DTAN  (ETAP''  -B23I  /  ( ( 1 .  ODOO^DT.AN  (ETaP  )  -  2 )  -'VS  1 ) 
SGN=DSIN(ALP1) /DABS (DS IN (aLPI ) ) 

A=A*SGN 

SIGK=(4.0D00*A**2-  (AN1-"'2*AN2*"2)  )  /  (  2.  ODOO'-'A-  (.ANi^DCOS  12.  ODOO 
&ETAP) +AN2*DSIN (2.0D00*ETAP) ) ) 

SIGK1=  SIGK2-SIGK 

Tl=  2. ODOO  ■'■aN2-SIGK*DSIN(2. ODOO*  ETAP) 
T2=SIGG2-2.0DOO*AN1*SIGK*DCOS(2.0DOO*ETAP) 

S I G 1 2  = .  5 ODOO  *DAT.AN  ( T 1  / T2 ) 

SIGGX2.0DOO*AN2-SIGK*DSIN(2.0DOO*ETAP)  )  /DSIN(2.0D00*SIG12) 

K1I=(SIGK1+SIGG1)/2.0D00 

K1II=(SIGK1-SIGG1)/2.0D00 

C.  . 

C.  . 

C.  . 

C  WRITE(9.11)  ETAP. KlI. Kill 
C  WRITE(9,in  SIGK,SIG12,SIGK1,SIGG1 
C  WRITE(9,8)  T1,T2 
C  8  F0RMAT(5X,3G14.7) 

c. . 
c. . 
c. .  . 

C...  PRINCIPLE  DIRECTIONS  OF  PINION  SURFACE  AT  POINT  M 

C. .  . 
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DO  15  1=1,3 

E1IH(I)=  DC0S(SIG12)*E2IH(I)-DSIN(SIG12) "E2IIH(I) 

E1IIH(I)=  DSIN(SIG12)  •'E2IH(I)-^DCOS(SIG12)*E2IIH(I) 

15  CONTINUE 

C  WRITE(9,11)  ElIH(l) ,E1IH(2) ,E1IH(3) 

C  WRITE(9,11)  ElIIH(l)  ,EinH(2)  ,E1IIH(3) 

C.  .  . 

C...  COINCIDE  THE  NORMALS  OF  CUTTER  AND  THE  PINION  SURFACES 

C.  .  . 

SMl=DSIN(GAMAl) 

CMI=DCOS(GAMAl) 

SP=DSIN(ALP1) 

CP=DCOSCaLP1) 

T1=-(XNH2+SP*SM1) 

T2=  CP*CM1 

IF  (JCH.EO.l)  THEN 

THF=DARCOS(Tl/T2) 

ELSE 

THF=DARCOS(Tl/T2) 

THF=360.0*CNST-THF 
END  IF 

Ba1=-CP*DSIN(THF) 

BA2=  CP'-'DSIN  (GAMaD  *DC0S  (THFl  -SP*DCOS  (GAMAD 
TT=-  (YNH2---*2*ZNH2'-*2) 

CSH=- (Ba1*YNH2+Ba2"2NH2) /TT 
SNH=(BA2*YNH2-BAr-2NH2) /TT 
PH  IH=2 . 0*DATAN2  (SNH .  ( 1 .  ODOO-^CSH  ■  :> 

C  WRITE(9,8)  THF,PHIH 

C 

C...  FIND  THE  PRINCIPAL  DIRECTIONS  OF  PINION  GENERATING  SURFACE 
C 

EFI (1)=-DSIN(THF) 

EFI(2)=  DCOS(THF) 

EFI(3)=  O.ODOO 
EFII(1)=  SP*DCOS(THF) 

EFII(2)=  SP''DSIN(THF) 

EFII(3)=-CP 

C  WRITE(9,11)  EFI  (1) ,EFI(2) .EFI (3) 

C  WRITE(9,il)  EFII (1) .EFII (2) ,EFII (3) 

C 

C...  FIND  THE  PINION  PRINCIPAL  DIRECTIONS  IN  SYSTEM  SMl 
C 

XX=  ElIH(l) 

YY=  DCOS(PHIH)*ElIH(2)-*-DSIN(PHIH)*ElIH(3) 

22=-DSIN(PHIH)  •'E1IH(2)-DC0S(PHIH)  -'ElIHCS) 

E1I(1)=  CM1’'XX-SM1*22 
E1I(2)=  YY 

E1I(3)=  SMl*XX+CMr-22 
C 

XX=  ElIIH(l) 

YY=  DCOS(PHIH) ''E1IIH(2)-DSIN(PHIH)*E1IIH(3) 

22=-DSIN  (PHIH)  •'Ell  IH  (2)  -^DCOS  (PHIH)  "Ell  IH  (3) 

ElII  (1)=  CM1*XX-SM1*22 
E1II(2)=  YY 

El  1 1  (3)  =  SMl*XX-'-CMl*2Z 
C 

C...  FIND  THE  UNIT  NORNM.^L  IN  SYSTEM  SMI 
C 

XX-  XNH2 

YY-  DC0S(PHIH)*YNH2*DSIN(PHIH)*2NH2 

I4q 
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ZZ— DSIN  (PHIH)  *YNH2+DC0S  (PHIH)  *ZNH2 
GN(1)-  CMl*XX-SMr^ZZ 
GN(2)»  YY 

GN(3)=  SMP'XX-t-CMP'ZZ 
C 

C...  EXPRESS  THE  POSITION  VECTOR  IN  SMI 
C 

XX*  XH2 

YY=  DCOS (PHIH) *YH2+DSIN (PHIH) "2H2 
2Z=-DSIN (PHIH) *YH2-DCOS (PHIH) "2H2 
RX=  CM1*XX-SM1*Z2 
RY*  YY 

RZ=  SM1*XX+CM1*Z2 
C  XX*-CP’''DCOS(THF) 

C  YY=-CP"DSIN(THF) 

C  ZZ=-SP 

C.  .  . 

C  WRITE(9,11)  Ell (1) .Ell U) .E1I(3) 

C  WRITE(9,11)  Ellld;  .EIII'.E;  .E1II'3; 

C  WRITE(9,11)  GN(1) ,GN(2) ,GN(3) 

C  WRITE(9,11)  xx.yy.zz 

C  WRITE(9,11)  RX.RY.RZ 

C.  .  . 

DO  20  1=1,3 
XX=-EFI (I) 

EFI(I)=  EFIKI) 

EFII (I)=XX 
20  CONTINUE 

WRITE  (9. 11)  Ell  (1)  ,  El  U2)  .Eli  (3; 

WRITE  (9,  11)  El  1 1  (1)  ,  El  1 1  (2)  ,£2  1 1  .'3; 

..  FIND  THE  ANGLE  FORMED  BETWEEN  PRINCIPAL  CURVATURES 

EFEl  (1) =  EFI (2) "Ell (3) -EFI (2) "El  I  (2) 

EFEl (2)=-EFI (1)"E1I (3)-EFI (3)*E1I  (1) 

EFEl  (3)=  EFI  (1)*E1I  (2)-EFI  (2)  "ElKD 
T1=0.0D00 
T2=0.0D00 
DO  30  1=1,3 
T2=EFI(I)*E1I(I)-T2 
T1=GN(I) *EFE1 (I)*T1 
30  CONTINUE 

SIGF1=2.0*DATAN2(T1, 1.0*T2) 

..  FIND  THE  CURVATURE  OF  PINION  GENERATION  SURFACE  AT  MEAN  POINT 

IF  (JCC.EQ.l)  THEN 
KFI=0.0 

B 1 2=0 . 5D00*  (K1 1-Kl 1 1) *DS IN (-2 . ODOO" S IGF  1 ) 
B11=KFI-K1I*DC0S(SIGF1)**2-Klir-DSIN(SIGF1) "*2 
TKK=  K 1 1 *DS I N (S I GF 1 ) **2-Kl 1 1 "DCOS (S IGF  1 ) **2 
KF 1 1  =  ( B 1 2 *  2 + B 1 1  *TKK)  / B  1 1 
WRITE(9,11)  SIGF1,KFII 

. .  FIND  THE  CUTTER  POINT  RADIUS  AND  ITS  CENTER 

SF=SG-'DCOS  (ALP2)  /DCOS  (ALPl) 

SF'DABS (RZ) /DCOS (ALPl) 

RCF=CP/DABS(KFII)-SF*SP 
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DO  AO  1-1,3 

P(I)-  GN(I)*CP-EFI (I)*SP 
AO  CONTINUE 

RCX-  RX-SF*EFI (1)+RCF*P(1) 

RCY-  RY-SF*EFI  (2)+RCF-'P(2) 

RCZ-  R2-SF*EFI  (3)-^RCF*P(3) 

C  WRITE(9,11)  RCF.SF 

C  WRITE(9,11)  RCX, RCY, RCZ 

ELSE 

KFI-l.O/RHO 

B12=0.5DOO*(K1I-K1II)*DSIN(-2.0DOO*SIGF1) 

B 1 1 -KF I -K1 1 *DCOS ( S I GF 1 ) **2-K 1 1 1 *DS I N ( S I GF 1 ) **2 
TKK=  K1 1-'DSIN  (SIGFl)  **2'^Klir'DCOS  (SIGFl)  **2 
KFII-(B12’''*2+B11*TKK)/B11 
C  WRITE(9,11)  SIGFl.KFII 

DBT—RZ 

RM=CP/DABS(KFII) 

ZO=-(DBT+RHO-'SP) 

XO=RM-RHO*CP 

RCF-XO+RHO'-DSQRT  (1.0-  (ZO/RHO)  **2) 

RCX=RHO*GN ( 1) -XO*DCOS (THF) -RX 
RCY=RHO=-GN  (2)  -XO*DSIN  (THF)  -RY 
RCZ=RHO*GN(3)-ZO+RZ 
C  WRITE(9,777)  XO.ZO 

C777  FOR.HAT(3X,'  XO,  XO  =  '  .  2  (2X,  Gl  A .  7)  ) 

C  WRITE(9,11)  RCF.RM.DBT 

C  WRITE(9,11)  RCX. RCY, RCZ 

END  IF 
C 

c 

C...  THE  FOLLOWING  IS  TO  FIND  THE  CUTTING  RATIO 
C 

CSMl-CMl 

SNMl-SMl 

T1X=EFI(1) 

TIY-EFI (2) 

TIZ-EFI (3) 

T2X-EFII  (1) 

T2Y-EFII  (2) 

T2Z-EFII (3) 

XN=GN(1) 

YN=GN(2) 

ZN=GN(3) 

RXC-RX 

RYC-RY 

RZC-RZ 

C 

C  WRITE(9,11)  RCX, RCY. RCZ 

C  WRITE(9,11)  TIX.TIY.TIZ 

C  WRITE(9,11)  T2X,T2Y,T2Z 

C  WRITE(9,11)  XN.YN.ZN 

C  WRITE(9,11)  RXC.RYC.RZC 

C.  . . 

C...  THE  FOLLOWING  IS  TO  DETERMINE  DELTA. EM. AND  IFM 

C.  . . 

M11-XN*T1Y-YN*T1X 
M12— CSM1*(YN*T1Z-ZN*T1Y) 

M21-XN*T2Y-YN*T2X 
M22— CSMl*  (YN*T2Z-ZN*T2Y) 

C  WRITE(9,11)  M11,M12,M21,M22 
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L11=(B12/B11*M11-M21)/KFII 

L12=(B12/B11*M12-M22)/KFII 

L21=-T2Z/T1Z*L11 
L22=-T2Z/T1Z-'L12-RYC*CSM1/T1Z 
DTT=B  12*KF  I  *T2Z-^B  1 1  *KFI  I*T1Z 
L11=-T1Z*(B11*M21-B12*M11) /DTT 

L12=(-B12*KFT-'RYC*CSM1-T1Z"(B11*M22-B12*M12)  ')  /DTT 
L  2 1 = T  2  Z  *  ( B 1 1  *  M  2 1  -  B 1 2  *  M 1 1)  /  DT  T 

L22=  ( -B 1 1  *KF  1 1  *R  YC*C SM 1  T22  -  (B 1 1  *M22- B 1 2*M  1 2 )  )  / DTT 

WRITE(9,11)  L11,L12,L21,L22 

X11  =  L21*T1X+L1P''T2X 

X12=L22*T1X+L12*T2X 

X21=L21-'T1Y+L11*T2Y 

X22=L22*T1Y+L12*T2Y 

X31  =  L21’'T1Z+L11*T2Z 

X32=L22*T1Z+L12*T2Z 

E11=YN*X11-XN*X21 

E12=YN-'X12-XN-'X22-ZN-X2r-CSM1-YN'''X31*CSM1 

E13=-(ZN*X22-YN*X32) "CSMl 

Y1 1=-XN* (RXC*SNM1-RZC"CSM1 ) -YN*RYC-'SNM1 

Y12=Y11"SNM1 

X13=X12-RYC*SNM1 

X23=X22+RXC'’SNM1-RZC*CS.M1 

X33=X32+RYC'-CSM1 

Y2 1=-XN*X2  I  -'SNM1^ZN'*X21  "CSMl- YN-  (XI  i  -SNMl  -X3 1  '■CSMi  ' 
Y22=-XN*X23*SNM1+YN'PSNM1"X13-CSM1*X33>-ZN*X23"CSM1 


THE  EFFECT  OF  SECOND  ORDER  RATIO  OF  ROLL  ON  A33 
T.Ml  =XN  ■■'X 1 1  -  YN'"'X2 1  -ZN"X3  i 
TM2=XN*X13*YN"'X23-ZN"X33 
C  WRITE (9, 163)  TM1,TM2 

C163  FORMAT (2X, 'TM1.TM2  ’ . 2  (2X . G14 . 7 ) ) 

221=C2*TM1 

ZZ2=C2*  (TM2-^SNMl*TMr) 

ZZ3=C2*SNM1*TM2 

C.  .  . 

C  Zl=KFir-Lll**2-Ell 

C  Z2=2.0D00*KFir''Lir'Ll2-EI2-Y21-YIl 

C  Z3=KFir--L12**2-E13-Y22-Y12 

Z I  =  K  F  r- •  L  2 1  *  -  ■  2 + K  F I  r  •  L 1 1  •  •  •  2  -  E 1 1  -  Z  Z I 

Z2=2.0*KFr'L21'-L22+2.ODOO-KFII"Lll*Ll2-E12-Y21-Yl 1* 

Z3=KFr'''L22**2+KFir-L12-*2-EI3-Y22‘Y'  2  +  ZZ3 

N11=KFII*L1I+M21 

NI2=KFH*L12+M22 

N21=KFr-L2UMll 

N22=KFI*L22*M12 

C  WRITE(9,11)  B12,Z1 .Mil .Nil 

AA=B12*Z1-N21*N11 
BB=B12*Z2-N21*N12-N22*N11 
CCC=B12*Z3-N22*N12 
^  WRITE  (9.1!)  AA.BB.CCC 

IF  (AA.GT. 0.000001)  GOTO  1949 

T1=-CCC/BB 

GOTO  1950 

1949  T1=(-BB+DSQRT(BB'-"2-4.0DOO"AA*CCC)  )  /  (2 . 0D00=''.A.-\) 

1950  FM1=T1+SNM1 
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CR1=FM1 

RAP=1.0/CR1 

VF3=X3r-Tl+X32+RYC*CSMl 

C.  .  . 

C.  .  . 

C.  .  . 

C...  THE  DETERMINATION  OF  EM  AND  DELTA 

C.  .  . 

EMI=(Xir-TI-^RYC*FMI-XI3)/FMl 
XGI=(X2I*T1-RXC*FM1-X23) / (FMI*CSM1) 

C.  .  . 

RCX=RCX--XGI*CSMI 

RCY=RCY-EM1 

RC2=RC2^XG1*SNMI 

C.  .  . 

VI=RCY 

H1=RCX 

XBI=RCZ 

SR  1=DS0RT  (  V 1  1  •'*2; 

Q1=-DARSINIV1/SR1) 

XB1=-XBI 

....  DETERMINE  THE  CAM  SETTING 

RAl=l .0/CRl 

R.AM=TNI/TNII*RAI 

PS  I  I'DATAN  (C2--'RAM/  ( R.\M-  1 .  OJ  ) 

RL'P=  15 . 0*DCOS  (PSI 1 )  /  (RAM-  1 . 0) 

RUI=RUP 

C  DELT=0.0 

C 

CALL  CAM 

C  WRITE(9,191)  RUP 

C  WRITE(9,191)  RUl 

C191  FORMAT (2X, 'RUl  DELT  =  ’ , 2 (2X.G1A. 7) ) 

C  WRITE(9,199)  RAl ,C2.D6,E24,FI20 

C  WRITE(9.199)  RAl , CPF , DPF . EPF . FPF 

C  199  F0RMAT(2X. 'C.AM' ,2X,5(2X,G1A.7),) 

C 

c. .  . 

c 

c  WRITE(9,25)  FMl 

C  25  F0RMAT(5X.’  FMl  =  ',G1A.7) 

C  WRITE(9,A4)  EM1,SR1,Q1 

C4A  F0RMAT(2X.  'EMl.SRl.Ql  '  ,  3  (3X,  G14 . 7)  :i 

C  WRITE(9,A5)  XGl ,XB1 , VI ,H1 

CA5  FORMAT (2X, ' XGl , XB 1 , V 1 .HI' , A (2X.G14. 7) ) 

C.  .  . 

IF  (KSIDE.EQ.0.0)  THEN 
WRITE(9, 131) 


$  2X,'*  OUTPUT  FOR  GEAR  CONVEX  SIDE  *'  J 

ELSE 

WRITE(9,331) 

331  FORMAT  (/2X,  '  **********************Vt******A***Vr*V:***vVt  .  _  ! 
"  2X,'*  OUTPUT  FOR  GEAR  CO.NCAVE  SIDE  / 

END  IF 
WRITE(9, 13) 
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FORMAT  (/2X,  '  *•'*******'*'•'************'******************  ’  j 
$  2X,’*  gear  CUTTER  SPICIFICATIONS  *’,/ 

&  2  X  '  ’’  *  *  *  *  *  *  *  *  *  *  *  *  *  *  ’**  *  *  *  *  *  ■*  *  *  *  iK  it  it  it  it  it  it  i;  it  it  it  it  f  y  ^ 

WRITE(9.115)  DC2,PW2,ALP2 

FORMAT(/2X,'  GEAR  CUTTER  DIAMETER  :  DC2  ='  ,G14.7,/ 

$  2X, '  CUTTER  POINT  WIDTH  :  PW  ='  ,G14.7,/ 

&  2X,  '  CUTTER  BLADE  ANGLE  :  PHI2  ='  .G14.7,//) 

WRITE(9,3) 

FORMAT (/2X,  '******-'-*^''********************************’  / 
$  2X,'*  BASIC  GEAR  MACHINE-TOOL  SETTINGS  ■”  , / 

^  2X, ’ ****************************-**********-- . _ /) 

WRITE  (9 ,4)  Q2 , SR2 , XG2 , XB2 , EM2 , GAMA2 , RAG 


FORMAT (/2X, ' 

Basic  cradle  angle  : 

Q2  = 

’  ,G14.7./ 

& 

2X,  ' 

RADIAL  SETTING  : 

SR2  - 

'  .G14.7,/ 

// 

2X,  ' 

machine  center  to  back 

:  XG2  =• 

'  ,G14.7./ 

// 

2X,  ' 

SLIDING  base  : 

XB2  = 

'  .G14.7,/ 

S 

2X,  ' 

BLANK  OFFSET  : 

EM2  = 

'  ,G14.7,/ 

// 

2X,  ' 

■MACHINE  ROOT  ANGLE  : 

GAMA2  = 

'  .G14.7,/ 

// 

2X.  ’ 

RATIO  OF  ROLL  : 

RAG  = 

'  ,G14.7,//) 

WRITE(9 

6) 

FORMAT(//2X, 

■  .r*********,c**********Vt 

& 

2X. 

BASIC  PINION  MACHINE-TOOL 

SETTINGS 

WRITE(9.7)  ALPl,  RCF.  Q1 . SR  1 , XGl , XB 1 . EMI , GAMAl , RAP 


FORMAT (/2X, ' 

BLADE  angle  : 

ALPl 

=  ’ 

.G14.7,/ 

& 

2X,  ' 

POINT  RADIUS  : 

RCF 

=  ’ 

,G14.7,/ 

& 

2X,  ' 

Basic  cradle  angle  : 

Q1 

=  ’ 

,G14.7./ 

S 

2X,  • 

R AD lAL  ScTiiNG  • 

SR  I 

=  ■ 

,  G 1 A  .  7  , 

L 

2X.  ' 

-MACHINE  CENTER  TO  BACK 

;  XGl 

=  ' 

,GU.7./ 

S 

2X,  ' 

SLIDING  Base  : 

XB2 

=  ' 

.Gi4. 7,/ 

& 

2X,  ' 

BLANK  OFFSET  : 

EMI 

=  ■ 

.G14.7,/ 

S 

2X,  ' 

MACHINE  ROOT  /„NGL£  : 

GA>LAi 

=  ■ 

,G1A. 7, / 

& 

2X,  ’ 

R.ATIO  OF  ROLL  : 

RAP 

=  ' 

,GU.7,/ 

IF  (JCC.EQ.2}  THEN 
WRITE  (9, 61) 

F0RMAT(//2X.  '^>******-*^''-':"*^= 

£.  2X,'”  COORDINATES  OF  THE  CENTER  O'"  IHE  .ARC  ' 

WRITE(9,71j  XO,ZO 

FORMAT  (/2X,'  RADIAL  COORDINATE  :  XO  =  '  .GU.',. 

&  2X,  '  AXIAL  COORDINATE  :  ZO  = '  ,G14.7,//) 

ELSE 

GOTO  1919 
END  IF 


1919  CONTINUE 

WRITE(9, 16) 

16  FORMAT  (//2X,  ' 

&  2X,’-CAM  SETTINGS  AND  COEFFICIENTS  OF  TAYLOR  SERIE 

WRITE(9,17)  PSIl,  RUP,  DELT,  RAl , C2 . D6 . E24 . F 1 20 


p 

0RMAT(/2X.  ' 

GUIDE  ANGLE  : 

PSIl  =  ' 

,G1A.7,/ 

& 

2X,  ' 

CAM 

PITCH  RADIUS  : 

RUP  =' 

,  G I  A  .  ,  ^  / 

& 

2X,  ' 

CAM 

SETTING  : 

DELT  =' 

,  G  *  A  .  7  , 

& 

2X,  ' 

1ST 

ORDER  COEFFICIENT 

RAl  =' 

,  G 1  A  .  7  ,  / 

& 

2X,  ' 

2ND 

ORDER  COEFFICIENT 

C2  =  ' 

.  G 1  A  .  7  ,  z 

& 

2X,  ' 

3RD 

ORDER  COEFFICIENT 

D6  =  ' 

,  G  i  ^ 

i. 

2X.  ' 

4TH 

ORDER  COEFFICIENT 

E  2  A  =  ’ 

,  G  1  A  .  7  .  ■ 

i, 

2X.  ' 

5TH 

ORDER  COEFFICIENT 

F120 

.  G ;  -  .  7  . 

C 

c. .  . 


CALL  TCA 
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c 

c. , . 

...  DEFINE  THE  INITIAL  POINT 


XI  (D-THIG 
XI  (2) =  0.000000 
XI (3)=THF 
XI(A)=0.0 
XI(5)=  0.00 

FIND  THE  INITIAL  CONTACT  POINT 

5555  N=5 

ERRREL=0. lD-10 
ITMAX=200 
PHI2P=PHI2P0 
IF  (JCC.EQ.l)  THEN 

CALL  DNEQNF(FCN,ERRREL.N,ITMAX.XI.X,FN0RM; 

ELSE 

CALL  DNEQNF(FCNR.ERRREL.N. ITMAX.XI . X, FNORM) 

END  IF 
PHI1P0=X(5) 

C 

C 

C 

c 

?k:2P1  =  PHI2P0-180.0"CN'ST/TN2-TL1''180.0-X:oT/  ;6.0*TN’2) 
PHI2P2  =  PHI2P0+180.0*CNST/TN2'TL2*180.0*CNST/  (6.0'-'TN2) 

r...<=  1 

PHI2P=PHI2P1 
333  CONTINUE 

IF  (JCC.EQ.l)  THEN 

CALL  DNEONF(FCN,ERRREL.N. ITMAX.XI. X, FNORM) 

ELSE 

CALL  DNEQNF(FCNR,ERRREL.N,IT.MAX.  XI.  X.  FNORM) 

END  IF 

C  XI(1)=X(1) 

C  XI(2)=X(2) 

C  XI(3)=X(3) 

C  XI(4)=X(A) 

C  XI(5)=X(5) 

C 

C...  find  the  transmission  error 
C 

C  ERRR=PHI2P-PHI2P0-TN1/TN2*(X(5)-PHI1P0) 

LTKR=PHI2P-PHI2P0+TN1/TN2*(X(5)-PHI1P0) 

ERR (KK) =3600. 0*ERRR/CNST 
PI2P(KK)=PHI2P 
C 

C...  computer  the  contact  path 

C 

xlc=  x2m 

rlc=  dsqrt  (y2m’’'*2+z2in**2) 

xcp(KK)=  xlc*dcos  (rgma2)  •*Tlc*dsin(rgma2)  +ox 
yep  (KK)  "-X 1  c’’'ds  in  (rgnia2)  *rlc*dcos  (rgma2)  ■•■oy 
C 

C...  COMPUTER  THE  PRINCIPAL  DIRECTIONS  AND  CURVATURES  OF  GE.\R 
C 


TH-X(l) 
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PH-X(2) 

ST»DSIN(TH) 

CT»DCOS(TH) 

SH-DSIN(PH) 

CS=DCOS(PH) 

SP=DSIN(ALP2) 

CP=DCOS (ALP2) 

SM='DSIN(G.^2) 

CM=DC0S (GAMA2) 

..  DEFINE  VECTORS  TO  COMPUTER  THE  SECOND  ORDER  PROPERTY  OF  GEAR 

ES(1)=-DSIN(TH-PH) 

ES(2)»  DCOS(TH-PH) 

ES(3)=  0.0 

EQ (1) =-SP*DCOS (TH-PH) 

EQ (2) =-SP*DSIN (TH-PH) 

E0(3)=-C: 

CN (1) =XNM 
CN(2)=YNM 
CN(3)=ZNM 
KS=CP/(RC2-SG*SP) 

KQ=0.0 
W1 (1)=-CM 
Wl(2)=  0.0 
W1  (3)=-SM 
W2(l)=  0.0 
W2(2)=  0.0 
W2(3)=-CR2 

VT1(1)=  yM-'SM-EM2=-SM 
VTl  (2)=-XM*SM-  (Z.M-XB2)*CM 
VTl  (3)=-YM*CM-EM2*C.M 
VT2(1)=  YM-'CR2 
VT2(2)=-XM*CR2 
VT2(3)=  0.0 
DO  110  1=1,3 
W12(I)=W1 (I)-W2(I) 

V12(I)=VT1 (I)-VT2(I) 

110  CONTINUE 
C 

c 

PI21=0.0 

CALL  CURVaI 

K2I=KF 

K2II=KH 

PHI2=PH/CR2 

sh2=dsin  (phi2) 

ch2=dcos  (phi2) 

xX=  C.M*el  (l)+SM*ef  (3) 

yY=  ef(2) 

zZ=-SM*ef  (l)+CM’‘ef  (3) 
ef  (l)=xx 

ef(2)=  CH2-'yY-SH2*zZ 
ef(3)=  SH2*yY+CH2’'zZ 
c 

xX-  CM*eh(l)+SM-'eh(3) 
yY”  eh(2) 

zZ=-SM-'eh(l)+CM*eh(3) 
eh (1)  XX 

eh(2)=  CH2-'yY-SH2’''zZ 
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eh(3)=  SH2*yY+CH2*zZ 

C.  .  . 

CHP=DC0S(PHI2P) 

SHP=DSIN(PHI2P) 

CMM=DCOS (GAMMA) 

SMM=DS IN (GAMMA) 

XX=  ef(l) 

YY=-ef  (2)  *CHP+ef  (3)  -'shp 
ZZ=-ef  (2)  ’''SHP-ef  (3)  *chp 
E2IH(1)=  XX''CMM+22*SMM 
E2IH(2)=  YY 

E2IH(3)— XX*SMM+2Z*CMM 

C.  .  . 

XX>  eh(l) 

YY=-eh (2) *CHP+eh (3) *shp 
ZZ=-eh (2) *SHP-eh (3) *chp 
E2IIH(1)=  XX*CMM+2Z-'SMM 
E2IIH(2)=  YY 
E2IIH(3)=-XX*SMM+Z2”CMM 
C 

C...  COMPUTER  THE  PRINCIPAL  DIRECTIONS  AND  CURVATURES  OF  PINION 
C 

TH1=X(3) 

PH1=X(4) 

STP=DSIN(TH1+PH1) 

CTP=DC0S(TH1+PH1) 

IF(JCC.EQ.l)  THEN 
SP1=DSIN(ALP1) 

CP1=DCOS(aLP1) 

ELSE 

SGN=ALP1 /DABS (ALP  1) 

ALP=SGN*ALP 

SP1=DSIN(ALP) 

CPl=DCOS(ALP) 

END  IF 

SM1=DSIN(GAMA1) 

CMl=DCOS (GAMAl) 

C 

C...  DEFINE  VECTORS  TO  COMPUTER  THE  SECOND  ORDER  PROPERTY  OF  PINION 
C 

ES(1)=-STP 

ES(2)=  CTP 

ES(3)=  0.0 

EQ(1)=  SPP''CTP 

EQ(2)=  SP1*STP 

E0(3)=-CP1 

CN(1)=XNM1 

CN(2)=YNM1 

CN(3)=ZNM1 

IF  (JCC.EQ.l)  THEN 

KS=CP1/(RCF>SF*SP1) 

KQ=0.0 

ELSE 

KS=DCOS (ALP) / (RHO*DCOS (ALP) +XO) 

KQ“1.0/RHO 
END  IF 
W1  (D-  CMl 
Wl(2)=  0.0 
Wl(3)=  SMI 
W2(l)=  0.0 
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W2(2)-  0.0 

W2(3)-  CRIT 

VTl  (D— YM1*SM1-EM1*SMI 
VT1(2)-  XM1*SM1-(ZM1-XB1)*CM1 
VTl  (3)-  YM1*CM1+EM1*CM1 
VT2(1)— YM1*CR1T 
VT2(2)-  XMI^CRIT 
VT2(3)-  0.0 
DO  210  1-1,3 
W12(I)-W1  (I)-W2(I) 

V12(I)-VT1 (I)-VT2(I) 

210  CONTINUE 
C 
C 

PI21-PCR1T 
CALL  CURVaI 

C  WRITE(9,12)  KF.KH.SIGSF 

KlI-KF 
KlII-KH 

C  PHIl-PHl/CRl 

SH1=DSIN(PHI1) 

CHl-DCOS(PHIl) 

XX-  CM1*EF(1)+SM1*EF(3) 
yY-  ef(2) 

Z2— SMl^EF  ( 1) +CMP''EF  (3) 
ef  ( 1) =xx 

EF(2)=  CHr''YY+SHl*ZZ 
EF(3)=-SHl-'YY-rCHl"ZZ 
c 

XX- 

yY-  eh (2) 

Z2--SM1“EH(1)+CM1*EH(3) 
eh  ( 1) -XX 

EH(2)=  CHP'YY+SHl-'ZZ 
EH(3)=-SH1*YY+CH1*Z2 

C.  .  . 

CH1P=DC0S(X(5)) 

SH1P=DS1N(X(5)) 

E1IH(1)=EF(1) 

E1IH(2)=  CH1P*EF(2)-SH1P*EF(3) 

E1IH(3)=  SH1P*EF(2)+CH1P’''EF(3) 

E1IIH(1)-EH(1) 

E1IIH(2)-  CH1P*EH(2)-SH1P*EH(3) 

E1IIH(3)-  SH1P*EH(2)'^CH1P*EH(3) 

DO  109  1-1,3 
E1IH(I)=-E1IH(I) 

E1IIH(I)--E1IIH(I) 

109  CONTINUE 

COMPUTER  THE  DIMENSION  AND  ORIENTATION  OF  THE  CONTACT  ELLIPSE 

GNH(1)=XNH2 
GNH(2)=yNH2 
GNH(3)=2NH2 
CALL  ELLIP 
AXl (KK)-A2L 
AX2(KK)=B2L 
ANGl (KK)=TAU1R 
ANG2(KK)=TAU2R 
C 
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KK-KK+1 

PHI2P-PHI2P+180.0*CNST/(TN2*6.0) 

IF(PHI2P.LE.  (PHI2P2+0.0001))  GOTO  333 

C.  . 

C.  .  . 

WRITE(9,441) 

441  FORMAT  (/  ,  '  ****■'***  ••****-^*-''**='***-'’'*-'**:-******;f-r*-*Vr«*->**;c*'it;r  ' 

&  /.'"TRANSMISSION  ERROR  IN  A  MESHING  PERIOD 

C 

DO  444  I=1,KK-1 
PI2P(I)=PI2P(I)/CNST 
WRITE(9,555)  PI2P(I) .ERR(I) 

555  F0RMAT(3X,3(G14.7,3X)) 

444  CONTINUE 
C 

WRITE(9,551) 

551  FORMAT  (/,  '  _ 

i,  /,'■•  CONTACT  PATH  FOR  A  PAIR  OF  TEETH  IN  MESH 

DO  666  I=1,KK-1 
WRITE(9,747)  XCP(I) . YCPII) 

747  F0RMAT(3X,2(G14.7,3X) ) 

666  CONTINUE 
C 

WRITE(9,661) 

it  / ,  ' DIMENSION  .\ND  ORIENTATION  OF  CONTACT  E'L^°S' 

DO  888  I=i.KK-l 

WRITE  (9, 889)  .\X1  i  li  ,.ANGl  '\)  .P-Xl^l)  ..ANG2(:) 

889  F0RMAT(3X.4(G14.7.3X) I 

888  CONTINUE 
C 

c 

c 

IFCCL.EQ.  I'j  GOTO  1111 
IF(JCL.EQ.3j  goto  1113 
c 

C...  V  and  H  check  for  toe  POSITION 

c 

HMT  =  WD-CC-3 . 0/4  .  O'TW'^  (DT.AN  (Fa)  -DT.AN  (RA)  ) 
DED2T=DED2-3.0/4.0‘-'FW*DTAN(Ra) 

TMCD=MCD-0. 25-TW 

XL=TMCD^^DCOS  (  PGM,a2)  *  (DED2T-HMT/2 . 0)  ^-'DS  IN  f  PGM.A2 ) 
RL=TMCD*DSIN(PG>Ua2)-  (DED2T-HMT/2.0)  ’'DCOS  (  PGM.A2) 

C 

C...  FIND  THE  MEAN  CONTACT  POINT  ON  THE  GEAR  SURFACE 
C 

ERRREL=0. ID-7 
N=2 

ITMAX=200 

IF  (JCH.EQ.l)  THEN 
XI  (1)=270.0"CNST-B2 
ELSE 

XI  (I)=B2 

C  XI  (1)=90.0’--CNST-B2 
END  IF 
XI  (2)=0.0 

CALL  DNEQ.NF  ( FCN 1  .  ERRR  EL ,  N ,  I TMAX ,  X I  .  X  .  FNORM; 
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TH=X(1) 

PH-X(2) 

ZYl-X(l) 

ZY2-X(2) 

N=3 

ERRREL=0. lD-10 
ITMAX=200 
XI (1)=0.0 
XI  (2)=THF 
XI(3)=  0.0 
IF  (JCC.EQ.l)  THEN 

CALL  DNEQNF  (FCNM ,  ERRREL ,  N ,  IT.MAX . XI , X .  FNORM) 

ELSE 

CALL  DNEQNF  (FCNMR.  ERRREL,  N.  IT.MAX.  XI  .X,  FNORM) 

END  IF 
PHI2P0=X(1) 

XI  (1)=ZY1 
XI f2)=2Y2 
XI  (3)=X(2) 

XI  (4)=X(3) 

XI (5)=PHI1P 
WRITE(9, 1A9) 

149  FORMAT  (//6X.  '  ******Y:Vr*;i,r::**::Vr*::*, :*****::.::;, r::::  •  _  / 
i.  6X.'"  V  and  H  CHECK  AT  TOE  POSITION  ”  '  ,  / 

WRITE(9,139)  V,H 

139  F0RM.4T(//4X.  '  V  =  '  .G14. 7,  • H  ='.014.7/  \ 

JCL=3 

GO  TO  5555 
C 

C...  V  .and  H  CHECK  FOR  HEEL  POSITION 
C 

Cl  1 13  HMH=WD-CCC-1 . 0/4. 0^'FW*  (DT.AN  t  FA^  -DT.ANA.Ra: 

C  DED2H=DED2-  1 . 0/4. 0*FW=-DTaN  !Ra) 

C  HMCO=MCD-0.25-'FW 

1113  HMH=WD-CC-0. 16-TW-MDTAN';FA/*DT.AN'Ra) 

DED2H=DED2-0.  16''FW'-DT.aN(RA) 

HMCD=MCD-0. 16'TW 

XL=HMCD'-'DC0S  (PGM.A2)  *  (DED2H-H.MH/2 . 0)  ■■DS  IN  i  PGM.^2) 

RL=HMCD"DSIN  (PGMA2)  -  (DED2H-H.MH/2. 0)  "DCCS  (  PG.M-^Z  ; 

ERRREL=0. lD-7 
N  =  2 

ITMAX=200 

IF  (JCH.EQ.l)  THEN 
XI  (1)=270.0--'CNST*B2 
ELSE 

C  XI  (1)=90.0--CNST-B2 
XI (1)=B2 
END  IF 
XI C2)=0.0 

CALL  DNEQNF(FCN1, ERRREL, N,  IT.MAX, XI, X.FNORM) 

TH=X(1) 

PH=X(2) 

ZY1=X(1) 

ZY2=X(2) 

C...  FIND  THE  V  AND  H  VALL’E  FOR  HEEL  POSITION 
N=3 

ERRREL=0. lD-10 
ITMAX=200 


160 


XI(l)-0.00 
XI (2)»THF 
XI(3)=  0.0 
C  XI(2)»THF+0.2 

C  XI(3)=-0.2 

IF  (JCC.EQ.l)  THEN 

CALL  DNEQNF(FCNM,ERRREL.NMTHAX.XI,X.FN0RM) 

ELSE 

CALL  DNEQNF(FCNMR.ERRREL,N, ITMAX.XI .X.rNCRM/ 

END  IF 
PHI2P0=X(1) 

XI  (1)=ZY1 
XI  (2)=ZY2 
XI  (3)=XC2) 

XI  (4)=X(3) 

XI  ('5)=PHI1P 

C  WRITEf9,ll)  PHI2P0.PHI1P 
WRITE'9, 159) 

159  F0RHaT(//6X.  ■ 

t,  6X.  ' V  and  H  CHECK  aT  HEEL  PCS  IT  1  - 

i  6X ,  ’  "  ' ''  "  "  '' . . 

WRITE{9,169)  V,H 

169  FOR.HAT  (//4X,  '  V  =  ■  ,G14.7,  ' H 

c. .  . 

JCL=1 
GOTO  5555 
nil  CONTINUE 

I F  '  KS  IDE .  Ev  •  0 '  GOTO  ipy',' 

STOP 
END 
C 

C...  FCNl  IS  TO  FIND  THE  ME.AN  CCNTACT  POINT 
C 

SUBROUTINE  FCNllX.F.N' 

IMPLICIT  REAL^-'S  (.^-H.Q-Z) 

INTEGER  N 

REALMS  X(N) , FIN; .mcd 

COMMON/Al  /CNST  .  TN 1 .  TN2  .  C  ,  FW .  G.AMMa  .  x  1  .  r  1  .  rr.c c 
COMMON/A3/B2  ,  RG.M.A2  .  FG.M.A2  ,  ?CM.a2  .  D2?. .  D2F .  aDD2  .  D 
COMMON /a4/ SR  2 . 02  .  RC2  .  PW2  .  X32  .  XG2  .  EM2  ,  Ga.M.A2  .  CF 
C0MM0N/A5 / SG , XM . YM . ZM . XNM , YNM . ZNM . X2M . Y  2M . 2  2M . XN 
i,XNH2  ,  YNH2  ,  ZNH2  .  XH2  .  YH2 . 2H2 
TH=X(1) 

PH=X(2) 

SP=DSIN(ALP2) 

CP=DC0S (ALP2) 

SM=DSIN(GAMA2) 

CM=DC0S(G.AMA2) 

STP=DSIN(TH'PH) 

CTP=DCOS(TH-PH) 

XNM=-CP*CTP 
YNM=-CP"STP 
ZNM=  SP 

A.A1=RC2*STP'SR2''DSINi-02-PH) 

Aa2=RC2-'CTP*SR2-DCOS  ' -Q2-PH'. 

AX=-EM2-'SM 
AY=  <B2'-'CM 
AZ=  E:-i2-CM 

c 

C.  .  .  FIND  SG 


u  V  \  'I  M  7 2 
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to  a 


oonon  non 


C 

Tl-  XNM* (AX-AAl* (SM-CR2) ) *YNM* (aY+AA2* (SM-CR2) ) -ZNM* (aZ+AA1  -XM) 
T2=-XNM*  (SM-CR2)  *SP*STP+ YNM*  (  (SM-CR2)  ’''SP’'CTP-CP*CM)  +ZNM*CM*SP*STP 
SG=T1/T2 

XM“  (RC2-SG*SP)  *CTP+SR2-'DC0S  (-Q2-PH) 

YM=  (RC2-SG'^SP) ''STP+SR2*LSIN(-Q2-PH) 

ZM=-SG*CP 
XM=-SG*SP*CTP^AA2 
YM=-SG-'SP*STP-^AA1 
ZM=-SG*CP 

xX=  CM-'XM-^SM*ZM-XG2-XB2"SM 
yY=  YM-t-EM2 

z2'-SM*XM+CM*ZM-XB2*CM 
XN-  CM*XNM-SM"ZNM 
YN=  YNM 

ZN  =  -SM’'XNM->-CM"ZMM 
PHI2=PH/CR2 
sh2=ds in (phi2) 
ch2=dcos  (,phi2) 

X2M=  xX 

Y2M=  CH2"yY-SH2*2Z 
Z2M=  SH2"yY-CH2*zZ 
XN2M=  XN 

YN2M=  CH2*YN-SH2«ZN 
ZN2M=  SH2-'YN*CH2*ZN 
F(1)=X2M-XL 

F  (2)  =  Yy-”'2-Zz""2-RL-”'2 
F  (2)  ='Y2M**2*22M*''2-RL**2 
RETURN 
END 

...  SUBROUTINE  CURVAl  IS  TO  COMPUTER  THE  CURVATURE  OF  THE 
GENERATED  SURFAFE 

SUBROUTINE  CURVaI 
IMPLICIT  REaL-8(a-H,0-Z) 

REAL-'8  KS,KQ,KF,KH 

DIMENSION  ESN(3) ,EQN(3) .W1VT2(3) ,WV12(3) ,W2VT1 (3) 

C0MM0N/a6/ES(3) ,EQ(3) ,CN(3) ,W1 (3) .W2(3) .W12(3) , VTl (3)  VT2(3) 
SV12(3) .KS,KQ,KF,KH,EF(3) .EH(3) ,SIGSF.PI21 

C.  .  . 

ESN  ( 1 )  =  CN  (2)  ■■■ES  (3)  -CN  (3)  *ES  (2) 

ESN  a)  =-  (CN  (1)  *ES  (3)  -CN  f  3)  •■ES  (1) ) 

ESN  (3) =  CN  ( 1 ) *ES (2) -CN  (2) *ES ( 1 ) 

C.  .  . 

EQN(1)=  CN(2)  •'EQ(3)-CN(3) -E0(2) 

EQN  (2)  =- (CN  (1 )  *EQ  (3) -CN  (3) -'EQ  (1) ) 

EQN  ( 3 )  =  CN  ( 1 )  * EQ  ( 2 )  -CN  ( 2)  ='EQ  ( 1 ) 

C.  .  . 

W1 VT2  ( 1 )  =  W1  (  2)  *VT2  (3)  -W1  (3)  ■'VT2  (2) 

W1VT2(2)=-(W1  (1)’''VT2(3)-W1(3)*VT2(1)) 

W1VT2(3)-  W1 (1)*VT2(2)-W1 (2)*VT2(1) 

C.  .  . 

W2VT1(1)=  W2(2)-'VT1  (3) -W2  (3)  "VTl  (2) 

W2VT1 (2) =- (W2 (1) *VT1  (3)-W2 (3) *VT1  (1) ) 

W2VT1(3)=  W2(1)*VT1 (2) -W2 (2) *VT 1 ( 1) 

C.  .  . 

WV12(1)=  W12(2)  •'V12(3)-W12(3)*V12(2) 

WV12(2)=- (W12(l) *V12(3)-W12(3)*V12(1)) 

WV 1 2  (3)  =  W 1 2  ( 1 )  -' V 1 2  (2)  -W1 2  (2)  " V 1 2  ( 1 ) 
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non  o  n  o 


V12S-0.0 
V12Q-0.0 
WNES-0.0 
WNEQ'0.0 
VWN-  0.0 
WlTN-0.0 
W2TN='0.0 
VT2N-0.0 

DO  1  1-1,3 

V12S-  V12(I)*ES(I)+V12S 
V12Q-  V12(I)*EQ(I)+V12Q 
WNES-  W12(I)*ESN(I)*WNES 
WNEQ*  W12(I)*EQN(I)^WNEQ 
VWN  =  CN(I)*WV12(I)^VWN 
W1TN-  CN(I)-'W1VT2(I)*W1TN 
W2TN-  CN(I)-'W2VT1  (I)-W2T.N 
VT2N=  CN(I)*VT2(I)-^VT2N 

1  CONTINUE 

COMPUTER  THE  CURVATURE  OF  THE  GENERATED  SURFACE 

A13— KS*V12S-WNES 
A23=-KQ"V12Q-WNEQ 

A33=KS''V12S"*2*KQ"V12Q''"2-VWN-W1TN*W27\+P:2r-VT2N/W2  (3) 
Tl=2.0D00"Al3*A23 
T2-A23'”-2-Al3**2- (KS-KQ) *a23 
SIG1F=0.5D00*DATAN2(T1 ,T2) 

KF-0. 50D00" (KS*KQ) -0. 5D00- fAl 3*^2 -a23  •^'2) /a33 
&-^A13'--a23/  (a33"DSIN(2.0D00-SIG1F}  ) 

KH=  KF-2.0D00^M3^-'a23/  (a33’-'DSI.\  :2.0D00=--srGlF;  ) 

SIGSF-SIGIF 
DO  2  1=1,3 

EF(n=  DC0S(SIG1F)*ES(I)-DSIN(.SIG1F)"EQ(I) 

EH(I)=  DSIN(SIG1F)-"ES(  n-DCOS(SIGlF''  '••EQC!) 

2  CONTINUE 
RETURN 
END 

..  FCN2  IS  TO  FIND  THE  INITIAL  GEAR  ROTATIONAL  ANGLE 

SUBROUTINE  FCN2(X,F,N) 

IMPLICIT  REAL*P  (a-H.O-Z) 

INTEGER  N 
REAL--8  X(N),F(N) 

COMMON/Al /CNST ,  TN 1  ,  TN2  , C  ,  FW ,  GAMM.^,  XL ,  RL ,  MCD 

COMMON /a5/SG . XM . YM , ZM , XNM , YNM , ZNM , X2M . Y2M , Z2M , XN2M , YN2M . ZN2M , 
i.XNH2  ,  YNH2  ,  ZNH2  .  XH2  ,  YH2  ,  ZH2 
CM=DCOS (gamma) 

SM=DS IN (GAMMA) 

CHP=DCOS(X(l) ) 

SHP=DSIN(X(1)  ) 

XX=  X2M 

YY=-Y2.M-XHP*Z2M*SHP 
ZZ=-Y2M*SHP-Z2M*CHP 
XH2=  XX*CM-^ZZ*SM 
YH2=  YY>-C 
ZH2=-XX-SM+ZZ*CM 
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o  o  o 


XX-  XN2M 

YY— YN2M*CHP+ZN2M*SHP 
ZZ— YN2M*SHP-ZN2M*CHP 
XNH2-  XX*CM+ZZ*SM 
YNH2=  YY 

ZNH2— XX’^SM+ZZ^CM 

C.  .  . 

R12-TN1/TN2 

V12X=-(YH2-C)*SM*R12 

V12Y-  XH2*SM*R12  +  (1 . 0+R12=-CM)  *ZH2 

V 1 2Z— YH2*  ( 1 . 0+R 1 2*CM)  -^C*CM*R  1 2 

F(1)=XNH2*V12X+YNH2*V12Y+ZNH2*V122 

RETURN 

END 

C.  .  . 

C...  THE  FOLLOWING  IS  THE  TCA  SUBROUTINE  FOR  CURVED  BL.ADE 

C.  .  . 

SUBROUTINE  FCNR(X,F,N) 

IMPLICIT  REAL'"8(a-H,0-Z) 
real "8  x (N) , f (N) 

DIMENSION  CH(3) ,P(3) ,E1EF(3) ,ESN(3) .EQN(3) ,W1VT2(3} ,WV12i3) , 
$W2VT1  (3)  ,EFIH(3)  ,EFIIH(3)  .RH(3)  .GNH(3)  ,E2IH(3)  ,E2IIH(,31  . 
i.ElIH(3)  ,E1IIH(3)  ,EFI(3)  ,EFII(3)  ,  El  I  (3)  .  El  1 1  (3)  .GN(3)  ,EFE1  ^3  )  . 
£.ERR(20)  ,XP(20)  ,YP(20) 

COMMON/Al/CNST, TNI ,TN2,C,FW, Gamma, xl, rl ,mcd 
C0MM0N/A2/B1  ,RGMAl,FGM,Al,PGMAl  ,D1R.D1F,.ADD1  .DEDl 
COMMON/A3/B2  ,  RGMA2  ,  FGM.A2 ,  PGMA2  ,  D2R  .  D2F ,  ADD2  .  DED2  .  WD ,  CC  .  D2P 
COMMON/AA/SR2 , Q2 , RC2 , PW2 , XB2 , XG2 , EM2 . GaMA2 , CR2 . AL?2 , PH  1 2 . PHI  2? 
COMMON/ a5/SG,XM,YM,2M,XNM,YNM,2NM.X2M.Y2M.Z2M,XN2M.YN2M,ZN2M, 
S<XNH2 ,  YNH2 ,  ZNH2 ,  XH2 ,  YH2 .  ZH2 

COMMON/a6/ES(3) ,EQ(3) ,CN(3),W1(3) ,W2(3) .W12(3) ,VT1-3' ,VT2i3) . 
$V12(3) ,KS,KQ,KF.KH,EF(3} ,EH(3) ,SIGSF,PI21 
COMMON/A7/SRl,Ql.Rcf ,PWl,XBl.XGl.EMl,GaM.Al .CRI.ALP: .PHIl .PHIIP 
COMMON/A8/Sf ,XM1, YMl.ZMl.XNMl, YNMl.ZNMl.XlM. YIM.ZIM, 

&.XN1M,  YN1M,ZN1M,XNH1 ,  YNHl  ,ZNHl  ,XH1 ,  YHl  .ZHl 
COMMON/a9/PHI2P0.OX,O2,XO,ZO.RHO.ALP, V.H.CRIT, PCRIT 
COMMON/a11/RAM,PSI1.C2.D6,E2A,F12O,CX6,DX2A,EX120.RU1 ,DELT.RUP, 
SRAl,CPF,DPF,EPF.FPF 
TH=X(1) 

PH=X(2) 

SP=DSIN(ALP2) 

CP=DCOS(ALP2) 

SM=DSIN(GAMA2) 

CM=DC0S(GAMA2) 

STP=DSIN(TH-PH) 

ctp=dcos(th-ph) 

XNM=-CP*CTP 
YNM=-CP*STP 
ZNM=  SP 

AAl=RC2''STP-t'SR2-DSIN(-Q2-P’''' 

AA2=RC2*CTP+-SR2*DC0S  (-Q2-PH) 

AX=-EM2*SM 
AY=  XB2*CM 
AZ=  EM2*CM 

. . .  FIND  SG 

Tl=  XNM*(AX-AAl*  (SM-CR2))*YNM-  (AY-AA2'MSM-CR2 )  ;  -ZNM- (AZ-.^Al  ^■'CM) 
T2=-XNM*(SM-CR2)  ■■SP*STP-YNM*  (  (SM-CR2  )  •■SP-^TP-CP^'CM  ’  -ZNM  ■'CM^'SP -STP 
SG-T1/T2 
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XM=  (RC2-SG*SP)  •'CTP+SR2*DC0S  (-Q2-PH) 

YM=  (RC2-SG’''SP)*STP-^SR2*DSIN(-Q2-PH) 

ZM“-SG’'CP 

C  Xrt=-SG*SP*CTP+AA2 

C  YM=-SG*SP*STP+AAl 

C  ZM=-SG’'CP 

xX=  CM*XM+SM'“ZM-XG2-XB2*SM 
yY=  YM+EM2 

zZ=-SM*XM+CM*ZM-XB2-'CM 
XN=  CM*XNM+SM*ZNM 
YN=  YNM 

ZN=-SM*XNM+CM*ZNM 
PHI2=PH/CR2 
sh2=ds in  (phi2) 
ch2=dcos  (phi2) 

X2M=  xX 

Y2M=  CH2"yY-SH2"z2 
Z2M=  SH2*yY-^CH2"z2 
XN2M=  XN 

YN2M=  CH2*YN-SH2*2N 
ZM2M=  SH2*YN+CH2*ZN 
CMM=DCOS (GaMKA) 

SMM=DS IN (Gamma) 

CHP=DC0S(PHI2P) 

SHP=DSIN(PHI2PJ 
XX=  X2M 

YY  =  -Y2M*CHP--Z2M*SHP 
ZZ=-Y2M*SHP-22M"CH? 

XH2=  XX’-CMM^ZZ'-SMM 
YH2=  YY*C-V 
2H2=-XX‘'SMM+ZZ''CMM 

C.  .  . 

XX=  XN2M 

yy=-yn2m''C::p-"2N2M"Shp 
ZZ  =  -YN2M"SHr-ZN2M-'CHP 
XNH2=  XX''CMM*ZZ"SMM 
YNH2=  YY 

ZNH2  =  -XX"SMM-^ZZ’XMM 

C...  DEFINE  THE  PINION  SURFACE 
C 

TH1=X(3) 

PH1=X(A) 

SM1=DSIN(GAM.'\1) 

CMl=DCOS (GAMAIJ 
STP=DSIN(TH1-"PH1) 

CTP=DC0S(TH1^PH1) 

C 

C...  FIND  CR1T,PF,PPF,PCRIT 
C 

DDD=DABS (PHI) 

IF (ODD. LE. 0.001)  GOTO  6 

PHIl  =  RAl"(PHl-CPF*PHl**2-DPF*PHl**3-EPF-THl*’-4-FPF*PHr--'5) 
PF=RAl’-  (1 .0-2.0*CPF’=PHl-3.0*DPF*PHl*--2 
$-4 . 0*EPF*PHr-*3-5 . 0"FPF*PH1**4) 

PPF=-RAr- (2.0’'CPF+6.0*DPF*PHl*12.0"EPF*PHr-*2*20.0*FPF"PHr-”3) 

CR1T=1.0/PF 

PCR1T--PPF/PF''"3 

GOTO  7 

6  PHIl=RAr-PHl 
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r>  o  o  o  o 


CRIT-CRI 

PCR1T-2.0*CPF/(Ra1**2) 

7  CONTINUE 
CR1T=CR1 
PCRl 7=0.000 

FIND  THE  NOMAL  OF  THE  EQUIDISTANCE  SURFACE 

XMO-  XO*CTP+SRl*DCOS (-Q1-PH1) 

YMO=  X0*STP+SR1*DSIN(-Q1+PH1) 

ZMO=  ZO 

V1X=-YM0*SMI-EM1*SM1 
VI Y-  XM0*SM1-(ZM0-XB1)*CM1 
VIZ-  YM0>''CM1  +  EM1*CM1 
V2X— YM0*CR1T 
V2Y=  XMO^CRIT 
V2Z=  0.0 
VX-V1X-V2X 
VY-V1Y-V2Y 
VZ-V1Z-V2Z 
TX— CTP 
TY— STP 
TZ-0.0 
FX-  STP 
FY— CTP 
FZ-0.0 

XNN=  FY’'VZ-FZ*VX 
YNN=  F2*VX-FX*VZ 
ZNN=  FX-'VY-FY*VX 
DDD=DSQRTi;XNN**2*YNN**2-2NN**2) 
XNMl-XNN/DDD 
YNMl-YNN/DDD 
ZNMl-ZNN/DDD 

DT=TX*XNM1  +  TY-'YNM1  +  TZ*ZNM1 
IF(DT.GE.O.O)  GOTO  10 
XNM1--XNM1 
YNMl  —  YNMl 
ZNMl—ZNMl 
10  CONTINUE 

XM1=  XMO-RHO*XNMl 
YM1=  YMO-RHO*YNMl 
ZM1=  ZMO-RHO*ZNMl 

ALP=DARCOS  (TX*XNM1  +  TY*YNM1-^TZ*ZNM1) 
xX=  CM 1 *XM 1 +SM 1 *ZM 1 -XG 1 -XB 1* SM 1 
yY-  YMl+EMl 

zZ— SM1*XM1+CM1*ZM1-XB1'-CM1 

XN1=CM1*XNM1+SM1*ZNM1 

YNl-YNMl 

ZNl  — SMr''XNMl+CMl*ZNMl 
C  PHIl-PHl/CRl 
shl=ds in  (phi  1) 
chl-dcos (phi  1) 

XIM-  xX 

YIM-  CHl*yY+SHl’''zZ 
ZlM— SHl*yY+-CHl*zZ 
XNIM-  XNl 

YNIM-  CHl*YNl+SHr’'ZNl 
ZNIM— SH1*YN1+CH1*ZN1 
PHI1P-X(5) 
shlP”dsin(phil’’) 
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chlP=dcos  (phi IPj 
XH1=  XlM+H 

YH1=  CH1P’'Y1M-SH1P*Z1M 
2H1=  SH1P*Y1M+CH1P*Z1M 
XNH1=  XNIM 

YNH1=  CH1P’'^YN1M-SH1P*Z.VIM 

ZNH1=  SH1P*YN1M+CH1P*ZN]M 

F(1)=XH2-XH1 

F(2)=YH2-YH1 

F(3)=ZH2-ZH1 

F(A)=XNH2-XNH1 

F(5)=ZNH2-ZNH1 

RETURN 

END 

C.  .  . 

C...  THE  FOLLOWING  IS  THE  SUBROUTINE  FOR  STRAIGHT  BLADE 

C.  .  . 

SUBROUTINE  FCN(X,F,N) 

IMPLICIT  REAL'-S  (a-H.O-Z; 
realms  x(N) , f (N) 

DIMENSION  CH(3) , P  (3) , ElEF  (3) ,ESN(3) ,EON(3) ,W1VT2(3) ,WV12(3)  , 
$W2VT1(3) ,EFIH(3) ,EFIIH(3) ,RH(3) ,GNH(3) ,£2IH(3) ,E2IIH(3) , 
i.ElIH(3)  ,E1IIH(3)  .EFI(3)  .EFII(3)  ,E1I  (3)  ,  El  1 1  (3)  ,GN(3)  ,EFE1  (3  )  , 
iERR(20) ,XP(20) ,YP(20) 

COMMON/a1/CNST,TN1 ,TN2,C,FW, GAMMA, xl,rl .mcd 
COMMON/A2/B1 ,RGMAl , FGMAl , PGMAl , DIR . D1 F . ADDl , DEDl 
COMMON/ A3/B2  ,  RGM.^2  ,  FGMA2  .  PGMA2 ,  D2R  ,  D2F ,  .ADD2  ,  DED2  .  WD .  CC  .  D2P 
COMMON/AA/SR2,Q2,RC2.PW2,XB2,XG2,EM2.GaMA2,CR2,ALF:.PHI2.?K:2? 
COMMON/aS /SG . XM , YM , ZM , XNM, YNM , 2NM , X2M , Y2M , Z2M , XN2M , YN2M , 2N2M . 
i.XNH2  ,  YNH2 ,  ZNH2 ,  XH2 .  YH2 ,  ZH2 

COMMON/ A6/ES  (3)  ,EQ(3)  .CN(3)  ,Wl(3)  ,W2(3)  ,W12  (3)  ,  VT1(3:)  ,VT2':3;  . 
SV12 (3) .KS.KQ.KF.KH.EF (3; ,EH(3) ,SIGSF,PI21 
COMMON/A7/SR1 ,Q1 .Rc£,PWl,XBl,XGl .EMI .GaMAl ,CR1 ,ALP1 ,PHI 1 .PHI  IP 
COMMO  N7  A8 / S  £ , XM 1 , YM 1 , ZM I , XNM 1 , YNM 1 , ZNM 1 , X 1 M , Y 1 M , Z 1 M , 

6.XN1M,  YNIM.ZNIM.XNHI  .YNHl.ZNHl.XHl ,  YHl.ZHl 
COMMON / a9 /PH I 2P0 , OX , OZ , XC , 20 . RHO , ALP , V , H , CR 1 T , PCR 1 T 
COMMON/aI 1 /RAM , PS  1 1 , C2 , D6 . E2A , FI  20 . CX6 . DX24 , EXl 20 , RU 1 . DELT , RUP 
SRAl ,CPF,DPF,EPF, FPF 
TH=X(1) 

PH=X(2) 

SP=DSIN(aLP2) 

CP=DC0S(ALP2; 

SM=DSIN(GAMA2) 

CM=DC0S (GAMA2) 

STP=DSIN(TH-PH) 

CTP=DC0S(TH-PH) 

XNM=-CP"CTP 
YNM=-CP*STP 
ZNM=  SP 

AAl=RC2*STP-SR2“DSIN(-Q2-PH) 

Aj\2=RC2*CTP+SR2*DCOS(-Q2-PH) 

AX=-EM2*SM 
AY-  XB2*CM 
AZ=  EM2”CM 

.  . .  FIND  SG 


Tl=  XNM*  (AX-AAl*  (SM -CRZ)  )  *YNM*  (AY'-.AA2*  (SM-CR2)  )  -^ZNM*  (AZ-.A.R1 '■CM' 
T2--XNM* (SM-CR2) *SP*STP*YNM* ( (SM-CR2) *SP*CTP-CP*CM) *ZNM*CM*SP^'STP 
SG-T1/T2 
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Xrt=  (RC2-SG’''SP)  ''^CTP*SR2*DC0S  (-Q2-PH) 
YM=  (RC2-SG''SP)*STP^SR2*DSIN(-Q2-PH) 
ZM=-SG-'CP 

C  XM=-SG*SP*CTP+AA2 

C  yM=-SG*SP*STP+AAl 

C  ZM=-SG*CP 

xX=  CiM--XM+SM*ZM-XG2-XB2*SM 
yY=  YM^EM2 

zZ=-SM--XM+CM-'ZM-XB2"CM 
XN=  CM'-XNM+SM*ZNM 
YN=  YNM 

ZN=-SM*XNM+CM*ZNM 
PHI2=PH/CR2 
sh2=ds in (phi  2) 
ch2=dcos (phi  2) 

X2M=  xX 

Y2M=  CH2-'yY-SH2*2Z 
Z2M=  SH2*yY^CH2'--2Z 
XN2M=  XN 

YN2M=  CH2-'YN-SH2*ZN 
ZN2M=  SH2*YN>CH2-'ZN 
CMM=DCOS (GAMKa) 

SMN=DS IN (GAMMA) 

CHP=DC0S(PHI2P) 

SHP=DSIN(PHI2P) 

XX=  X2M 

YY=-Y2M"CHP-22M’'SHP 
2Z=';_'r-'SHP-Z2M"CHP 
XH2=  XX’-'CMM-  ZZ’-SMM 
YH2=  YY+C-V 
C  YH2=  YY-C-V 

Z  H  2 = -  XX  '  S  I'lM  -  Z  Z  -  C  MM 

C  .  . 

XX=  XN2M 

YY  =  -YN2M--CHP-^2N2M*SHP 
ZZ=-YN2M-'SHP-ZN2M’-CH? 

XNH2=  XX*CMM*ZZ-'SMM 
YNH2=  YY 

ZNH2  =  -XX-'SMM^ZZ-CMM 
C 

C...  DEFINE  THE  PINION  SURFACE 
C 

TH1=X(3) 

PH1=X(4) 

SP1=DSIN(-ALP1) 

CP1=DC0S(-ALP1) 

SM1=DSIN(GAMA1) 

CMl  =  DCOS(GAMAn 
STP  =  DSIN(TH1-^PH1) 

CTP=DC0S(TH1-^PH1) 

XNMl  =  -CPr'CrP 
YNM1=-CP1''STP 
ZNM1=  SPl 

ABl=RrF"STP*SRl”DSIN(-Ql*PHl) 

AB2=RCF-'CTP*SR1"DC0S(;-Q1+PH1) 

AXX=-EM1*SM1 

AYY=  XB1*CM1 

A2Z=  EMI  ■•■cm; 

C 

C...  FIND  SF,CR1T,PF,PPF,PCR1T 
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c 

DDD=DABS (PHI) 

IFCDDD.LE. 0.001)  GOTO  6 

PHll=RAr''(PHl-CPF*PHl  '•■2-DPF^''PHl**3-EPF*PHr-^4-FPF''PHr"'5  i 
PF=RAl  (1.0-2. 0*CPF"PHl-3 . 0*DPF--'PH1 
$-4.0*EPF-'PHP-*3-5.0^TPF-^PHP'='4) 

PPF=-RAr'(2. 0-'CPF*6.0"DPF*PHl-12.0*i-'’F*PHl-^^2-20.0^TPF'-pHr^^^3) 

rRlT=1.0/PF 

PCR1T=-PPF/PF""3 

GOTO  7 

6  PHI1=RA1*PH1 
CR1T=CR1 

PCRlT=2.n*CPF/  (RAl*^'2J 

7  CONTINUE 

C  CR1T=CR1 

C  PCR1T=0.000 

T 1  =  XNM 1  ( AXX- AB  1  ( SM  i  - CR  1 T )  .)  - 
i.YNMr-  (ayy-^ab2-'(smi-crit)  )  -ZN'Mi'=;.\z2'ABi -cm:  ) 
T2=-XNM1-'(SM1-CRIT)  ^^S?1"STP* 

&.YNM1*  (  (SMl-CKlT)  -'SP  P'CTP-C? ;  "CM  1 )  -ZNM 1 '•'CMl '•^SP  r'STP 
SF=1 ' /T2 
C 

XM1=  (RCF-SF--''S?1)"CT?-SR1"DC0S(-01-PH1) 

YM1=  (RCF-SF"SPl)*ST?-SRr-DSIN(-Ql'PHl) 

ZM1  =  -SF'"CP1 

xX=  CMl-XMl'SMl-ZMl  -XG1-XB1-'-Sm; 
yY=  YMl-EMl 

z2  =  -SMl"-XMl*CMl"ZMl-XBl"CM2 
XN 1 =CM 1 -XNM 1 -SM 1 "ZNM 1 
YN1=YNM1 

ZN 1  =  -  SM  1  -'XNM  1  -  CM :  ■■  Z  NM 1 
C  PHIl=Phl/CRl 

sh l=ds i n ( ph i i ; 
chl  =  dcos  (phi  1 ) 

X1M=  xX 

Y1M=  CHl-yY-SHl'-zZ 
ZlM=-SHl-yY-CHl‘--zZ 
XN1M=  XNl 

YN1M=  CH1-YN'1-SH1"ZN1 
ZN  1  M=  -SH  1  YN 1  -CH 1  '■ZN  1 
C  WRITE (9, 111)  XlM.YlM.ZlM 

PHI1P=  X(5) 
shlP'^ds  in  (phi  IP) 
chlP  =  dcos (phi  IP) 

XH1=  XIM-H 

YH1=  CHlP-'YlM-SHiP-ZlM 
ZH1=  SHlP-'-YlM-^CHlP-ZlM 
XNH1=  XNlM 

YNH 1 =  CH  1 P- YN  1 M- SH  1 P "ZN i M 

ZNH 1  =  SH  1 P  ■■■  YN  1 M  -  CH  1  P  "  ZN  1 M 

F(1)=XH2-XH1 

F(2) =YH2-YH1 

F(3)=ZH2-ZH1 

F(A)=XNH1-XNH1 

F(5)'YNH2-YNH1 

RETURN 

END 


C...  Aa  OOC'IO 

C,  ,  ,  .  00020 

C...  "  SUBROUTINE  ELLI?  IS  TO  DETERMINE  THE  SIZE  AND  OP  I N  T  •.  7 1 0’N  '  AA  00030 
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*  OF  THE  contact  ELLIPSE 


SUBROUTINE  ELLIP 
IMPLICIT  REAL-'8(A-H,0-Z) 

REAL*8  KS,KQ,K1I ,K1II ,K2I.K2II 

DIMENSION  R0(3) ,ETA2(3) ,ZETa2(3) ,E1E2(3) ,ETA(3) ,ZETa(3) 
DIMENSION  E1IH(3) ,E1IIH(3) ,E2IH(3) .E2IIH(3) ,GNH(3) 

COMMON/Al /CNST , TN 1 , TN2 , C , FW , GAMMA , X 1 , r 1 . racd 
COMMON/A3/B2 , RGMA2 , FGMa2 , PGMA2 , D2R , D2F , ADD2 , DED2 , WD , CC , D2P 
COMMON/A4/SR2 , Q2 , RC2 , PW2 , XB2 , XG2 , EM2 , GaMA2 , CR2 , ALP2 , PH  1 2 , PHI2P 
COMMON /a5 / SG  ,  XM  ,  YM  ,  ZM ,  XNM ,  Y.NM .  ZNM ,  X2M ,  Y2M ,  Z 2M ,  XN2M .  YN2M .  ZN'2M , 
5.XNH2 ,  YNH2 ,  ZNH2  ,  XH2 ,  YH2  ,  ZH2 
COMMON /a9 / PH I 2P0 , OX , OZ , XO , ZO , RHO ,ALP,V,H.CR1T,PCR1T 
COMMON/AIO/KII ,K1II,K2I,K2II,DEL,E1IH.E1IIH,E2IH,E2IIH.CNH. 
&a2P,B2P, TAU1R,TAU2R,A2L,B2L 
CNST=DARCOS (-1 -ODOO) /180.00 

C.  .  . 

E1E2(1)=  E1IH(2)  *E2IH(3)-E1IH(3)  ■■E2IHv2) 

E1E2(2)=-  (ElIH(l)  ••E2IH(3;-E1IH(3)  "EZlHll)  ) 

E1E2(3)=  ElIH(l)  *£2IH(2)-E1IH(2)*E2IH(U 

C.  .  . 

T1=0.0 
T2=0.0 
DO  1  1=1,3 

Tl=  E1IH(I)*E2I«(I)-T1 
T2=  GNH(I)-E1E2(I)-T2 
;  CONTINUE 

C.  .  . 

T1=T1‘1.0D00 

S I G 1 2=  2 . ODOO "DAT  AN  2  I T  2 , T 1 ) 

C.  .  . 

SK1=  KII+KIII 
SK2=  K2I+K2II 
SG1=  KII-KIII 
SG2=  K2I-K2II 

C.  .  . 

Tl  =  SGl-SG2--DCOS  (2 . 0D00''SIG12) 

T2  =  SG2*DSIN(2.0D00'''SIG121 

T3=DSQRT(SGl"*2+SG2*"2-2.0D00"SGr-SG2''DC0S ;2. ODOO-SIG12; ) 

C.  .  . 

TX=T2/T3 
TY=T1/T3+1 .ODOO 
ALP12=DAT.AN2(TX,TY) 

C.  .  . 

C...  THE  DIRECTION  AND  LENGTH  OF  THE  .AXES  OF  CONT.ACT  ELLIPSE 

C.  .  . 

r  DFT =0.C0700D00 

AL=0. 25D00-- (SK1-SK2-T3) 

BL=0.25D00*(SK1-SK2*T3) 

C  WRITE(9.5)  SIGI2,AL.BL 

C  5  F0RMAT(2X,  '  SIG12.  AL.BL  =  '  ,  3  2X  .  Gi  4 . 7)  ) 

AL=DABS (AL) 

BL=DABS IBL) 

a2L  =  2.0D00''-DS0RT(DEL/aL) 

B2L=2.0D00"DSQRT (DEL/BL) 

C.  .  . 

C.  .  . 

DO  2  1=1,3 

eta(i)  =  DCC'S (ALPI2 ' 'El  Tr"  : : -Ds :n ; AL?:2^  ^^El  I  IH (  I  i 


AA 

OOOAO 

AA 

00050 

AA 

00060 

AA 

00070 

AA 

00080 

AA 

00090 

AA 

00100 

OOllO 

A--\ 

00120 

AA 

00130 

AA 

00140 

.AA 

00150 

AA 

00160 

AA 

00170 

-AA 

00180 

AA 

00190 

AA 

00200 

002 1 

■AA 

00220 

.AA 

00230 

AA 

00240 

.AA 

00250 

•AA 

00260 

.AA 

00270 

AA 

00280 

.AA 

00290 

AA 

00300 

AA 

00310 

.AA 

00320 

.AA 

00330 

.A. 

00340 

-AA 

00350 

.AA 

00360 

.AA 

00370 

AA 

00380 

AA 

00390 

.AA 

00400 

■AA 

00410 

•AA 

00420 

AA 

00430 

.AA 

o 

O 

o 

00450 

-AA 

00460 

AA 

00470 

AA 

00480 

■AA 

00490 

AA 

00500 

AA 

005  10 
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n  n  o 


ZETA(I)=  DSIN(ALP12)  "El  IH  (I) -^DCOS  (ALP12)  ='ElIIH(i; 

AA 

00520 

2 

CONTINUE 

■AA 

00530 

C. . . 

.AA 

00540 

C. . . 

DETERMINE  THE  PROJECTION  OF  CONTACT  ELLIPS  IN  .AXI.AL  SECTION 

AA 

00550 

C..  . 

AA 

00560 

CHP=DC0S(PHI2P) 

•AA 

00570 

SHP=DSIN(PHI2P) 

.AA 

00580 

C.  . . 

AA 

00590 

CMM=DCOS (GAMMA) 

AA 

00600 

SMM=DS IN (GAMMA) 

.AA 

00610 

c. .  . 

.AA 

00620 

XX=  ETA(1)*CMM-ETA(3)*SMM 

.AA 

00630 

YY=  ETA(2) 

•AA 

00640 

ZZ=  ETA(l) *SMM+ETA(3)*CMM 

A.iA 

00650 

ETA2(1)=  XX 

.AA 

00660 

ETA2 (2) =-YY*CHP-ZZ"SHP 

.AA 

00670 

ETA2(3)=  YY*SHP-2Z“CHP 

-AA 

00680 

c. .  . 

.AA 

00690 

XX=  ZETA(l)  *CMM-ZETA(3) -'SMM 

.AA 

00700 

YY=  ZETA(2) 

-AA 

00710 

ZZ=  ZETA(l)  •■SMM-*-ZETA(3)*CMM 

.AA 

00720 

ZETA2(1)=  XX 

.AA 

00730 

ZETA2 (2) =-YY*CHP-ZZ"SHP 

-AA 

00740 

ZETa2(3)=  YY*SHP-Z2*CHP 

.AA 

00750 

c. .  . 

•AA 

00760 

R  0  ( 2 )  =  Y  2  M  /  D  S  Q  R  T  ( Z  2  M  *  "  2  -  Y  2  M  "  ■■  ■  2 ) 

.AA 

00770 

R0(3)=Z2M/DSQRT(Z2M-"-"2-Y2M**2.; 

.AA 

00780 

R0(1)=0.0D00 

.-v.’k 

c. .  . 

-AA 

00800 

T1 1=0.0D00 

•AA 

00810 

T12=O.ODOO 

.AA 

00820 

DO  3  1=1,3 

.Am. 

00830 

TI2=  ETA2(I)*R0( I)-Ti2 

.-uA 

00840 

T11=ZETA2CI)"R0(I)*T!1 

.AA 

00850 

3 

CONTINUE 

AA 

00860 

c. .  . 

-AA 

00870 

TAUl  =DAT.AN2  f  T1  1  ,  ZETa2  ( i  )  ■ 

.AA 

00880 

TAU2  =  DAT.AN2  (T  1 2  ,  ETA2  i  1 ) 

.AA 

00890 

c. .  . 

.AA 

00900 

a2P=a2L*ZETA2  (  D  /DCOS  (.TALI') 

.AA 

00910 

B2P=B2L*ETA2 (1) /DCOS (TaU2) 

•AA 

00920 

c. .  . 

.AA 

00030 

T.AU1R=  (TAU1-RGMA2)  /CNST 

009^0 

TAU2R=  ('TAU2-RGM.A2)  /CNST 

AA 

00950 

RETURN- 

AA 

00960 

END 

AA 

00970 

c. .  . 

C...  THE  FOLLOWING  IS  THE  V-H  CHECK  PROGR.aM  FOR  CURVED  BL.\DE 

C.  .  . 

SUBROUTINE  FCNMR (X . F , N) 

IMPLICIT  REAL-''8  ( A-H,0-2) 
real  '■'S  x  (N  )  ,  f  (N ) 

COMMON/ Al  /CNST  ,  TN 1  .  TN2  ,  C  .  FW .  G.AMM.A .  x  1  .  r  I  ,  me  d 

COMMON  /a5  /  SG  .  XM  ,  YM  .  ZM .  XNM ,  YNM ,  ZNM ,  X2M .  Y2M .  Z2M  .  XN2M  ,  YN  2M  .  Z'.  2M . 
^XNH2,YNH2,ZNH2,XH2, YH2.ZH2 

COMMON/A7/SR1 .01 , Re f , ?Wl , XB 1 . XGl , EMI .GaMAl .CR 1 , ALP  1 . PH  I  1 . PH  I  1 P 
COMMON/A9/PHI2PO,OX.OZ.XG,ZO.RHO.ALP. V.H.CRIT, PCRIT 
COMMON/ A 1 1 /RAM . PS  I  1 , C2 . D6 . E2d , FI  20 . CX6 , DX2A . EX  1 20 , RU 1 , DELT . RUP , 
SRAl.CPF.DPF.EPF.FPF 
CM  =  DCOS  cga-mma; 
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u  u 


SM-DS IN (gamma) 

C  CHP-DCOS(PHI2PO) 

C  SHP=DSIN(PHI2P0) 

CHP=DCOS (X(l)) 

SHP-DSIN(X(1)) 

XX=  X2M 

YY=— Y2M*CHP+Z2M*SHP 
Z2=-Y2M"SHP-Z2M*CHP 
XH2-  XX*CM+ZZ*SM 
YH2=  YY-^C 
ZH2— XX*SM+ZZ*CM 

C.  .  . 

XX=  XN2M 

YY=-YN2M"CHP+ZN2M"SHP 
Z2=-YN2M-'SHP'ZN2M*CHP 
XNH2=  XX*CM+ZZ*SM 
YNH2=  YY 

ZNH2=-XX*SM^ZZ''CM 

C...  DEFINE  THE  PINION  SURFACE 
C 

TH1=X(2) 

PH1=X(3) 

SM1=DSIN(GAMA1) 

CMl=DCOS(GAMAl) 

STP=DSIN(TH1*PH1) 

CTP=DC0S(TH1*PH1) 

FIND  CR1T,PF,PPF,PCR1T 

DDD=DABS(PH1) 

IF(DDD.LE.0.00n  GOTO  6 

PHI1=Ra1*  (PHl-CPF-THl--'-2-DPF"PHl-'-'3-EPF''PHr--4- 
PF=RAr''  (1 .0-2.0*CPF-PHl-3.0"DPF"PHr"-2 
S  -  A  .  0  *  E  P  F  *  P  H 1  •  •  *  3  -  5 . 0  •  •  F  P  F  -  •  P  H 1  •  •  *  ) 

PPF=-RAr-' (2 . 0''CPF’-6  .  O-'DPF'-PHl- 12. 0'"£PF”?HI  "2*2 
CR1T=1.0/PF 
PCR1T  =  -PPF/PF"''3 
GOTO  7 

6  PHIl  =  RArTHl 
CR1T=CR1 

PCRlT=2.0*CPF/(RAr-'-'2) 

7  CONTINUE 
C 

C  FIND  THE  NOMAL  OF  THE  EQU IDI ST.ANCE  SURFACE 

C 

XMO=  XO'XTP*SRr-DCOS(-Ql*PHn 
YMO=  XO-'STP*SRl*DSIN(-Ql-PHl) 

ZMO=  ZO 

V1X=-YM0"SM1-EM1”SMI 

VI  Y=  XMO-'SMl-  (ZMO-XBl  )  *CM1 

V1Z=  YMO-'CMl^EMl-'CMl 

V2X=-YM0’'CR1T 

V2Y=  XMO-'CRIT 

V2Z=  0.0 

VX=V1X-V2X 

VY-V1Y-V2Y 

V2-V1Z-V2Z 

TX'-CTP 

TY--STP 


FPF”PH  1  ■■''5 ' 

0. 0''FPF"PHI  ““3) 
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TZ=0.0 
FX=  STP 
FY—CTP 
FZ=0.0 

XNN=  FY-'VZ-FZ’''VX 

YNN=  FZ*VX-FX’'VZ 

ZNN=  FX*VY-FY-'VX 

DDD=  DS  QR  T  (XNN  •'  *  2  Y  NN  ■  ■  2  *  ZNN  *  ■'  2  ) 

XNM1»XNN/DDD 

YNM1  =  YNN/l''DD 

2NM1=ZNN/D0D 

DT=TX*XNM1-'  ^Y*YNM1-^TZ*ZNM1 
IF(DT.GE.O.O)  GOTO  10 
XNM1=— XNMl 
YNMl^-YNMl 
ZNM1=-ZNM1 
10  CONTINUE 

XM1=  XM0-RH0*XNM1 
YM1=  YM0-RH0*YNM1 
2M1=  2M0-RH0*ZNM1 

ALP=DARC0S(TX*XNM1^TY’'YNM1-T2*ZNM1) 
xX=  CM  1  -'XN 1  *SM  1  *2M  1  -XG 1  -XB 1  *SM  1 
yY=  YMl+EMl 

zZ = -  SM 1  ’'XM 1  ^  CM  1  -'EM  1  -  XB  1  *CM  1 
XN 1  =CM  r'''XNM  1  SM 1 "  Z  NM 1 
YN1=YNM1 

ZN1  =  -SM1-'XNM1*CM1*ZNM1 
C  PHri=PHl/CRl 

shl  =  Ms in  (phi  1) 
chl'dcos  (phi  1) 

X1M=  xX 

Y1M=  CHP'yY-SHl-zZ 
ZlM=-SHl*yY*CHr-z2 
XN1M=  XNl 

YNIM”  CHr-YNl  +  SHl-'ZNl 
ZNlM=-SHr-YNl*CHl-'ZNl 

C.  .  . 

C 

TT=YN1M**2-ZN1M**2 
SH1P=(-ZN1M*YNH2*YN1M--'ZNH2;’  /tt 
CH1P=  (  YNlM-YNHZ-t-ZNlM^ZNHZ)  /TT 
PHI  1P=2.0D00*DaTA,\2  (SHIP.  (1 .  OdOO*CH1P  ) ) 

C.  .  . 

C.  .  . 

XH1=  XIM 

YH1=  CHIP’-YIM-SHIP-'ZIM 
ZHl-  SHIP'-'YIM-CHIP-'ZIM 
XNH1=  XNIM 

YNH1=  CH1P-YN1M-SH1P'=ZN1M 
ZNH1=  SH1P-YN1M*CH1P''ZN1M 
V  =  -  (YH2-yHl) 

H=XH2-XH1 

C.  .  . 

F(1)=2H2-ZH1 

F(2)=XNH2'XNH1 

C  F(2)=YNH2’”'2*ZNH2*’^2-TT 

C 

C 

C.  .  . 

R12=TN1/TN2 
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V12X=0.0-YH2’'SM’'R12 

V12Y=  ZH1-^R12'MXH2*SM+ZH2*CM) 

V12Z=-YH1-R12*YH2*CM 
C  V12X=-(YH2-(C-V))*SM*R12 

C  V12Y=  XH2-'SM''R12*  (1 .0^R12*CM)  *ZH2 

C  V  12Z»- YH2-'  ( 1 . 0+R  12*CM)  -  (C-V)  ='CM*R  12 
F  ( 3 )  =XNH2 * V 1 2X-*-  YNH2 *V  1 2  Y-ZNH2 -'V 1 2Z 
RETURN 
END 

C.  .  . 

C...  THE  FOLLOWING  IS  THE  V-H  CHECK  SUBROUTINE  FOR  STRAIGHT  BL.ADE 

C.  .  . 

SUBROUTINE  FCNM(X,F.N) 

IMPLICIT  REALMS (a-H,0-Z) 
realms  x (N) , f (N) 

COMMON/Al/CNST,TNl ,TN2.C, FW.GAMMA.xl , rl .mcd 

COMMON/A5 /SG , XM , YM , ZM , XNM. YNM , ZNM . X2M , Y2M , Z2M , XN2M . YN2M . ZN2M , 
&>XNH2 ,  YNH2  ,  ZNH2  ,  XH2  ,  YH2  .  ZH2 

COMMON/A7/SR1  ,Ql,Rcf  .PWl  ,XB1  .XGl  .EMI  .GaM.-kl  .CRl  .ALPl  .PHI  I  ,PKI  1? 
COMMON/a9/PHI2PO.OX.OZ,XO.ZO,RHO.ALP, V.H.CRIT, PCRIT 
COMMON/Al 1/RAM. PSI 1 .C2.D6.E24.F 120.CX6,DX24, EX  120. RUi .DELT.RUP, 
$RAl,CPF,DPF.EPF,FPF 
CM=DCOS (GAMMA) 

SM=DSIN(GAMM/.) 

CHP=DCOS(X(l)) 

SHP=DSIN(X(1)) 

XX=  X2M 

YY  =  -Y2M"CHP-Z2M''SHP 

ZZ=-Y2M*SHP-Z2M'^CHP 

XH2=  XX"CM^ZZ"SM 

YH2=  YY*C 

2H2  =  -XX-''SM-Z2^--CM 

C.  .  . 

XX=  XN2M 

YY  =  -YN2M''CHP-ZN2M"SHP 
ZZ=-YN2M*SHP-ZN2M'THP 
XNH2=  XX*CM-2Z''SM 
YNH2=  YY 

ZNH2  =  -XX''SM*ZZ'^CM 

C.  .  . 

C 

C...  DEFINE  THE  PINION  SURFACE 
C 

TH1=X(2) 

PH1=X(3) 

SP1=DSIN (-ALPI) 

CPl=DCOS (-ALP1) 

SM1=DSIN (GAMAl) 

CMl=DCOS  (GA.M.A1) 

STP-DSIN(THl-PHli 
CTP=DCOS (THl^PHi ) 

XNM1'-CP1‘''CTP 
YNM1  =  -CPP'STP 
ZNM1=  SPl 

AB1=RCF-STP*SR  r-'DSIN  (-01*PH1) 

AB2  =  RCF"CTP*SRP^DC0S  (.-Ql-PHl  ' 

.4XX=-EM1’'SM1 
AYY=  XBl'-'CMl 
AZZ=  EMl’XMl 
C 
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C...  FIND  SF.CRIT.PF.PPF.PCRIT 
C 

PHIl“RAr'(PHl-CPF*PHl’-*2-DPF*PHI**3-EPF-'PHl**4-FPF*PHl**5) 

PF=RA1*(1.0-2.0’-CPF*PH1-3.0*DPF*PH1**2 

$-A.0'-''EPF*PHr''*3-5.0*FPF*PHl**A) 

PPF=-RAr''  (2. 0’''CPF+6. 0*DPF*PH1-12. 0*EPF*PH1  **2*20. 0’'FPF*PH1“’'3J 
CR1T=1.0/PF 
PCR1T=-PPF/PF**3 
C 

Tl=  XNM1*(aXX-.AB1*  (SMI -grit;)* 
i-YNMl*  (AYY+AB2*(SM1-CR1T))*ZNM1*(AZZ*ABP-CM1) 
T2=-XNM1*(SM1-CR1T) -'SPl-'STP* 

&  YNM 1  *  (  (SM 1  -CR 1 T)  *S P 1  *CTP-CP  1  *CM I )  *ZNM  1  ’‘CM  1  *S P 1  *S TP 
SF=T1/T2 

XM1=  (RCF-SF*SP1) *CTP-SR1*DC0S(-Q1*PH1) 

YM1=  (RCF-SF*SP1) -'STP'SRl-'DSINC-Ql  +  PHl) 

ZM1=-SF*CP1 

xX=  CM  1  *XM  1  *  SM 1  •- ZM 1  -XG 1  -XB 1  SM 1 
yY=  YMl+EMl 

zZ=-SMl*XMl*CMl*ZMl-XBi*CMl 

XN1=CM1*XNM1+SM1*ZNM1 

YN1=YNM1 

ZN1=-SM1*XNM1+CM1*ZNM1 
C  PHI1=PH1/CR1 

sh l  =  ds in  (phi  1) 
chl  =  clcos  (phi  1) 

X1M=  xX 

Y1M=  CHl*yY*SHr-'2Z 
ZlM=-SHl*yY-CHl*zZ 
XN1M=  XNl 

YN1M»  CHl*Y.Nl*SHl*2Ni 

ZNlM=-SHr''YNl*CHl*ZNl 

TT*YN1M**2+ZN1M**2 

SH1P= (-ZN1M*YNH2*YN1M*ZNH2) /TT 

CH1P=(  YN1M*YNH2*ZN1M*Z‘-'H2)  /TT 

PHI1P=2.0DOO*DATAN2(SH1P,  ( 1  .ODOO*CHlP'i  ) 

C  SH1P=DSIN(PHIIP) 

C  CH1P=DC0S (PHI  IP) 

XH1=  XIM 

YH1=  CH1P*Y1M-SH1P*Z1M 
ZHl'  SHlP*YlM*CHlP*ZiM 
XNHl-  XNIM 

YNH 1  =  CH 1 P  YN 1 M-SH 1 P  ^^ZM  M 
ZNHl”  SH1P*YN1M*CH1P^^ZNIM 
V=- (YH2-YH1) 

H=XH2-XH1 

C.  .  . 

F(1)=ZH2-ZH1 

F(2)=XNH2-XNH1 

C  F  (2) =YNH2**2*ZNH2**2-TT 

C 

C 

c. .  . 

R12-TN1/TN2 

V12X“0.0-YH2*SM*R12 

V12Y-  ZH1*R12*  (XH2*SM*ZH2-'CM) 

V12Z--YH1-R12*YH2*CM 
C  V12X-- (YH2- (C-V) ) *SM*R12 

C  V12Y=  XH2*SM '^12*  (1 .0*R12^^CM)  *ZH2 

C  V12Z  =  -YH2*(1 .0*R12  -CM) *  iC-Vl*CM"R12 
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F(3)=XNH2*V12X+YNH2*V12Y+ZNH2-'V12Z 

RETURN 

END 

C 

C  SUBROUTINE  CAM  IS  FOR  THE  COEFFICIENTS  OF  GENERATION  MOTION 
C 

SUBROUTINE  CAM  .\A  00070 

IMPLICIT  REALMS  (A-H,  0-2)  .A,A  00080 

COMMON/a11/RAM,PSI1,C2,D6,E24.FI20,CX6.DX2A,EX120.RU1 .DELT.RUP, 

SRAl ,CPF,DPF,EPF,FPF 
T 1  =  1 . 0-^ 3 . 0*C 2 *DTAN  ( P S 1 1 ) 

£.-  ( 1 . 0-RAM)  **3“  (RUl  **3/ 15 . 0*-2-DELT)  /  ( 15 . 0"DC0S  (PS  1 1 ) ) 
T2=1.0+(RUl+DELT)/(15.0"DCOS(PSIi) ) 

CX6=T1/T2 

Tl=  6 . 0’-C2-^DCOS  (PSIi)  -  (A.  0-'CX6-3 . 0--C2-^^2-  1 . 0)  ^'DS  IN  (PS  1 1) 
i,-6 . 0*C2''  ( 1 . 0-RAM)  (RU1**3/15 . 0-^*3-DELT/  15 . 0) 

T2=  DCOS(PSI1)*(RU1-DELT)/15.0 
DX2A=  T1/T2 

T 1=  ( 10 . 0-CX6-  15 . 0-'C2'^'2-  1 . 0)  *DCOS  (PS  1 1 ; 
i.  -  ( 5 . 0  "DX2  A-  1 0 . 0----C2*CX6-  1 0 . 0'''C2 )  ”DS  IN  '  PS  1 1 ) 

-  ( 10. 0*CX6*  ( 1 . 0-RAM)  *--'2 

i,  -  15 . 0'-'C2**2*  ( 1 . 0-R.AM)  )  (RU 1  *"3/15 . 0^'*3-DELr/  15 .0) 
i.  -  (1 . 0-RAM)  •'•*5  {RU  1  ■•■■•'5/  1 5 . 0**5-DELT/  15.0) 

T2  =  DCOS(PSI  D-  (RU1‘DELT)/15.0 

EX120=T1/T2 

D6=CX6-3.0*C2**2 

E2A  =  DX2A-C2’'(  15.0*C2'^--^2-10.0-''CX6! 

rl2G-EXl  20- 1 5 . 0'X2"Da2  A-  1C5.G’C2””2"'CA6-Ci'’2  -iG.'j'CA6''.i 

C?F=C2/2.0 

DPF=O6/6.0 

E?F=E2A/2A.O 

F?F=F120/120.0 

RETURN 

END 
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